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SUMMARY 


An  acousto-electromagnetic  land  mine  detection  technique  is  being  investigated. 
A  two-dimensional  finite-difference  model  for  elastic  waves  has  also  been  developed,  but 
it  is  a  purely  linear  model.  Strong  nonlinearities  are  typical  of  the  soils  in  which  mines 
are  buried.  The  purpose  of  this  thesis  is  to  characterize  these  nonlinearities  for  the 
propagation  of  high  frequency  seismic  waves  (30  -  2000  Hz)  in  moist,  compacted  sand 
so  that  the  parameters  used  in  acousto-electromagnetic  land  mine  detection  may  be 
improved  and  the  nonlinearities  may  be  incorporated  in  the  computer  model. 

The  frequency  response  of  the  soil  model  was  recorded  as  a  function  of  drive 
amplitude  and  propagation  distance.  The  amplitude  response  of  the  soil  model  was 
recorded  as  a  function  of  frequency  and  propagation  distance.  The  fundamental  and  first 
five  harmonics  were  saved  for  each.  Three  elastic  wave  transducers  (shakers)  were 
characterized  so  that  source  nonlinearities  could  be  compared  to  propagation  path 
nonlinearities.  Characterization  of  the  shakers  included  foot  motion  under  unloaded  and 
sand-loaded  conditions. 

The  source  and  propagation  path  produced  nonlinearities  as  shown  by  harmonic 
generation  in  accelerometers  mounted  to  the  shaker  foot  and  radar  measurements  of  the 
soil  surface  displacement.  Frequencies  in  the  100  -  600  Hz  band  propagated  best  while 
frequencies  above  600  Hz  attenuated  rapidly.  Once  the  shaker  foot  to  sand  coupling  was 
changed  results  did  not  repeat  with  the  same  precision  as  when  it  was  left  alone. 
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CHAPTER  I 


INTRODUCTION 

Responsible  for  approximately  26,000  injuries  or  deaths  per  year,  an  estimated 
100  million  land  mines  lay  in  countries  throughout  the  world  [1].  Modem  land  mines  are 
constmcted  primarily  of  plastic,  making  detection  with  conventional  metal  detectors 
unreliable.  This  has  motivated  many  researchers  to  examine  new  techniques  of  mine 
detection.  One  such  technique  interrogates  for  the  presence  of  mines  using  high 
frequency  (30-2000  Hz)  seismic  waves.  Radar  is  used  as  a  non-contact  measure  of  the 
mine’s  seismic  signature.  This  eliminates  the  signal  to  noise  and  signal  to  reverberation 
problems  that  have  been  encountered  in  attempts  to  employ  conventional  pulse-echo 
techniques  to  the  seismic  detection  of  mines.  The  technique  requires  relatively  high 
frequency  seismic  excitations  in  order  to  interrogate  the  shallow  burial  depths  and  the 
small  dimensions  typical  for  landmines  (<12  inches)  [2]. 

Mine  detection  operations  commonly  occur  in  unconsolidated  soils.  Typically 
there  are  strong  nonlinearities  found  in  these  conditions.  The  purpose  of  this  research 
was  to  experimentally  characterize  these  features  using  moist,  compacted  sand  as  a  soil 
model.  Propagation  responses  were  measured  as  a  function  of  frequency,  amplitude,  and 
propagation  distance  over  the  ranges  of  interest  for  mine  detection.  Analysis  of  the 
generated  data  yielded  conclusions  about  the  following  subjects:  (a)  frequency  threshold 
dependency  on  the  amplitude  of  excitation  and  propagation  distance,  (b)  saturation 
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threshold  dependency  on  the  frequency  of  excitation  and  propagation  distance,  (c)  source 
and  near  source  nonlinearities  versus  propagation  path  nonlinearities,  (d)  differences 
between  various  acoustic  transducer  arrangements  (transducer  type,  frequency  response, 
motion  behavior,  sources  of  nonlinearities,  etc.)  and  (e)  possible  exploitation  of 
nonlinearities  to  enhance  characteristics  of  the  incident  signal  (level,  bandwidth, 
directivity,  duration,  etc.). 
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CHAPTER  II 


BACKGROUND 

General 

Finding  new  and  better  ways  of  detecting  land  mines  has  become  increasingly 
important  in  recent  years  due  primarily  to  the  large  numbers  actively  employed 
throughout  the  world  and  the  shift  towards  mines  manufactured  with  plastic  casings. 

Mine  warfare  has  also  gained  public  attention  through  recent  attempts  at  an  international 
ban  on  landmines,  changes  in  U.S.  landmine  usage  policy,  and  U.S.  deployments  to 
densely  mine-laden  countries  such  as  Bosnia. 

Current  technology  relies  heavily  upon  electromagnetic  detection  techniques.  The 
most  prevalent  is  common  metal  detection,  and  for  good  reason.  This  technique  is  very 
reliable  and  easy  to  employ  due  to  the  significant  differences  in  the  electromagnetic 
properties  of  the  metal  mine  and  the  ground.  Even  with  newer  mines  encased  in  plastic, 
there  is  still  a  detectable  albeit  fainter,  electromagnetic  signature  that  may  be  found  with 
this  equipment.  Ground  Penetrating  Radar  (GPR)  also  exploits  the  mine-ground 
electromagnetic  properties  successfully.  One  of  the  problems  with  these  techniques 
however  is  that  they  will  also  detect  every  soda  can,  coin,  bolt,  and  any  other  piece  of 
elctromagnetically  significant  trash  in  the  ground.  This  yields  a  large  percentage  of  false 
alarms  when  attempting  to  clear  an  area  of  mines. 
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As  researchers  attempt  to  improve  mine  detection  capabilities  by  finding  a  way  to 
augment  electromagnetic  exploitation  techniques,  many  have  turned  to  the  development 
of  acoustic  detection  techniques.  The  biggest  advantage  of  acoustic  detection  is  the  large 
contrast  between  the  acoustic  properties  of  the  ground  and  those  of  the  mine,  regardless 
of  metal  or  plastic  casing.  Buried  clutter  such  as  rocks,  sticks,  or  man-made  objects  also 
exhibit  significantly  different  acoustical  properties  than  landmines.  An  example  of  how 
different  the  ground  and  mine  acoustic  properties  are  is  the  shear  wave  velocity,  which  is 
approximately  20  times  greater  through  the  mine’s  explosives  than  it  is  through  the 
surrounding  soil  [2]. 

A  mine  is  also  a  complex  structure  that  includes  a  flexible,  smooth  plastic, 
wooden,  or  metal  casing;  a  quantity  of  explosive  materials;  a  firing  mechanism  (trigger, 
firing  pin,  and  volatile  initiating  charge);  and  air  pockets  (usually  surrounding  the  firing 
mechanism).  Each  of  these  components  may  vibrate  under  the  action  of  forces  inherent 
in  the  mine  without  an  externally  applied  force.  The  frequency  at  which  this  vibration 
occurs  is  called  the  natural  frequency.  If  a  frequency  of  excitation  coincides  with  any  of 
the  natural  frequencies  of  the  mine,  resonance  occurs.  At  resonance,  the  amplitudes  of 
motion  may  become  very  large  [3].  Because  the  mine  structure  is  so  complex  the 
probability  of  achieving  resonance  is  higher  than  it  would  be  for  a  simple,  homogeneous 
material  or  structure.  The  exaggerated  displacement  amplitudes  associated  with 
achieving  resonance  can  be  used  to  detect  the  location  of  the  mine  when  buried. 

Much  of  the  current  research  involves  pulse-echo  techniques  such  as  the  use  of 
echo  location  (direct  excitation  and  reception  of  seismic  waves)  [4],  or  detection  by 
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spectral  analysis  of  surface  waves  (SASW)  [5],  Some  of  the  problems  that  these 
techniques  have  encountered  are  practical  implementation  issues,  low  signal  to  noise 
ratios,  inability  to  differentiate  between  mines  and  debris  of  similar  acoustic  reflectivity 
[6],  and  significant  residue  hiding  the  object  reflection  when  incident-reflected  signal 
subtraction  is  used  [7],  Also,  in  the  case  of  SASW,  incident  surface  waves  arrive  at  the 
receivers  almost  simultaneously  with  the  reflected  waves  from  a  shallowly  buried  object 
such  as  a  mine. 

A  new  mine  detection  technique  that  has  shown  potential  for  minimizing  the 
previously  mentioned  pulse-echo  problems  utilizes  an  elastic  wave  and  an 
electromagnetic  sensor  [2],  An  elastic  wave  transducer  (shaker)  has  been  used  to 
generate  waves  in  the  ground.  An  electromagnetic  radar  is  used  to  measure  the  surface 
displacement  as  these  elastic  waves  travel  through  the  ground.  Because  the  mine  has 
very  different  mechanical  properties,  as  previously  mentioned,  the  unique  resonances  of 
the  mine  and  the  reflection  and  scattering  of  the  waves  cause  the  ground  above  it  to 
vibrate  differently  from  the  surrounding  ground.  The  radar  identifies  where  the  ground 
particle  displacements  are  different,  thus  identifying  where  the  mine  is  buried. 

This  technique  has  resolved  several  practical  implementation  issues  and  has 
shown  results  with  excellent  signal  to  noise  ratios.  Experiments  have  also  been 
performed  with  various  types  of  buried  clutter  such  as  rocks  and  sticks.  In  all  cases  the 
clutter  has  not  been  detected  which  indicates  that  this  method  would  greatly  reduce  the 
number  of  false  alarms  in  a  de-mining  operation.  The  technique  has  detected  mines 
buried  12  inches  deep  and  has  also  been  able  to  detect  mines  underneath  a  groundcover 
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(pine  straw).  In  addition  to  this  procedure  showing  promise  for  use  on  its  own,  another 
advantage  would  be  its  use  in  conjunction  with  other  detection  systems.  Because  the 
acousto-electromagnetic  detection  technique  senses  different  phenomena,  it  may  be 
employed  with  conventional  detection  techniques  such  as  metal  detection,  ground- 
penetrating  radar,  or  infrared  detection  in  order  to  synergistically  improve  the  chances  of 
successful  mine  detection. 

The  tests  of  this  system  are  being  done  in  a  large  “sandbox”.  Sand  was  chosen  as 
the  soil  medium  for  its  workability  when  burying  and  digging  up  mines,  and  its 
relevance.  Moist  sand  is  the  most  common  soil  encountered  by  Marines  when 
conducting  amphibious  landings  so  it  was  also  of  practical  interest.  Figure  2. 1  is  a 
picture  of  the  experimental  setup  showing  this  sandbox. 

Within  the  sandbox  there  is  a  120  cm  by  80  cm  area  referred  to  as  the  scan  region. 
The  scan  region  is  the  entire  area  used  to  take  surface  displacement  measurements  when 
looking  for  mines.  The  lead  edge  of  the  scan  region  is  approximately  30  cm  from  the 
shaker  foot.  The  radar  is  mounted  on  a  three-dimensional  positioner  that  moves  it  around 
the  scan  region.  The  x  axis  of  the  scan  region  is  defined  as  perpendicular  to  the  length  of 
the  rectangular  shaker  foot  and  the  y  axis  is  defined  as  parallel  to  the  length  of  the 
rectangular  shaker  foot.  Surface  displacement  measurements  are  taken  every  cm  in  the  x- 
direction  and  every  two  cm’s  in  the  y-direction  for  a  total  of  4,961  measurements.  This 
two-dimensional  scan  may  be  processed  in  time  in  order  to  create  a  time  progression,  or 
movie,  of  the  wave  propagating.  Different  colors  represent  the  magnitude  of  surface 
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Figure  2. 1  Photograph  of  the  Experimental  Setup 


Surface  Wave 
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displacements  thereby  showing  how  the  wave  propagates  and  identifying  where  the  mine 
is  located. 

The  shaker  generates  the  elastic  waves  utilizing  a  30  Hz  -  2000  Hz  chirp  of  3.6 
second  duration  and  a  rectangular  shaker  foot  that  sits  parallel  to  the  lead  edge  of  the  scan 
region.  These  relatively  high  frequencies  are  necessary  in  order  to  interrogate  the  soil  for 
smaller  anti-personnel  mines,  which  can  be  as  small  as  two  inches  in  diameter.  The 
shaker  produces  compressional,  shear,  and  surface  waves  within  the  sandbox. 

Figure  2.2  illustrates  the  different  types  of  waves  produced.  The  picture  comes 
from  a  computer  simulation  of  the  waves  [8].  It  is  a  cross-sectional  view  of  an  elastic 
half  space  with  a  point  source.  The  lower  half  represents  the  sand  containing  the  waves 
propagating  in  three  dimensions  (x,z,t).  The  half  space  is  bounded  on  top  by  a  free 
surface  which  separates  the  sand  from  the  air.  The  boundary  condition  at  this  surface  is 
pressure  =  0.  This  causes  an  impedance  mismatch  at  the  boundary  creating  a  pressure 
release  surface  (pcsand »  pcair). 

The  various  waves  are  defined  by  their  particle  motion.  The  compressional  wave 
propagates  by  means  of  pure  translational  particle  motion  (particle  volume  changes  but  is 
irrotational).  This  wave  travels  the  fastest  (250  m/s  for  the  model  in  Figure  2.2)  and  has 
the  lowest  surface  normal  displacement.  The  shear  wave  propagates  by  means  of  pure 
rotational  particle  motion  (particles  remain  equivoluminal).  It  travels  slower  than  the 
compressional  wave  but  slightly  faster  than  the  surface  wave  (87  m/s  for  the  model  in 
Figure  2.2).  The  particle  motion  of  the  surface  wave  is  elliptical  in  that  it  contains  both 
translational  and  rotational  components  of  motion.  It  travels  at  a  speed  very  close  to,  but 


8 


slightly  less  than,  the  shear  wave  and  is  confined  to  approximately  one  wavelength  from 
the  free  surface.  It  propagates  cylindrically  while  the  compressional  and  shear  waves 
propagate  spherically  [9],  The  shear  and  surface  waves  cannot  be  separated  within  the 
sandbox  because  the  dimensions  are  such  that  the  waves  do  not  propagate  far  enough  for 
the  waves  traveling  at  almost  identical  speeds  to  separate. 

Another  way  to  view  the  surface  displacements  is  by  using  a  waterfall  plot,  which 
is  also  called  a  seismogram.  Waterfall  plots  are  generated  by  reading  the  surface 
displacements  with  a  one-dimensional  scan.  The  radar  reads  the  surface  displacement  at 
a  point  in  front  of  the  shaker  and  records  the  information  in  the  time-domain.  A  fast 
Fourier  transform  (FFT)  algorithm  then  takes  the  measurement  into  the  frequency- 
domain.  During  post-processing,  this  is  convolved  with  the  FFT  of  a  differentiated 
Gaussian  pulse  and  then  taken  back  to  the  time-domain  with  an  inverse  fast  Fourier 
transform  (IFFT)  algorithm.  The  result  is  a  pulse  that  represents  the  waveform  in  the 
time-domain.  The  radar  then  moves  to  the  next  position  along  the  axis  of  interest  and  the 
process  is  repeated.  Once  all  of  the  pulses  are  recorded,  they  are  plotted  one  above  the 
other  so  that  the  bottom  line  is  the  data  taken  at  the  closest  position  and  the  top  line  is  the 
data  taken  at  the  furthest  position.  The  y  axis  of  the  waterfall  plot  represents  the  distance 
from  the  shaker. 

Figure  2.3  shows  six  different  waterfall  plots  side  by  side  [2].  For  all  of  the 
waterfall  plots,  121  measurements  were  taken  along  the  x-axis  (perpendicular  to  the 
shaker  foot)  at  one  cm  intervals.  The  first  plot  was  done  with  no  mine  present  which  is 
also  referred  to  as  a  clean  scan.  The  other  plots  have  mines  present  in  the  vicinity  of  the 
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Figure  2.3  Waterfall  Plots  (Clean  Scan  and  with  Mines) 
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Figure  2.4  Computer  Simulation  Detecting  Buried  Mine 
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shaded  area.  The  pulse  amplitude  is  larger  in  these  regions  due  to  the  greater  surface 
displacements  caused  by  the  mine  resonances.  These  figures  illustrate  the  value  of  one¬ 
dimensional  scans  and  also  show  what  the  waveform  looks  like  as  it  propagates, 
attenuates,  and  disperses  through  the  sand  when  there  is  no  mine  present  and  how  the 
waveform  is  affected  in  the  presence  of  a  mine. 

Efforts  to  model  these  waves  and  their  interaction  with  mines  are  also  underway. 
In  order  to  do  this  a  two-dimensional  finite-difference  model  for  elastic  waves  has  been 
developed.  The  model  develops  a  first  order  stress-velocity  formulation  from  the 
equation  of  motion  and  the  stress-strain  relations.  The  system’s  equations  are  then 
discretized  using  centered  finite-differences.  The  outward  traveling  waves  are  absorbed 
by  a  perfectly  matched  layer  surrounding  the  discretized  solution  space  [8].  The  results 
of  this  modeling  have  proven  to  be  fairly  accurate. 

Figure  2.4  is  a  waterfall  plot,  for  a  mine  buried  at  x  =  90  cm,  generated  by  this 
computer  model  [8].  Although  the  two-dimensional  finite-difference  model  yields 
reasonably  accurate  results,  it  does  not  take  into  consideration  nonlinearities  present  in 
the  system.  This  figure  shows  the  mine  resonance  and  a  reflected  wave.  The  results  of 
the  computer  model  are  much  smoother  than  an  actual  one-dimensional  scan  however. 
The  same  levels  of  dispersion  are  not  present  and  the  reflected  pulse  is  much  more 
noticeable  in  the  computer  model.  The  model’s  accuracy  could  be  improved  if 
nonlinearities  were  taken  into  consideration. 
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Literature  Review 


Current  literature  on  wave  propagation  in  unconsolidated  soil  is  incomplete  for 
the  mine  detection  system  discussed.  There  are  two  main  reasons  for  this.  The  first  is 
inadequate  information  about  the  frequency  range  of  interest.  Seismology  research  is 
concerned  primarily  with  very  low  frequencies.  The  acousto-electromagnetic  form  of 
detection  however,  utilizes  a  bandwidth  of  30-2000  Hz  which  is  regarded  as  noise  by 
many  seismologists.  The  reason  for  this  bandwidth  is  twofold.  One,  frequencies  less 
than  100  Hz  are  extremely  susceptible  to  dispersion  and  experience  much  greater 
attenuation  between  5  and  10  meters  [4].  The  wavelengths  associated  with  the  surface 
wave  for  frequencies  lower  than  100  Hz  are  also  not  suitable  for  detecting  smaller  mines. 
Frequencies  higher  than  this  do  not  propagate  far  enough  to  be  useful.  Two,  the  purpose 
of  the  launched  ground  wave  is  to  excite  the  resonances  of  the  mine.  The  greater  the 
bandwidth  of  the  incident  signal,  the  more  modes  of  vibration  in  the  mine  that  are 
excited. 

The  second  reason  that  the  current  literature  is  incomplete  is  the  depth  at  which 
the  mine  is  buried.  Seismological  research  has  not  been  concerned  with  shallow  depths 
(<  12  inches).  Underground  obstacle  or  cavity  detection  has  focussed  much  deeper. 
Some  current  methods  of  detection  encounter  problems  when  applied  this  close  to  the 
surface  as  was  the  case  with  the  SASW  method. 

Although  current  seismological  research  is  insufficient  to  support  this  method  of 
mine  detection,  it  does  yield  valuable  insight  into  the  particle  behavior  observed  in  the 
sandbox.  The  equations  of  motion  that  apply  to  homogeneous,  isotropic,  elastic  media 
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can  be  written  in  terms  of  stresses  as  shown  in  Equation  2. 1  which  applies  to  the  x 
direction  only: 

p(5^u/5t^)  =  (dox/dx)  +  (dxxy/dy)  +  (dxxz/dz)  Eqn.  2. 1  [9] 

In  order  to  express  the  right-hand  side  of  this  equation  in  terms  of  displacements, 
relationships  for  the  stresses  in  terms  of  Earners  constant,  the  shear  modulus,  cubical 
dilatation,  and  shear  strain  are  used.  Combining  these  with  relationships  for  strain  and 
rotation  in  terms  of  displacement  yield  Equation  2.2: 

pid^u/dt^)  =  (>i+G)(5E/ax)  +  GV^u  Eqn.  2.2  [9] 

There  are  two  solutions  for  this  equation  of  motion.  The  first  applies  to  the 
propagation  of  a  dilatational  wave  (compressional  wave  or  P-wave).  The  second  applies 
to  a  distortional  wave  (shear  wave  or  S-wave).  Because  the  experiments  in  the  sandbox 
are  in  an  elastic  half-space,  a  third  type  of  wave  is  also  present.  This  wave  is  the 
Rayleigh  or  surface  wave.  A  surface  wave  is  confined  to  the  surface  of  the  elastic  half¬ 
space  and  contains  both  x  and  z  direction  displacements.  If  potential  functions  and  Y 

are  chosen  to  correspond  to  dilatation  and  rotation  of  the  wave  respectively,  then 
Equation  2.3  is  the  equation  of  motion  for  the  surface  wave  derived  from  Equation  2.2. 
p(d/dK)id^^/dtY  p(a/az)(a^T/at^)=(>.+2G)(a/ax)(v^«))+G(a/az)(v^'P)  Eqn.  2.3  [9] 
The  compressional  wave  has  the  highest  velocity.  The  shear  waves  and  surface 
waves  are  slower  and  their  velocities  are  very  similar.  The  surface  wave  travels 
approximately  94  %  as  fast  as  the  shear  wave.  The  particle  motion  of  the  compressional 
wave  is  in  the  x  direction  and  the  shear  wave  particle  motion  is  in  the  z  direction.  If  the  x 
and  z  components  of  the  surface  wave  are  added  the  particle  motion  is  in  a 
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counterclockwise  direction  as  the  wave  propagates  to  the  right.  The  percent  of  total 
energy  carried  by  each  of  these  three  waves  was  measured  for  a  circular  footing  and 
found  to  be  distributed  as  such:  surface  wave  67%,  shear  wave  26%,  and  compressional 
wave  7%  [9],  As  the  waves  propagate,  the  compressional  and  shear  wave  amplitudes 
decrease  as  f  *  due  to  spherical  spreading.  The  surface  wave  amplitude  decreases  as 
due  to  cylindrical  spreading. 

Biot  developed  stress-strain  relationships  for  wave  propagation  in  porous 
saturated  solids.  He  discovered  that  there  is  only  one  type  of  shear  wave  that  propagates 
through  the  elastic  structure  because  there  is  no  structural  coupling  between  the  elastic 
structure  and  the  fluid.  This  is  because  there  is  no  shearing  stiffness  in  the  fluid.  There 
are  two  compressional  waves  however.  One  propagates  through  the  elastic  medium  and 
the  other  through  the  fluid.  They  are  coupled  by  the  stiffnesses  and  motions  of  the  elastic 
medium  and  fluid.  The  compressional  wave  propagating  in  the  fluid  is  the  fastest.  It  is 
faster  than  if  it  were  traveling  in  fluid  alone  due  to  a  pushing  effect  by  the  elastic 
medium.  The  next  fastest  is  the  compressional  wave  in  the  elastic  medium.  This  is 
slower  than  if  the  medium  were  dry  because  of  the  drag  caused  by  the  water  in  the  pores 

[9]. 

The  horizontal  water  table  affects  wave  propagation  also.  If  the  water  table  is 
close  enough  to  the  surface,  then  reflected  and  refracted  waves  from  this  boundary  can 
become  a  factor.  It  can  also  influence  the  wave  velocities  depending  on  whether  the 
measurements  are  taken  below  or  above  the  water  table.  The  amount  of  air  in  the  pores 
of  the  half-space  makes  a  large  impact  on  the  wave  velocity.  Richart  gives  an  example  in 
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which  a  0.1%  increase  in  air  bubbles  slows  wave  speed  from  4800  ft/sec  to  1204  ft/sec. 
This  is  important  to  note  because  the  sand  being  used  for  the  mine  detection  experiments 
is  not  saturated  completely  by  water.  The  water  table  is  a  couple  feet  below  the  surface 
while  the  sand  at  the  surface  and  through  the  depths  pertaining  to  mine  detection  is  damp. 
There  is  plenty  of  air  in  the  pores  of  the  sand  because  after  it  is  watered  the  sand  has 
room  to  drain  to  the  lower  water  table. 

Some  relevant  material  was  found  done  by  Sabatier  [10,  1 1]  in  the  area  of 
acoustic-seismic  coupling.  Within  this  research  are  conclusions  about  the  frequency  and 
depth  dependence  of  attenuation.  In  general,  Sabatier  found  that  higher  frequencies  in 
the  ground  attenuate  faster  than  lower  frequencies.  For  example,  200  Hz  attenuates  at  a 
rate  of  0.06  -  0.1  dB/cm  but  1200  Hz  attenuates  at  a  rate  of  0.16  -  0.22  dB/cm  [10]. 
Sabatier  also  did  tests  with  a  speaker  source  in  the  air  and  a  buried  microphone  to 
determine  the  effects  of  frequency  and  depth  on  the  wave  attenuation.  He  found  that  at  5 
cm  below  the  surface  waves  at  1  Hz  were  attenuating  at  approximately  2  dB/cm  while 
waves  of  1000  Hz  at  the  same  point  were  attenuating  at  approximately  14  dB/cm.  The 
same  test  was  done  10  cm  below  the  surface.  This  time  waves  of  1  Hz  were  attenuating 
at  approximately  3  dB/cm  and  waves  of  1000  Hz  were  attenuating  at  30  dB/cm  [1 1]. 
These  results  showed  that  wave  attenuation  is  dependent  on  both  frequency  and  depth 
and  that  the  effects  of  frequency  and  depth  are  related.  Waves  attenuate  at  a  more  rapid 
rate  the  deeper  they  travel  and  this  rate  of  attenuation  is  greater  for  higher  frequencies 
than  for  lower  frequencies. 
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Sabatier  has  also  employed  acoustic  techniques  to  the  detection  of  buried  objects 
[12].  The  system  used  in  this  study  included  a  sound  source  that  was  above  the  ground. 
Acousto-seismo  coupling  was  relied  upon  to  transmit  a  wave  through  the  ground.  The 
receiver  was  a  microphone  that  recorded  the  reflected  signal.  This  signal  contained  the 
reflection  from  the  surface  and  a  reflection  of  smaller  amplitude  from  a  buried  object. 

The  tests  proved  effective  for  objects  buried  less  than  five  cm  deep.  Some  results  that 
could  be  a  problem  when  applying  this  to  mine  detection  were  also  found.  The  type  of 
porous  media  made  a  difference  as  to  how  pronounced  the  reflected  signal  from  a  buried 
object  was,  a  significant  signal  is  reflected  from  a  hole  (false  alarms),  and  smearing  of  the 
surface  reflection  and  buried  object  reflection  occurred. 

Objective 

The  objective  of  this  thesis  was  to  characterize  the  nonlinearites  of  the 
propagation  path  so  that  the  results  may  then  be  used  to  refine  the  experimental  procedure 
for  acousto-electromagnetic  mine  detection  and  be  applied  to  the  finite-difference  time- 
domain  computer  model.  An  example  of  an  observed  nonlinearity  that  prompted  this 
research  was  a  graph  of  displacement  versus  input  voltage  in  which  the  curve  rose 
linearly,  as  expected,  until  it  began  to  saturate.  When  the  voltage  continued  to  increase, 
the  curve  eventually  began  to  rise  again  until  it  reached  a  second  point  of  saturation. 

When  this  same  test  was  done  at  a  different  point  in  the  box,  the  voltage  where  the 
saturation  effects  were  seen  was  different.  This  indicated  that  the  phenomenon  was  a 
result  of  nonlinearities  in  the  propagation  path  and  not  the  source. 
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Another  indicator  of  nonlinearities  in  the  system  included  several  harmonics 
being  produced  in  the  frequency  response  of  the  surface  displacements.  How  much  the 
propagation  path  contributed  to  these  nonlinearities  versus  how  much  the  source  or 
source  to  sand  coupling  contributed  was  unknown.  The  goal  was  to  characterize  the 
nonlinearities  of  the  propagation  path  but  this  could  not  be  done  without  considering  the 
entire  system.  Nonlinearities  may  occur  in  the  power  amplifier,  impedance  mismatches 
between  the  power  amplifier  and  the  shaker  resulting  in  oscillating  transfer  functions,  the 
motion  of  the  shaker  foot,  the  shaker  foot  to  sand  coupling,  along  the  propagation  path,  or 
at  the  receiver.  These  sources  of  nonlinearities  must  be  isolated  to  determine  which  are 
routinely  encountered,  the  relative  impact  of  these  nonlinear  effects,  how  these  effects 
may  be  reduced  or  eliminated,  and  which  of  these  effects  should  be  considered  in  the 
computer  model. 

The  data  collected  fell  into  two  major  categories:  1)  frequency  response  as  a 
function  of  drive  amplitude  and  distance  between  the  source  and  receiver  and  2)  surface 
displacement  as  a  function  of  drive  frequency  and  distance  between  the  source  and 
receiver.  The  frequency  response  data  would  show  how  the  drive  amplitude  of  the 
incident  wave  and  the  propagation  path  affect  the  surface  displacement  as  a  function  of 
frequency  for  a  fixed  distance.  The  surface  displacement  data  would  show  how  the 
frequency  and  the  propagation  path  affect  the  surface  displacement  as  a  function  of  the 
drive  amplitude  for  a  fixed  distance.  In  both  cases,  the  fundamental  driving  frequency 
and  its  harmonics  were  evaluated. 


17 


CHAPTER  m 


INSTRUMENTATION  AND  EQUIPMENT 

Figure  3.1  shows  the  general  configuration  of  the  experimental  arrangement  for 
Experiment  Two.  The  only  difference  in  Experiment  One  was  that  accelerometers  were 
not  used.  Table  3.1  lists  the  major  components  of  the  experimental  setup,  the 
manufacturers  of  the  equipment,  and  the  equipment  models.  This  configuration  is  the 
same  one  used  when  scans  are  performed  for  the  acousto-electromagnetic  detection  of 
mines. 

The  sandbox  is  approximately  4.5  m  long,  4.5  m  wide,  and  1.5  m  deep.  It  is 
wedge  shaped  as  seen  in  Figure  2.1.  The  end  where  the  shaker  sits  is  approximately  five 
feet  deep.  This  depth  extends  across  half  of  the  sandbox  and  then  slopes  up  to  the  top  of 
the  far  end.  It  contains  about  50  tons  of  packed  sand  with  a  water  table  that  averages  two 
feet  below  the  surface.  The  dimensions  are  such  that  measurements  in  the  scan  region  do 
not  record  reflections  from  the  sides. 

The  elastic  wave  transducers  are  20  pound  and  100  pound  shakers.  They  make 
contact  with  the  sand  through  the  use  of  a  shaker  “foot”  which  can  be  different  sizes  and 
shapes.  The  radar  is  a  homodyne  type  mounted  on  an  XYZ  positioner.  The  data 
acquisition  and  positioning  is  automated  using  Lab  VIEW  code  to  control  the  generated 
signal,  the  position  of  the  radar,  the  various  circuitry,  and  the  collection  and  processing  of 
data.  The  major  component  details  are  listed  below. 
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Figure  3.1  Experimental  Setup 


NAME 

MANUFACTURER 

MODEL 

REMARKS 

Data  Acquisition  Card 

National  Instruments 

PCI-MIO- 

16E-1 

1.25  MS/s  (single 
channel) 

Multi-channel  filter 

Krohn-Hite 

3944 

Filters  out  <  30  Hz 

Stanford  Research  System 

SR560 

Stanford  Research  System 

SR560 

Power  supply 

Top  ward 

ESeuHH 

Accelerometers 

Power  supply 

Topward 

Radar 

Amplifier 

Crown 

CE2000 

Modified 

Shaker 

Vibration  Test  Systems 

VG  100-6 

1001b 

Power  meter 

Hewlett-Packard 

437B 

Radar 

Home-made 

Accelerometers 

Kistler 

Sensitivity;  3.4- 
3.6mV/g 

Accelerometers 

PCB  Piezotronics 

352C67 

Sensitivity: 

109.5mV/g 

Accelerometers 

PCB  Piezotronecs 

352B22 

Sensitivity;  9.3- 
10.6mV/g 

Table  3.1  Experimental  Component  Details 
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Software 


The  collection  of  data  and  operation  of  the  instrumentation  was  coded  in 
Lab  VIEW.  The  Lab  VIEW  codes  are  contained  in  Chapter  XIII  -  Appendix  E.  The  data 
processing  and  analysis  was  done  with  MATLAB.  The  various  MATLAB  codes  are 
found  in  Chapter  XII  -  Appendix  D. 

Data  Acquisition  Card 

The  data  was  collected  with  the  use  of  a  National  Instruments  Data  Acquisition 
Card  (DAC).  The  PCI-MIO-16E-1  model  card  used  had  the  following  analog  input 
characteristics:  16  single-ended  or  8  differential  channels,  12  bit  resolution,  and  a 
maximum  single  channel  sampling  rate  of  1.25  MS/s.  It  had  a  1.6  MHz  bandwidth  for 
small  (-3dB)  signals  and  a  1  MHz  bandwidth  for  large  (1%  total  harmonic  distortion) 
signals.  The  DAC  had  the  following  analog  output  characteristics:  2  voltage  channels, 
12  bit  resolution,  and  a  maximum  single  channel  update  rate  of  IMS/s.  Digital 
input/output  had  8  input/output  channels. 


Radar 

The  homodyne  type  radar  measured  surface  displacements  by  executing  a  phase 
comparison.  It  operated  with  a  power  of  1  W  and  had  a  sensitivity  of  1  nm.  The  spot 
size  was  approximately  2  cm  in  diameter.  This  spot  size  would  limit  the  accuracy  of 
measurements  above  3000  Hz  due  to  the  small  wavelengths.  A  power  meter  was 
mounted  above  the  radar  to  monitor  whether  or  not  the  radar  was  operating  within  its  5 
dBm  to  15  dBm  optimal  range  (-5  dBm  to  5  dBm  as  read  on  the  power  meter). 
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Positioner 


The  positioner  is  mounted  approximately  1.5  m  above  the  surface  of  the  sand.  It 
is  capable  of  moving  in  the  x,  y,  or  z  directions  at  various  rates  and  ranges.  With  the 
experimental  arrangement  currently  employed,  the  positioner  is  limited  to  190  cm  in  the  x 
direction.  It  can  easily  range  120  cm  in  the  y  direction  and  30  cm  in  the  z  direction  which 
in  no  way  limits  current  data  collection  processes.  All  data  for  Experiments  One  and 
Two  were  taken  along  the  x  axis  with  the  waveguide  positioned  1.3  to  1.8  cm  above  the 
sand  surface. 


Shaker 

The  signal  that  drives  the  shaker  is  produced  in  Lab  VIEW  and  sent  through  the 
break-out  box  to  a  Crown  CE2000  amplifier.  The  signal  is  then  sent  to  a  shaker  that 
makes  contact  with  the  sand  through  the  use  of  a  shaker  foot.  All  measurements  for 
Experiment  One  were  made  with  a  20  pound  shaker  using  a  rectangular  foot  measuring 
21.6  cm  in  length  by  1.3  cm  in  width  (surface  area  =  28.1  cm^).  During  Experiment  Two, 
three  shaker-foot  combinations  were  utilized.  The  first  was  the  same  shaker  and  foot 
used  in  Experiment  One.  The  second  was  the  20  pound  shaker  and  a  round  foot  with  a 
diameter  of  10.2  cm  (surface  area  =  81.7  cm^).  The  third  was  a  VG  100-6  shaker, 
capable  of  producing  a  force  of  100  pounds,  and  a  rectangular  foot  measuring  30.2  cm  in 
length  and  3 .2  cm  in  width  (surface  area  =  96.6  cm^).  A  blower  cools  the  shaker 
throughout  its  operation.  Figure  3.2  is  a  picture  of  the  different  shakers  and  feet.  Figure 
3.3  is  a  picture  of  the  accelerometer  placement  for  Experiment  Two. 
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Figure  3 .2  lOOIb  Shaker  with  Rectangular  Foot  Mounted  (left  -  view  from  side), 
201b  Shaker  with  Circular  Foot  Mounted  (right  top  -  view  from  bottom).  Small 
RectangularFootfor201b  Shaker  (right  bottom  -  view  from  side) 


Figure  3 .3  Accelerometer  Placement  on  the  1 001b  Shaker  with  Rectangular  Foot 
(Bottom  View)  -  3  PCB  352B22  accelerometers  to  record  vertical  acceleration  and 
1  PCB  3  52C67  accelerometer  to  record  horizontal  acceleration 


CHAPTER  IV 


EXPERIMENT  ONE 

The  objective  of  this  experiment  was  to  measure  the  surface  displacements  of  the 
sandbox  as  a  function  of  drive  amplitude,  drive  frequency,  and  propagation  distance.  The 
resulting  data  would  be  used  to  determine  both  the  frequency  responses  and  amplitude 
responses  for  the  fundamental  frequency  and  harmonics.  There  were  also  two  secondary 
objectives.  First,  data  would  be  collected  in  order  to  separate  the  compressional  and 
surface  waves.  Second,  data  would  be  taken  with  a  different  shaker  foot  to  try  to  alter  the 
relative  content  of  pressure  wave  and  surface  wave.  As  a  result  of  collecting  and 
processing  this  data,  many  ways  to  improve  the  data  collection  were  found.  This  lead  to 
Experiment  Two  which  is  discussed  in  Chapter  VI. 

Data  was  taken  the  following  way  in  order  to  determine  the  frequency  and 
amplitude  response  of  the  sand.  Continuous  wave  (CW)  signals  from  33  Hz  to  2002  Hz, 
at  1 1  Hz  increments,  were  used  at  a  given  amplitude  and  position.  It  was  already  known, 
from  the  mine  detection  experiments,  that  the  waves  propagated  through  the  box  in  less 
than  0.07  second.  Therefore  a  frequency  increment  of  14  Hz  (1/0.07)  provided  ample 
resolution  to  document  the  wave  propagation.  It  was  also  known  that  60  Hz  and  its 
harmonics  were  very  large  in  the  noise  floor.  By  choosing  1 1  Hz  as  the  frequency 
increment  however,  of  the  180  frequencies,  only  three  coincided  with  a  multiple  of  60  Hz 
and  the  first  would  not  occur  until  660  Hz.  This  reduced  the  impact  of  the  60  Hz  noise. 
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Once  the  frequencies  had  been  swept  through,  the  amplitude  increased  linearly 
and  the  frequencies  swept  through  again.  A  total  of  24  amplitudes  were  used  ranging 
from  0.06  volts  to  0.96  volts  (input  to  power  amplifier)  at  a  given  setting  on  the  power 
amplifier.  All  voltages  are  peak  unless  otherwise  specified.  Once  the  frequencies  had 
swept  through  the  given  range  for  each  of  the  24  amplitudes,  the  radar  was  moved  to  a 
different  location  and  the  procedure  was  repeated.  This  process  was  done  at  the 
following  locations  in  the  sandbox:  x  =  40  cm,  x  =  80  cm,  and  x  =  120  cm,  all  along  the 
X  axis  (origin  located  26  cm  from  shaker  foot).  The  idea  behind  collecting  the  data  like 
this  was  that  both  the  frequencies  and  amplitudes  used  were  dense  enough  to  generate 
both  the  frequency  response  and  the  amplitude  response  from  the  same  data.  The  details 
of  this  experiment,  including  the  experimental  design,  data  collection  procedure,  and 
results,  may  be  found  in  Chapter  IX  -  Appendix  A. 

Several  lessons  were  learned  when  this  data  was  processed.  First,  the  voltages 
ranging  from  0.06  to  0.96  were  only  meaningful  for  producing  total  system  transfer 
functions  for  a  given  power  amplifier  setting.  They  could  also  be  normalized  by  the  data 
for  0.06  V  in  order  to  study  relative  effects.  It  would  have  been  much  more  useful 
however  to  know  what  force  the  shaker  foot  applied  to  the  sand.  That  way  frequency 
response  as  a  function  of  the  shaker  force  would  be  known.  This  involved  measuring  the 
current  input  to  the  shaker  (could  be  used  to  calculate  power  amplifier  transfer  function) 
and  using  the  shaker  specifications  to  realize  what  force  was  applied  given  the  input 
voltage  from  Lab  VIEW.  This  was  done  in  Chapter  VI  -  Experiment  Two. 
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The  second  lesson  learned  was  that  the  method  for  processing  the  data  produced  a 
significant  leakage  of  the  fundamental  and  harmonics  into  the  surrounding  frequencies. 
This  was  evidenced  by  the  noise  level  rising  with  the  increasing  amplitude  of  the  incident 
signal.  A  three  and  a  half  second  incident  signal  was  used  so  that  any  ring-up  or  ring- 
down  transients  could  be  cut  out  of  the  data.  Only  two  seconds  from  the  middle  of  the 
signal  was  saved  but  this  was  not  being  cut  out  at  an  integer  number  of  periods  for  each 
frequency.  When  this  piece  of  the  signal  was  taken  into  the  frequency-domain,  the 
additional  amount  past  the  integer  number  of  cycles  contained  frequency  components 
different  than  the  continuous  wave.  These  components  showed  up  in  the  noise  floor. 
Figure  4. 1  shows  this  rising  noise  level  for  two  different  amplitudes.  This  problem  was 
corrected  for  Experiment  Two  as  described  in  Chapter  VI. 

The  third  lesson  learned  was  that  the  amplitudes  increased  to  a  point  where 
dynamic  fluidization  occurred.  This  happened  around  amplitudes  in  the  vicinity  of  0.82 
V.  Once  the  shaker  foot  would  bury  itself  in  the  sand,  it  would  be  removed,  the  sand 
would  be  repacked,  and  the  shaker  would  be  placed  back  on  the  ground.  Although  this 
was  not  an  anticipated  problem,  once  the  shaker  was  moved  the  results  were  not 
repeatable.  The  shaker  foot  to  sand  contact  was  different  every  time  the  shaker  was 
placed  in  the  sandbox  and  therefore  measurements  needed  to  be  taken  without  moving 
the  shaker  until  the  data  collection  was  complete.  This  precluded  the  use  of  the  high 
amplitudes  that  caused  dynamic  fluidization  or  a  different  shaker  foot  was  needed. 

The  problem  of  burying  the  shaker  foot  was  best  seen  in  the  amplitude  response 
measurements  such  as  the  one  shown  in  Figure  4.2.  The  curve  begins  relatively 
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smoothly.  Everywhere  there  is  an  x  on  the  curve  is  where  the  shaker  had  been  picked  up 
in  order  to  repack  the  sand  under  it.  These  places  yielded  very  large  jumps  in  the  curve 
because  the  new  shaker  foot  to  sand  contact  and  the  sand  particle  matrix  under  the  foot 
was  significantly  different.  The  method  of  collecting  the  data  for  Experiment  Two 
corrected  this  problem  as  described  in  Chapter  VI. 

The  fourth  lesson  learned  involved  the  range  of  amplitudes  being  used.  The  lower 
end  of  the  amplitudes  was  not  low  enough  to  ensure  that  the  data  collection  was 
beginning  in  a  linear  region  prior  to  becoming  nonlinear.  There  were  also  not  enough 
data  points  to  generate  a  smooth  amplitude  response.  On  the  other  hand,  there  were 
many  more  amplitudes  than  required  to  generate  the  fi-equency  response  for  increasing 
amplitudes.  This  also  changed  the  way  data  was  collected  for  Experiment  Two  as 
described  in  Chapter  VI. 

The  next  phase  of  Experiment  One  was  to  collect  data  in  order  to  separate  the 
compressional  wave  from  the  surface  wave.  Data  was  taken  in  the  same  format  as 
described  above  but  it  was  taken  on  the  x  axis  for  x  =  190  cm.  The  intent  was  to  let  the 
waves  propagate  far  enough  that  they  would  separate  into  distinguishable  pressure  and 
surface  waves  in  the  time  domain.  From  here  the  waves  would  be  zeroed  out  one  at  a 
time  while  the  other  was  taken  back  into  the  frequency  domain.  The  same  frequency 
responses  and  amplitude  responses  would  be  produced  for  each  individual  component  to 
measure  the  relative  contribution  of  each. 

Two  problems  arose  during  this  procedure.  The  first  was  a  linear  assumption 
used  when  taking  the  data  into  the  time  domain.  The  frequency  response  for  33  Hz  to 
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2002  Hz  was  convolved  with  the  FFT  of  a  differentiated  Gaussian  pulse  as  described  for 
the  mine  detection  data  processing  in  Chapter  n  -  Background.  This  was  then  taken  into 
the  time-domain  with  an  IFFT  algorithm.  This  processing  method  did  not  take  into 
account  that  any  given  frequency  may  have  been  producing  harmonics  that  would  show 
up  in  another  frequency  bin.  Because  of  the  significant  nonlinearities  in  the  system  this 
assumption  was  not  going  to  be  reasonable. 

The  second  problem  came  in  trying  to  determine  at  what  point  the  surface  wave 
ended  and  the  pressure  wave  began.  Although  the  waves  had  propagated  through  a 
considerable  distance,  significant  dispersion  seemed  to  occur.  This  caused  the  lagging 
compressional  wave  fronts  to  remain  close  to  the  leading  surface  wave  fronts.  In  addition 
to  this  the  compressional  wave  had  attenuated  to  a  point  where  it  was  difficult  to  see. 

This  was  a  result  of  having  taken  data  at  x  =  190  cm  instead  of  doing  a  one  dimensional 
scan  out  to  190  cm  and  then  using  a  waterfall  plot  to  trace  the  progress  of  the 
compressional  wave.  The  results  of  this  phase  of  Experiment  One  led  to  the  conclusion 
that  the  compressional  and  surface  waves  could  not  be  separated  to  determine  their 
relative  contribution  to  the  frequency  response  of  the  sand. 

Experiment  One  ended  at  this  point.  The  problems  encountered  prompted  a 
significant  redesign  of  the  data  collection  and  processing  procedure.  The  new 
experimental  procedure  needed  to  make  use  of  the  lessons  learned  in  Experiment  One  to 
fix  the  leakage  during  data  acquisition,  increase  the  dynamic  range  of  amplitudes  used  to 
generate  the  amplitude  response,  and  prevent  the  shaker  from  causing  dynamic 
fluidization  or  having  to  be  moved  during  the  data  collection.  A  study  of  various  shakers 
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and  shaker  feet  was  performed  to  try  to  fix  this  last  problem.  The  results  of  this  study  are 
found  in  Chapter  V  -  Acoustic  Transducers.  The  lessons  learned  from  Experiment  One 
and  the  study  of  shakers  were  incorporated  into  Experiment  Two  in  order  to  correct  the 
above  mentioned  problems.  The  way  this  was  done  is  addressed  in  Chapter  VI  - 
Experiment  Two. 
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CHAPTER  V 


ACOUSTIC  TRANSDUCERS 

One  of  the  biggest  problems  with  Experiment  One  was  the  fact  that  the  shaker 
foot  was  experiencing  dynamic  fluidization  at  the  moderately  high  amplitudes.  It  would 
also  bury  in  a  very  short  amount  of  time  at  frequencies  less  than  100  Hz.  This  was  due 
mainly  to  the  small  surface  area  of  the  foot  as  opposed  to  the  duration  of  the  experiments. 
It  was  also  assumed  that  a  significant  amount  of  the  nonlinearites  in  the  system  were 
being  produced  either  by  the  shaker  foot  or  due  to  the  shaker  foot  to  sand  contact.  The 
objective  of  the  experiments  was  to  determine  characteristics  of  nonlinearities  in  the  sand 
however.  In  order  to  do  this,  nonlinearities  produced  at  the  shaker  foot  needed  to  be 
quantified  and  minimized. 

The  above  mentioned  problems  prompted  a  study  of  the  acoustic  transducers. 

The  purpose  of  the  study  was  to  characterize  the  motion  of  different  shaker  feet, 
determine  propagation  characteristics  for  each  of  their  radiated  waves,  and  select  the  one 
with  the  fewest  nonlinearities  to  conduct  Experiment  Two.  Three  shaker  -  shaker  foot 
combinations  were  examined  in  detail.  The  first  was  the  arrangement  used  in  Experiment 
One  consisting  of  the  20  pound  shaker  with  the  small  rectangular  foot  (surface  area  = 

28.1  cm^).  The  second  was  the  20  pound  shaker  with  a  circular  foot  (surface  area  =  81.7 
cm^).  The  third  was  a  100  pound  shaker  with  a  large  rectangular  foot  (surface  area  =  96.6 
cm^).  This  last  combination  was  the  one  chosen  to  conduct  Experiment  Two. 
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In  order  to  characterize  the  behavior  of  the  different  shaker  feet,  data  was  taken  in 
a  three-step  process.  The  first  step  was  to  look  for  resonance  of  the  foot.  The  frequency 
response  and  amplitude  response  of  the  foot  under  no  load  (shaker  turned  upside  down  so 
that  foot  is  in  air)  was  measured.  This  was  done  utilizing  four  accelerometers  three  of 
which  were  PCB  352B22’s  and  one  of  which  was  a  PCB  352C67.  During  the  tests 
conducted  under  no  load  the  power  amplifier  was  put  on  one  gain  setting  lower  than  that 
of  the  tests  done  loaded.  This  prevented  the  foot  from  bottoming  out  the  suspension  of 
the  shaker.  The  second  step  was  to  look  for  sand  loaded  resonances.  The  same 
measurements  were  taken  under  loaded  conditions  (shaker  foot  placed  on  sand  as  it  was 
during  the  experiments).  All  four  accelerometers  were  left  in  position  from  step  one  to 
step  two.  The  third  step  was  to  measure  the  shaker  -  shaker  foot  combination’s 
propagation  behavior  by  taking  a  sample  amplitude  response  using  the  radar  for  a  given 
frequency  at  two  different  positions  in  the  sandbox  (x  =  10  cm  and  x  =  40  cm). 

For  each  of  the  frequency  response  data  sets,  frequencies  between  33  and  2002  Hz 
were  measured  at  33  Hz  increments.  This  was  done  for  0.5  V,  1.0  V,  2.0  V,  4.0  V,  and 
8.0  V.  For  each  of  the  amplitude  response  measurements,  amplitudes  were  swept 
logarithmically  from  0.03  to  8.3  volts  and  then  back  down  through  the  same  amplitudes 
for  a  total  of  120  measurements.  This  was  done  for  99  Hz,  198  Hz,  396  Hz,  792  Hz,  and 
1584  Hz.  On  these  graphs  the  upward  sweep  was  plotted  with  a  solid  line  and  the 
downward  sweep  was  plotted  with  a  dotted  line  in  order  to  check  for  any  hysteresis 
effect.  The  shaker  used  in  Experiment  Two  was  selected  based  on  these  tests.  At  that 
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time  the  input  voltage  was  related  to  the  shaker  foot  force  as  described  in  Chapter  VI  - 
Experiment  Two. 

201b  Shaker  (Rectangular  Foot') 

The  four  accelerometers  were  arranged  in  the  following  fashion  for  the  test  of  this 
foot.  The  three  PCB  352B22  accelerometers  were  placed  on  the  bottom  of  the  foot  at  the 
left  edge,  center,  and  right  edge.  This  placement  looked  just  like  Figure  3.3  for  the  100  lb 
shaker  and  large  rectangular  foot.  Their  sensitivities  were  10.6  mV/g,  9.3  mV/g,  and  10.1 
mV/g  respectively.  The  PCB  352C67  accelerometer  was  placed  on  the  side  of  the  foot, 
in  the  center,  to  capture  any  horizontal  motion.  Its  sensitivity  is  approximately  1 1  times 
(109.5  mV/g)  that  of  the  PCB  352B22  accelerometers. 

Figure  5.1  shows  the  unloaded  response  of  the  20  lb  shaker  with  rectangular  foot. 
For  the  unloaded  test  in  air  the  motion  of  the  center  of  the  foot  increased  linearly  with 
increasing  amplitude.  As  (a)  shows  however,  the  amplitude  was  substantially  reduced  at 
1584  Hz.  This  was  indicative  of  a  null  in  the  frequency  response  located  around  1575  Hz 
as  shown  in  (b).  This  figure  also  shows  a  resonance  centered  at  1785  Hz  for  the  motion 
in  the  center  of  the  foot.  The  harmonics  remain  in  the  noise  floor  (over  40  dB  less  than 
the  fundamental)  for  the  entire  range  of  amplitudes,  which  indicates  a  linear  behavior 
when  the  foot  is  unloaded.  The  ends  of  the  foot  are  not  experiencing  the  same  null  at 
1575  Hz  however.  This  is  indicated  in  Figure  (c)  as  the  left  edge  normalized  by  the 
center  is  30  to  45  times  greater  for  1584  Hz  than  it  is  for  the  other  frequencies.  The  foot 
is  fimctioning  as  a  dynamic  vibration  absorber  when  it  resonates  in  this  mode.  A 
comparison  of  the  left  and  right  edges  normalized  by  the  center  is  shown  in  (d).  There  is 
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Figure  5.1  -  201b  Shaker  with 
Rectangular  Foot  in  Air  (a)  Amplitude 
response  of  center  for  5  frequencies  (b) 
Frequency  response  of  center  for  5 
amplitudes  (c)  Amplitude  response  of  the 
left  edge  normalized  by  the  center  (d) 
Comparison  of  the  left  and  right  edges 
normalized  by  the  center  for  2.0  V  (e) 
Phase  of  the  left  and  right  edges 
normalized  by  the  center  for  2.0  V 
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a  symmetric,  relative  resonance  at  1575  Hz  with  an  associated  phase  shift  shown  in  (e). 
The  horizontal  accelerometer  did  not  register  anything  within  20  dB  of  the  values  of  the 
vertical  motion. 

Figure  5.2  shows  what  happened  when  the  sand  loaded  resonances  were 
examined.  As  expected,  the  loaded  response  of  the  shaker  in  sand  was  much  different. 

The  frequency  response  for  the  motion  in  the  center  of  the  foot  changed  as  shown  in  (a). 

It  was  no  longer  simply  related  to  the  drive  level  due  to  several  new  resonances.  The 
fundamental  and  two  harmonics  for  the  center  accelerometer  and  left  edge  accelerometer 
are  shown  in  (b)  and  (c)  respectively.  The  motion  became  nonlinear  as  indicated  by  the 
number  of  harmonics  that  were  produced  in  the  loaded  situation.  Just  how  complicated 
this  foot  motion  became  when  the  shaker  was  placed  on  the  sand  is  shown  in  (d).  Not 
only  did  the  shaker  exhibit  several  new  resonances  in  the  frequency  response,  but  the  left 
edge  could  vibrate  more  or  less  than  the  center  depending  on  which  frequency  is 
examined.  To  complicate  matters  further,  (e)  shows  that  the  flapping  of  the  ends  is  not 
symmetrical  across  the  foot.  Vertical  displacements,  bending  about  the  center  (in  phase 
end  displacements),  and  rocking  about  the  center  (out  of  phase  end  displacements)  were 
mode  shapes  that  all  appeared  to  be  present. 

Finally,  the  amplitude  response  data  was  taken  with  the  radar.  These 
measurements  were  taken  along  the  x  axis  at  x  =  10  cm  and  x  =  40  cm  on  the  positioner 
(origin  26  cm  from  shaker  foot).  This  data,  shown  in  Figure  5.3,  begins  at  amplitudes 
that  are  down  in  the  noise  floor.  They  rise  above  the  floor  once  the  input  voltage  reaches 
approximately  0.1  V.  The  curves  appear  to  rise  linearly  although  there  is  more  deviance 
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Figure  5.2  -  201b  Shaker  with 
Rectangular  Foot  on  Sand  (a) 
Frequency  response  of  center  for  5 
amplitudes  (b)  Frequency  response  of 
the  center  showing  fundamental  and 
first  two  harmonics  at  8.0  V  (c) 
Frequency  response  of  the  left  edge 
showing  fundamental  and  first  two 
harmonics  at  8.0  V  (d)  Amplitude 
response  of  the  left  edge  normalized  by 
the  center  for  5  frequencies  (e) 
Comparison  of  the  left  and  right  edges 
normalized  by  the  center  for  2.0  V 
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from  a  straight  line  at  x  =  40  cm  than  at  x  =  10  cm.  This  was  expected  since  the 
nonlinearities  of  the  sand  would  have  more  effect  on  the  signal  the  further  it  propagated 
in  the  sandbox.  Saturation  was  not  seen  on  the  high  end  of  the  amplitudes.  Dynamic 
fluidization  limited  drive  levels  before  signal  saturation  occurred. 

To  summarize,  the  small  foot  used  in  Experiment  One  had  extremely  complicated 
motion.  It  exhibited  several  resonances  and  harmonics  were  generated  at  the  foot.  It  did 
not  move  up  and  down  rigidly  as  assumed,  but  instead  flexed  in  a  nonsymmetrical  way 
while  also  rocking  back  and  forth.  These  factors  combined  to  make  the  20  pound  shaker 
with  small  rectangular  foot  a  poor  choice  for  Experiment  Two. 

201b  Shaker  TCircular  Footl 

The  four  accelerometers  were  arranged  in  the  following  fashion  for  the  test  of  this 
foot:  three  PCB  352B22  accelerometers  were  equally  spaced  on  the  bottom  of  the  foot 
around  the  edge,  and  the  PCB  352C67  accelerometer  was  placed  on  the  bottom  in  the 
center.  The  same  three  tests  done  on  the  20  lb  shaker  with  small  rectangular  foot  were 
again  done  for  this  foot. 

The  unloaded  test  produced  results  similar  to  those  for  the  rectangular  foot.  Once 
again  the  1584  Hz  frequency  did  not  increase  as  fast  as  other  frequencies  when  the 
amplitude  increased.  Figure  5.4  shows  the  results  of  the  unloaded,  round  foot  tests.  The 
frequency  response  at  the  center  of  the  round  foot  is  shown  in  (a).  The  frequency 
response  measured  by  one  of  the  edge  accelerometers  is  shown  in  (b).  It  remains  similar 
to  the  center  frequency  response  up  to  about  1000  Hz  at  which  point  the  behavior  is  very 
different.  The  measured  resonances  of  1225  Hz  and  1575  Hz  at  the  edge  indicate  that  the 
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Figure  5.4  -  201b  Shaker  with  Round 
Foot  in  Air  (a)  Frequency  response  of 
center  for  5  amplitudes  (b)  Frequency 
response  of  edge  1  for  5  amplitudes  (c) 
Comparison  of  edge  1  and  2  normalized 
by  edge  3  for  2.0  V  (d)  Phase  of  edge  1 
and  2  normalized  by  edge  3  for  2.0  V  (e) 
Amplitude  response  of  edge  1 
normalized  by  edge  3  for  5  frequencies 
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mode  shape  for  these  frequencies  have  anti-nodes  at  the  edge.  According  to  (c)  and  (d) 
the  motion  was  not  uniform  around  the  foot.  Two  of  the  edge  accelerometer  responses 
normalized  by  the  third  edge  are  shown  in  (c).  The  curve  remains  close  to  one  for 
frequencies  less  than  1000  Hz,  which  is  also  the  range  where  the  center  motion  was 
similar  to  the  edge  motion.  Above  1000  Hz  the  edge  motion  is  not  uniform.  The  motion 
is  complicated  in  that  two  edge  accelerometer  responses  normalized  by  the  same  edge 
accelerometer  yield  very  different  results.  There  is  edge  motion  occurring  out  of  phase 
from  the  motion  at  another  point  on  the  edge  as  shown  in  (d).  This  appears  to  be  a 
saddle-shaped  mode,  although  three  evenly  spaced  accelerometers  on  the  edge  did  not 
confirm  it.  Finally,  (e)  shows  that  1584  Hz  does  not  increase,  with  an  increase  in 
amplitude,  at  the  same  rate  around  the  edge. 

The  same  tests  were  then  performed  on  the  round  foot  as  it  was  sitting  on  the 
sand.  Figure  5.5  shows  the  results  of  this  test.  The  new  frequency  response  of  the  center 
motion  is  seen  in  (a).  Once  again  there  are  several  resonances  of  this  shaker  —  shaker  foot 
configuration.  The  same  two  edges  normalized  by  the  third  edge  were  checked  again  to 
see  how  the  loaded  conditions  affected  the  complicated  edge  motion  found  in  unloaded 
conditions.  Once  again  the  motion  at  one  part  of  the  edge  is  very  different  from  motion 
at  the  other  parts  of  the  edge  as  indicated  in  (b)  and  (c).  Not  only  is  the  foot  bending, 
there  are  frequencies  for  which  the  bending  is  in  phase  and  frequencies  for  which  the 
bending  is  out  of  phase.  The  response  of  this  foot  as  the  amplitude  increased  was  also 
examined  by  plotting  an  amplitude  response  of  an  edge  normalized  by  an  edge.  This  is 
shown  in  (d).  Although  the  motion  of  the  foot  is  very  complicated  due  to  the  difference 
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in  relative  magnitudes  around  the  edge,  it  is  made  even  more  complicated  because  this 
relativity  changes  with  increasing  amplitude.  Finally,  the  harmonics  being  produced  by 
the  foot  are  shown  in  (e)  for  the  center  accelerometer  response.  There  are  a  significant 
number  of  harmonics  manifested  above  the  noise  floor  throughout  the  range  of 
fi-equencies.  This  indicates  that  the  circular  foot  motion  is  nonlinear. 

The  round  shaker  foot  was  then  used  to  check  amplitude  response  for  396  Hz  at  x 
=  10  cm  and  x  =  40  cm  on  the  x  axis  of  the  sandbox.  These  curves  are  shown  in  Figure 
5.6.  The  curves  rise  out  of  the  noise  floor  once  the  amplitude  reaches  0. 1  volts.  At  x  = 

10  cm  the  saturation  curve  is  rising  linearly  as  expected.  Instead  of  beginning  to  saturate 
however,  the  curve  rises  rapidly  at  2.5  volts.  The  slope  of  this  curve  then  decreases 
further  up  in  amplitude.  The  curve  at  x  =  40  cm  is  very  nonlinear.  It  appears  to  reach 
saturation  prior  to  one  volt  but  then  increases  again.  In  both  of  the  saturation 
measurements  hysteresis  is  evident  as  the  amplitudes  sweep  back  down.  This  may  be  a 
result  of  front  edge  versus  back  edge  arrivals  but  when  the  amplitudes  come  back  down 
the  curves  become  straighter  which  would  indicate  that  some  packing  of  the  sand  had 
occurred. 

To  summarize  the  results  for  the  20  pound  shaker  and  round  foot,  the  behavior 
was  extremely  complicated  and  nonlinear.  Several  resonances  were  present  under  loaded 
conditions.  The  foot  began  bending  for  fi’equencies  above  1000  Hz.  This  bending  was 
not  symmetrical  around  the  edges.  It  occurred  with  different  magnitude  and  changed  in 
and  out  of  phase  as  the  amplitude  increased.  The  motion  was  too  complicated  for  three 
accelerometers  on  the  edge  to  accurately  document  the  different  mode  shapes  present. 
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The  one  advantage  that  the  round  foot  had  over  the  rectangular  foot  was  its  increased 
surface  area  which  could  facilitate  longer  measurements  before  experiencing  dynamic 
fluidization  of  the  sand. 


1001b  Shaker  tRectaneular  FooO 

The  four  accelerometers  were  arranged  in  a  similar  manner  for  the  test  of  this  foot 
as  they  were  for  the  rectangular  foot  used  on  the  20  pound  shaker  (see  Figure  3.3).  The 
three  PCB  352B22  accelerometers  were  placed  on  the  bottom  of  the  foot  at  the  left  edge, 
center,  and  right  edge.  Their  sensitivities  were  10.6  mV/g,  9.3  mV/g,  and  10.1  mV/g 
respectively.  The  PCB  352C67  was  once  again  placed  on  the  side  of  the  foot,  in  the 
center,  to  capture  any  horizontal  motion. 

The  first  test,  which  was  done  unloaded,  yielded  results  similar  to  those  for  the  20 
pound  shaker  with  rectangular  foot.  Figure  5.7  shows  these  results.  The  frequency 
response  of  the  center  accelerometer  is  shown  in  (a)  and  the  amplitude  response  is  shown 
in  (b).  Just  as  in  the  20  pound  shaker  test,  this  rectangular  foot  exhibited  a  null  and  one 
resonance.  Both  of  these  occurred  at  lower  frequencies  than  in  the  20  pound  shaker  case 
however.  The  null  for  the  100  pound  shaker  test  in  air  is  at  1250  Hz  and  the  resonance  is 
centered  around  1485  Hz.  This  resonance  shows  up  in  the  1584  Hz  curve  of  (b).  The  left 
edge  normalized  by  the  center  amplitude  response  is  shown  in  (c).  A  comparison  of  the 
left  and  right  edge  acceleration  normalized  by  the  center  is  shown  in  (d).  The  motion  is 
very  uniform  across  the  shaker  foot  with  the  exception  of  a  resonance  around  1250  Hz 
due  to  a  dynamic  vibration  absorber  effect.  A 180  degree  phase  shift  at  this  same 
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Figure  5.7  -  1001b  Shaker  with 
Rectangular  Foot  in  Air  (a)  Frequency 
response  of  center  for  5  amplitudes  (b) 
Amplitude  response  of  center  for  5 
frequencies  (c)  Amplitude  response  of  left 
edge  normalized  by  the  center  for  5 
frequencies  (d)  Comparison  of  left  and 
right  edge  normalized  by  the  center  for 
2.0  V  (e)  Phase  of  left  and  right  edge 
normalized  by  the  centerfor  2.0  V 
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frequency  is  apparent  in  (e).  According  to  (d)  and  (e)  however,  the  motion  of  the  two 
ends,  while  different  than  the  center,  is  almost  identical  with  respect  to  each  other. 

Figure  5.8  contains  the  results  of  the  sand  loaded  test  of  the  100  pound  shaker. 
When  it  was  tested  in  the  sand  its  frequency  response  exhibited  fewer  resonances  than  the 
20  pound  shaker  with  rectangular  foot.  The  frequency  response  for  the  center  motion  is 
shown  in  (a).  The  frequency  response  showed  a  null  around  100  Hz  but  this  was  not  a 
repeatable  result  according  to  all  of  the  measurements  taken  during  Experiment  Two. 
There  was  a  resonance  centered  around  1585  Hz.  The  motion  of  this  foot  was  nonlinear 
as  indicated  by  the  harmonics  shown  in  (b).  If  these  curves  are  compared  with  those  for 
the  small  rectangular  foot,  one  can  see  that  the  level  of  harmonics  being  produced  in  the 
100  pound  shaker  arrangement  is  not  as  significant  as  for  the  20  pound  shaker 
arrangement.  The  relative  motion  of  the  ends  is  shown  in  (c)  and  (d).  This  indicates  that 
although  there  is  still  some  flexing  of  the  left  and  right  ends,  it  is  symmetrical  across  the 
length  of  the  foot.  This  motion  was  less  complex  than  the  motion  of  the  20  pound  shaker 
with  rectangular  foot.  Unlike  the  20  pound  shaker  however,  the  100  pound  shaker 
exhibited  significant  motion  in  the  horizontal  direction.  The  center  frequency  response 
normalized  by  the  horizontal  response  is  shown  in  (e).  Seventy  percent  of  this  frequency 
range  contains  magnitudes  in  the  horizontal  direction  that  are  within  20  dB  of  the  vertical 
magnitudes. 

Finally,  the  amplitude  response  was  measured  at  x  =  10  cm  and  x  =  40  cm  with 
the  radar.  This  is  shown  in  Figure  5.9.  In  addition  to  producing  larger  displacements  in 
the  sand,  this  shaker  -  shaker  foot  combination  produced  a  more  linear  amplitude 


45 


Vortical/Horlzontal  Acceleration  Edge/Center  Acceleration  Acceleration  (m/s^) 


Frequency  (Hz) 


46 


response.  The  measurements  for  396  Hz  at  x  =  10  cm  and  x  =  40  cm,  as  shown  in  (a)  and 
(b)  respectively,  begin  in  an  obvious  linear  region  and  progress  to  what  appears  to  be  the 
beginning  of  saturation.  The  curves  also  have  less  of  a  hysteresis  effect  than  do  the 
previous  arrangements. 

The  100  pound  shaker  with  rectangular  foot  was  chosen  for  Experiment  Two. 

The  reason  for  this  decision  was  threefold.  First,  this  arrangement  produced  the  least 
amount  of  nonlinearities  at  the  source.  Fewer  harmonics  were  produced  in  the  foot,  and 
the  foot  motion  was  more  rigid  than  the  other  two  arrangements  tested.  Second,  the 
rectangular  foot  on  the  100  pound  shaker  had  sufficient  surface  area  to  support  the  shaker 
during  extended  experiments.  This  would  allow  for  taking  complete  sets  of  data  without 
moving  the  shaker  or  burying  the  shaker  foot.  Third,  the  100  pound  shaker  and  large 
rectangular  foot  had  fewer  resonances  in  the  sand  loaded  condition.  Ideally,  a  foot  with 
no  resonances  would  have  been  used  but  time  did  not  permit  designing  and  testing 
another  foot. 
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CHAPTER  VI 


EXPERIMENT  TWO 

Procedures 

This  section  includes  two  major  topics.  The  first  is  how  the  experiments  were 
designed.  It  describes  what  the  techniques  and  methods  for  gathering  data  were  and  why 
they  were  chosen  over  others.  The  second  major  topic  discusses  the  details  of  each 
individual  test  run  in  the  laboratory.  It  documents  the  conditions  that  were  unique  to  each 
test  so  that  they  may  be  considered  during  the  evaluation  of  the  data. 

Design  of  Experiments 

The  objective  of  Experiment  Two  was  to  generate  frequency  and  amplitude 
responses  for  the  sand,  correcting  the  mistakes  from  Experiment  One.  In  particular,  the 
input  voltage  was  related  to  the  shaker  foot  force,  the  data  collection  software  was  fixed 
to  prevent  leakage  from  influencing  the  noise  floor,  the  range  of  amplitudes  for  the 
saturation  curves  were  increased,  and  the  tests  were  done  without  moving  the  shaker. 

The  first  problem  that  needed  to  be  fixed  was  the  fact  that  the  input  voltages  from 
Experiment  One  were  meaningless  without  knowing  what  the  power  amplifier  was  doing. 
This  was  fixed  by  measuring  the  current  between  the  power  amplifier  and  the  shaker  with 
two  different  current  probes  to  check  accuracy.  The  force  of  the  shaker  foot  was  related 
to  the  current  into  the  shaker  by  a  10  Ibs/amp  approximation  given  in  the  specification 
sheet  for  the  1001b  shaker.  A  plot  of  the  frequency  response  for  shaker  force  per  input 
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voltage  is  shown  in  Figure  6. 1 .  The  two  current  probes  measured  the  same  thing  so 
Figure  6. 1  shows  the  data  from  just  one  of  them.  Five  different  amplitudes  were  tested 
(0.5  V,  1.0  V,  2.0  V,  4.0  V,  and  8.0  V)  and  each  time  the  shaker  force  per  input  voltage 
curve  was  the  same.  All  of  the  amplitude  response  tests  and  the  first  two  frequency 
response  tests  of  Experiment  Two  were  done  with  the  power  amplifier  on  Gain  Setting  1. 
The  third  and  fourth  frequency  response  tests  of  Experiment  Two  were  done  with  the 
power  amplifier  on  Gain  Setting  2.  Gain  Setting  1  is  approximately  7  dB  greater  than 
Gain  Setting  2. 

The  second  problem  dealing  with  the  data  acquisition  software  was  fixed  by  using 
buffers.  In  Experiment  One  a  4.096  second  input  was  used,  approximately  three  and  a 
half  seconds  of  which  was  a  continuous  wave  signal  and  approximately  half  second  of 
which  was  settling  time.  There  were  65536  points  recorded,  so  there  were  16,384  points 
in  1.024  seconds  of  the  signal.  The  first  and  last  16,384  points  were  windowed  out  to 
eliminate  the  ring-up  and  ring-down  transients.  The  middle  32,768  points  were  used  for 
the  FFT.  These  points  comprised  exactly  2.048  seconds.  When  a  33  Hz  signal  was 
generated,  67.584  cycles  of  this  waveform  were  contained  in  the  middle  32,768  points 
saved.  Neither  the  beginning  nor  the  end  coincided  with  the  point  between  two  cycles. 
This  held  true  for  frequencies  other  than  33  Hz  also. 

For  Experiment  Two  the  use  of  buffers  guaranteed  an  integer  number  of  cycles. 
The  Lab  VIEW  program  created  buffers  containing  2048  points  in  each.  This  number 
always  remained  the  same.  When  1 1  Hz  was  generated  there  was  one  cycle  of  this 
waveform  in  each  buffer  (2048  pts/cycle).  When  22  Hz  was  generated  there  were  two 
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Figure  6.2  Buffer  Technique  of  Taking  Data 
(22  Hz  shown  in  figure) 
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cycles  per  buffer  (1024  pts/cycle).  There  were  three  cycles  per  buffer  for  33  Hz  and  so 
on  so  that  there  was  always  2048  points  per  buffer  and  every  buffer  contained  an  integer 
number  of  cycles.  The  signal  generation  rate  used  was  22522.5234  pts/sec  so  each  buffer 
was  0.09093  second  long.  This  meant  that  the  actual  frequencies  used  did  not  match  the 
requested  frequencies  exactly.  For  example,  when  99  Hz,  198  Hz,  396  Hz,  792  Hz,  and 
1584  Hz  were  requested  for  the  amplitude  response  tests,  the  actual  frequencies  used 
were  98.98  Hz,  197.95  Hz,  395.90  Hz,  791.81  Hz,  and  1583.61  Hz.  The  rounded  off 
frequencies  are  used  throughout  the  discussion  of  the  results. 

Once  Lab  VIEW  generated  one  buffer  based  on  the  requested  frequency,  19  more 
buffers  were  generated  as  shown  in  Figure  6.2.  The  first  three  buffers  and  the  last  buffer 
were  windowed  out  to  eliminate  any  ring-up  and  ring-down  transient  effects.  The 
remaining  16  buffers  were  used  to  analyze  the  results.  This  ensured  that  the  number  of 
points  used  for  the  FFT  was  always  a  power  of  two  (32,768  in  this  case)  and  the  number 
of  cycles  was  an  integer  value.  The  highest  frequency  requested  was  2002  Hz.  This 
frequency  had  1 1  points  per  cycle  which  was  enough  to  prevent  aliasing. 

The  third  problem  from  Experiment  One  that  needed  to  be  corrected  was  the 
range  of  amplitudes  used  to  generate  the  amplitude  response.  Previously,  the  amplitudes 
were  not  low  enough  to  ensure  a  beginning  in  the  linear  region  of  the  curve.  Experiment 
Two  would  take  advantage  of  the  widest  range  of  amplitudes  possible.  The  data 
acquisition  card  utilized  would  limit  this.  The  maximum  voltage  that  could  be  input  from 
the  Lab  VIEW  program  was  10  volts.  Based  on  experience,  the  minimum  input  voltage 
that  registered  above  the  noise  floor  was  somewhere  around  0.025  volts.  The  starting 
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voltage  was  selected  as  0.03  volts.  The  previous  range  of  amplitudes  for  Experiment  One 
was  a  difference  of  24  dB  from  the  low  end  to  the  high  end.  By  starting  at  0.03  volts,  and 
increasing  logarithmically,  60  measurements  could  be  made  covering  about  a  50  dB 
dynamic  range.  The  highest  voltage  would  be  about  8.3  volts. 

The  fourth  problem  was  selecting  a  gain  setting  on  the  power  amplifier  that  would 
allow  8.3  volts  to  come  in  but  the  current  sent  to  the  shaker  would  not  drive  it  at  such  a 
level  that  the  foot  buried  into  the  sand.  Several  sample  amplitude  responses  were 
generated  in  order  to  find  this  setting  using  the  0.03  to  8.3  volt  range  selected  for  the 
card.  Once  the  setting  on  the  amplifier  was  selected  several  things  were  checked.  First, 
the  beginning  of  the  measurements  had  to  be  in  the  linear  regime  of  the  amplitude 
response.  Second,  the  power  amplifier  setting  had  to  be  high  enough  so  that  8.3  volts 
would  result  in  the  curve  showing  signs  of  saturation.  Finally,  the  setting  on  the  power 
amplifier  had  to  be  such  that  dynamic  fluidization  did  not  occur  during  a  24  hour  test. 

The  combination  of  voltages  and  the  gain  settings  seen  in  Figure  6.1  allowed  all  of  these 
criteria  to  be  met. 

When  preparing  for  Experiment  Two,  lessons  from  Experiment  One  and  from 
pre-experiment  tests  were  applied  to  the  design  of  the  experiment.  This  resulted  in  four 
additional  changes  that  were  incorporated  into  Experiment  Two.  These  changes  were 
replacing  the  calibration  runs  (discussed  in  Appendix  A  —  Experiment  One  Details)  with 
a  four-accelerometer  test,  the  number  of  tests  done  at  one  time,  how  often  the  sand  had  to 
be  reconditioned,  and  the  duration  of  the  incident  signal. 


53 


In  Experiment  One  a  calibration  scan  was  done  between  each  measurement  so 
that  after  the  data  was  processed  there  was  data  taken  at  the  same  position,  amplitude, 
and  frequency  for  each  measurement  to  compare  to  each  other.  It  was  realized  during 
Experiment  One  however  that  the  frequency  responses  and  saturation  curves  being 
generated  were  actually  repeatable  as  long  as  dynamic  fluidization  did  not  occur  and  the 
shaker  was  not  lifted  up  and  placed  back  down.  Because  of  this,  it  was  determined  that 
the  really  important  information  was  what  kind  of  foot  motion  was  being  generated  for 
that  particular  test  due  to  that  unique  shaker  foot  to  sand  contact. 

To  answer  this  a  test  was  run  before  every  frequency  response  or  amplitude 
response  data  group  was  taken.  The  test  recorded  the  responses  of  four  accelerometers 
placed  in  the  same  configuration  described  in  Chapter  V.  Prior  to  the  frequency  response 
measurements,  the  four  accelerometer  test  recorded  60  frequencies  (33  Hz  to  1980  Hz  at 
33  Hz  increments)  at  five  different  amplitudes  (0.5  V,  1.0  V,  2.0  V,  4.0  V,  and  8.0  V). 
Prior  to  the  amplitude  response  measurements,  the  four  accelerometer  test  recorded  60 
amplitudes  (0.03  V  to  8.3  V  increasing  logarithmically)  at  five  different  frequencies  (99 
Hz,  198  Hz,  396  Hz,  792  Hz,  and  1584  Hz).  In  addition  to  this,  two  accelerometers  were 
recorded  at  every  point  throughout  both  of  the  measurements.  By  doing  this,  for  every 
piece  of  information  collected,  one  could  go  back  to  see  what  the  shaker  foot  was  doing. 
It  could  then  be  determined  how  much  of  the  results  could  be  attributed  to  the  motion  of 
the  shaker  foot  and  how  much  could  be  attributed  to  the  propagation  path  in  the  sand. 

One  of  the  specific  things  that  the  calibration  tests  of  Experiment  One  were 
designed  for  was  to  check  how  much  the  sand  drying  affected  the  data  being  collected. 
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Now  that  a  shaker,  power  amplifier  setting,  and  input  voltage  had  all  been  selected  so  that 
dynamic  fluidization  did  not  occur,  the  method  for  checking  drying  effects  was  changed. 
For  Experiment  Two,  once  a  frequency  response  measurement  was  taken  the  procedure 
was  repeated.  This  allowed  for  a  comparison  of  the  results  with  nothing  changing  but  the 
moisture  content  in  the  sand.  Every  amplitude  response  set  was  also  done  twice  without 
moving  anything.  Another  feature  was  also  added  to  the  amplitude  response 
measurements.  Instead  of  sweeping  up  in  amplitude  and  then  going  to  the  next 
frequency,  the  test  swept  up  in  amplitude  and  then  back  down  the  same  way.  This 
provided  information  for  any  hysteresis  that  might  be  present  and  pinpointed  places 
where  the  particular  frequency,  amplitude,  and  time  had  caused  the  shaker  foot  to  sand 
contact  to  change  significantly. 

Prior  to  Experiment  One  it  was  thought  that  ten  hours  was  the  maximum  duration 
that  tests  could  be  run  before  halting  to  rewet  and  repack  the  sand.  After  several  weeks 
of  collecting  data  it  was  found  that  the  conditions  in  the  sandbox  remained  almost  the 
same  for  much  longer.  After  wetting  and  packing  the  sand  prior  to  a  test  the  factors 
affecting  the  wave  propagation  change  the  most  during  the  first  hour.  It  is  during  this 
time  that  the  moisture  in  the  sand  settles  into  some  quasi-equilibrium  state.  After  two 
hours  very  little  change  occurs  for  the  next  36  to  48  hours.  The  goal  for  Experiment  Two 
then  was  to  be  able  to  get  two  frequency  response  or  two  amplitude  response 
measurements  done  within  24  hours. 

The  last  big  change  from  Experiment  One  came  after  the  buffers  had  been  used  to 
run  some  sample  experiments.  Because  the  buffers  fixed  the  leakage  problem  the  noise 
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floor  remained  at  a  constant  level.  As  the  amplitudes  were  swept  up  an  excellent  signal 
to  noise  ratio  was  achieved.  Because  of  this  improved  ratio  many  harmonics  that  were 
lost  in  noise  previously  were  now  seen.  For  example,  at  x  =  10  cm  (origin  was  26  cm 
from  the  shaker  foot)  for  4.0  V,  the  first  two  harmonics  were  approximately  20  dB  above 
the  noise  floor  in  the  100  -  600  Hz  band.  Depending  on  the  frequency,  amplitude,  and 
position  in  the  sandbox,  up  to  a  dozen  harmonics  could  be  discerned  in  the  frequency 
spectrum.  In  Experiment  One,  two  harmonics  were  saved,  but  for  Experiment  Two,  five 
harmonics  would  be  saved  because  five  harmonics  could  often  be  seen  above  the  noise 
floor.  Another  advantage  of  the  improved  signal  to  noise  ratio  was  that  the  input  signal 
did  not  have  to  be  as  long  as  it  was  for  Experiment  One.  The  signal  was  reduced  in 
length  from  3.6  seconds  to  1.45  seconds  and  the  same  half  second  settling  time  was  left  at 
the  end.  When  the  accelerometer  data  was  being  taken  the  signal  was  reduced  to  0.36 
seconds.  This  saved  a  great  deal  of  measurement  time. 

The  final  experimental  design  involved  two  major  tests.  One  of  the  lessons 
learned  from  Experiment  One  was  that  for  plotting  frequency  response,  24  amplitudes 
was  much  more  detail  than  necessary.  Similarly,  for  plotting  amplitude  response,  180 
frequencies  was  far  more  than  necessary.  This  prompted  the  use  of  two  separate  tests  as 
opposed  to  the  plan  for  Experiment  One  which  was  to  take  data  dense  in  frequencies  and 
dense  in  amplitudes  and  use  the  same  set  of  data  to  plot  either  frequency  or  amplitude 
response.  By  breaking  it  into  two  separate  tests  the  necessary  information  was  captured 
while  recording  less  than  30%  of  the  information  collected  in  Experiment  One. 
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The  first  major  test  was  to  measure  the  frequency  response.  This  test  lasted 
approximately  25.5  hours.  It  began  with  the  four  accelerometer  measurements  already 
discussed.  From  there,  two  complete  frequency  response  sets  were  taken.  Each  set 
began  at  33  Hz,  0.5  V,  and  x  =  10  cm.  For  all  of  the  measurements  mentioned  here  on, 
the  origin  (x  =  0)  was  26  cm  from  the  shaker  foot.  The  frequency  then  increased  by  1 1 
Hz  increments  up  to  2002  Hz.  The  amplitude  then  increased  to  1.0  V  and  the  frequencies 
were  swept  through  again.  The  amplitude  was  increased  to  2.0  V,  4.0  V,  and  finally  8.0 
V  with  180  frequencies  checked  at  each  drive  level.  All  of  this  was  done  at  five  different 
locations  in  the  sandbox.  These  locations  were  x  =  10  cm,  20  cm,  40  cm,  80  cm,  and  160 
cm.  Throughout  the  entire  process  the  response  from  two  of  the  accelerometers  was 
recorded  in  addition  to  the  radar  reading  in  the  sandbox.  This  constituted  one  complete 
set  of  frequency  response  measurements.  As  soon  as  one  was  complete  the  entire 
procedure  was  repeated  for  a  second  set. 

The  second  major  test  was  to  measure  the  amplitude  response.  This  test  lasted 
approximately  17.25  hours.  It  began  with  the  four  accelerometer  measurements  just  as 
the  frequency  response  data  did.  From  there,  two  complete  amplitude  response  sets  were 
taken.  Each  set  began  at  0.03  V,  99  Hz,  and  x  =  10  cm.  The  amplitude  then  increased  by 
approximately  0.83  dB  60  times  up  to  8.3  V.  The  amplitudes  then  follow  the  same 
sequence  coming  back  down.  After  this,  the  frequency  increased  to  198  Hz  and  the 
amplitudes  were  swept  through  again.  The  frequency  continued  to  increase  to  396  Hz, 
792  Hz,  and  1584  Hz  as  the  120  amplitudes  were  measured  each  time.  All  of  this  was 
done  at  five  different  locations  in  the  sandbox.  These  locations  were  the  same  as  for  the 
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frequency  response  data  sets.  Once  again  the  response  from  two  of  the  accelerometers 
was  recorded,  in  addition  to  the  radar  reading  in  the  sandbox,  throughout  the  entire 
process.  This  constituted  one  complete  set  of  saturation  curve  measurements.  As  soon  as 
one  was  complete  the  entire  procedure  was  repeated  for  a  second  set. 

Data  Collection 

A  total  of  four  frequency  response  tests  (8  sets)  and  five  amplitude  response  tests 
(10  sets)  were  conducted.  The  general  procedure  for  any  given  test  was  the  same.  The 
sand  was  completely  saturated  with  water.  The  actual  water  table  remained  46  to  50  cm 
below  the  surface  of  the  sandbox.  After  the  sand  was  watered  down  it  was  packed  with  a 
hand  tamper  and  allowed  to  sit  for  a  minimum  of  two  hours  before  the  data  was  collected. 
Normally  after  about  one  hour  the  surface  was  given  another  light  mist,  repacked  and  left 
alone  for  another  two  to  three  hours.  Once  this  was  done  the  shaker  was  put  on  the  sand 
and  data  collection  commenced. 

Several  things  were  checked  at  the  beginning  of  the  measurements  for  a  relative 
comparison  of  the  conditions.  These  things  included  position  of  the  radar  waveguide  and 
the  radar  power  reading.  Table  6. 1  summarizes  the  pertinent  data  for  all  of  the  frequency 
response  tests  and  Table  6.2  summarizes  the  pertinent  data  for  the  amplitude  response 
tests.  Each  of  the  power  readings  at  the  origin  found  in  Table  6. 1  and  6.2  were  +/-  0.01 
dBm.  The  power  reading  for  any  given  frequency  response  or  amplitude  response  test 
remained  within  a  3  dB  range  throughout  the  entire  test.  Also,  lower  amplitudes  were 
used  for  the  fourth  frequency  response  test.  These  amplitudes  were  0.015625  V,  0.03 125 
V,  0.0625  V,  0.125  V,  and  0.25  V. 
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Test 

1 

2 

3 

4 

Date  Started 

7  Feb  00 

11  Feb  00 

13  Feb  00 

Time  Soaked 
and  Packed 

1030 

1230 

0840 

0350 

Time  Misted 
and  Repacked 

1240 

1300 

1930 

0550 

Start  Time 

1453 

1408 

2123 

0710 

Accelerometer 

Placement 

Center  vertical 
&  center 
horizontal 

Center  vertical 
&  buried  3  in. 
below  foot 

Center  vertical 
&  buried  5. 5 in. 
below  foot 

Center  vertical 
&  buried  5.5in. 
below  foot 

Gain  Setting 

1 

1 

2 

2 

Power 

Reading  at 
Origin 

10.17  dBm 

12.12  dBm 

8.24  dBm 

9.17  dBm 

Waveguide 
Distance  from 
Sand  (x=0) 

2  cm 

1.8  cm 

1.5  cm 

1.3  cm 

Table  6.1  -  Experimental  Procedure  Data  for  Frequency  Response  Tests 


Test 

1 

2 

3 

4 

5 

Date  Started 

25  Jan  00 

4  Feb  00 

8  Feb  00 

9  Feb  00 

10  Feb  00 

Time  Soaked 
and  Packed 

0830 

0630 

1630 

1400 

1000 

Time  Misted 
and  Repacked 

1300 

1700 

1515 

1110 

Start  Time 

1630 

1239 

1922 

1643 

1519 

Accelerometer 

Center 

Center 

Center 

Center 

Center 

Placement 

vertical  & 

vertical  & 

vertical  & 

vertical  & 

vertical  & 

center 

buried  3  in. 

buried  5.5 

buried  5.5 

buried  5.5 

horizontal 

below  foot 

in.  below 
foot 

in.  below 
foot 

in.  below 
foot 

Gain  Setting 

1 

1 

1 

1 

1 

Power 

Reading  at 
Origin 

10.04  dBm 

11.65  dBm 

9.95  dBm 

10.45  dBm 

9.55  dBm 

Waveguide 
Distance  from 
Sand  (x=0) 

2  cm 

1.8  cm 

1.8  cm 

1.5  cm 

1.5  cm 

Table  6.2  -  Experimental  Procedure  Data  for  Amplitude  Response  Tests 
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The  two  accelerometers  that  were  recorded  in  conjunction  with  the  radar 
measurements  were  of  the  following  types  and  locations.  The  first  frequency  and 
amplitude  response  tests  were  done  with  two  Kistler  accelerometers  mounted  in  the 
center  of  the  foot.  One  was  on  top  of  the  foot  to  measure  vertical  acceleration  and  the 
other  was  on  the  side  to  measure  horizontal  acceleration.  The  results,  which  will  be 
discussed  later,  prompted  the  subsequent  tests  to  be  done  with  an  accelerometer  buried 
under  the  shaker  foot.  The  second  frequency  and  amplitude  response  tests  were  done 
with  two  PCB  352C67  accelerometers.  One  was  attached  to  the  bottom  center  of  the  foot 
to  measure  vertical  acceleration  and  the  other  was  buried  in  the  sand  three  inches  below 
the  foot.  The  third  amplitude  response  measurement  was  repeated  with  the  buried 
accelerometer  5.5  inches  below  the  shaker  foot.  The  placement  of  the  accelerometers  for 
the  third  and  fourth  frequency  response  tests  and  the  fourth  and  fifth  amplitude  response 
tests  were  done  with  the  same  accelerometers  in  the  same  place.  They  were  PCB  352B22 
accelerometers,  one  of  which  was  mounted  on  the  bottom  center  of  the  foot  to  measure 
vertical  acceleration  and  the  other  was  buried  5.5  inches  below  the  shaker  foot. 

Results 

Before  the  results  of  the  frequency  and  amplitude  response  data  are  presented,  one 
should  note  where  the  noise  floor  was  for  these  experiments.  Figure  6.3  shows  the  noise 
floor  recorded  for  three  different  measurements.  The  first  was  recorded  with  Gain 
Setting  1  on  the  power  amplifier  and  a  drive  amplitude  of  8.0  V,  the  second  was  recorded 
with  Gain  Setting  2  on  the  power  amplifier  and  a  drive  amplitude  of  8.0  V,  and  the  third 
was  recorded  with  Gain  Setting  2  on  the  power  amplifier  and  a  drive  amplitude  of  0.25  V. 
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Figure  6.3  -  Noise  Floor  Measured  at  First  Position  (x  =  10  cm)  for  Fifth  Amplitude 
Tested  During  First  Iteration  (a)  Test  2  -  Gain  Setting  1,  Amplitude  =  8.0  V  (b)  Test  3  - 
Gain  Setting  2,  Amplitude  =  8 .0  V  (c)  Test  4  -  Gain  Setting  2,  Amplitude = 0. 25  V 


These  three  graphs  show  that  despite  different  drive  levels,  the  noise  floor  did  not  change. 
This  was  expected  and  showed  that  the  leakage  problem  from  Experiment  One  was 
corrected.  The  noise  floor  remained  as  shown  in  Figure  6.3  for  all  of  the  measurements 
taken  in  Experiment  Two. 

Frequency  Response  Data 

Four  complete  frequency  response  tests  were  taken.  Each  test  consisted  of  two 
iterations.  Each  iteration  included  180  frequencies  taken  at  five  drive  levels  at  five 
different  locations  as  previously  described.  Frequency  Response  Test  1  recorded  data 
that  indicated  a  coding  error  made  the  first  drive  level  0.05  V  instead  of  the  desired  0.5  V. 
For  this  reason,  the  first  frequency  response  test  is  not  used  to  describe  the  results. 

Figure  6.4  shows  the  surface  displacement  as  a  fimction  of  frequency.  The 
fundamental  frequency  is  plotted  for  five  drive  amplitudes  taken  from  Frequency 
Response  Test  2  (first  iteration).  Graphs  (a)  though  (e)  are  measurements  taken  at  the 
five  locations  in  the  sandbox.  Figure  6.5  shows  the  same  data  for  Frequency  Response 
Test  2  (second  iteration).  These  iterations  confirmed  two  things.  First,  the  two  sets  of 
data  were  similar,  as  expected,  because  the  shaker  was  not  moved  between  the  two 
iterations  of  this  test.  Second,  the  drying  of  sand  did  not  significantly  alter  results 
throughout  the  26  hours  required  to  take  all  of  the  data  shown  in  Figures  6.4  and  6.5. 
Therefore,  repeatability  of  results  may  be  achieved  during  a  26  hour  period  if  the  shaker 
is  not  moved. 

There  were  certain  characteristics  common  to  the  data  shown  in  Figures  6.4  and 
6.5.  The  waves  attenuated  as  they  propagated  in  the  sandbox.  This  was  shown  by  the 
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decreasing  surface  displacements  from  x  =  10  cm  (a)  to  those  at  x  =  160  cm  (e).  The 
amount  of  attenuation  was  frequency  dependent  however.  In  general,  higher  frequencies 
attenuated  faster  than  lower  frequencies.  The  largest  surface  displacements  occurred  in 
the  100  -  600  Hz  band,  but  at  x  =  160  cm,  the  frequencies  between  400  Hz  and  600  Hz 
had  attenuated  more  than  those  between  100  Hz  and  400  Hz.  Frequencies  above  600  Hz 
did  not  propagate  well  as  shown  by  the  decreasing  slope  above  600  Hz. 

These  figures  also  show  that  the  frequency  response  became  more  nonlinear  with 
increasing  drive  amplitude  as  expected.  By  looking  at  the  measurements  taken  at  x  =  10 
cm,  one  can  see  that  the  surface  displacement  doubled  as  the  drive  amplitude  doubled  in 
the  100  -  600  Hz  band.  Above  600  Hz  however,  there  was  a  point  at  which  doubling  the 
drive  amplitude  did  not  result  in  a  doubling  of  the  surface  displacement.  The  amplitude 
where  this  occurred  became  lower  and  lower  as  the  frequency  increased. 

Another  common  result  seen  in  Figures  6.4  and  6.5  was  that  the  frequency 
response  varied  more  as  the  drive  amplitude  increased.  This  can  best  be  seen  by  looking 
at  the  data  recorded  at  x  =  10  cm  (a).  The  smaller  amplitude  curves  are  smoother  for  a 
wider  range  of  frequencies.  As  the  amplitude  increased,  nulls  in  the  frequency  response 
appeared.  More  nulls  were  present  the  higher  the  amplitude  went.  The  amount  of 
variability  increased  during  propagation  as  seen  by  the  increasing  number  of  dips  in  the 
frequency  response  when  comparing  a  given  amplitude  in  (a)  to  those  of  (b),  (c),  (d),  and 

(e). 

Figure  6.6  shows  the  surface  displacements  versus  frequency  for  Frequency 
Response  Test  3  (first  iteration).  The  fundamentals  are  plotted  for  five  drive  amplitudes 
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taken  at  five  positions.  This  data  was  taken  with  the  power  amplifier  on  Gain  Setting  2. 
The  data  was  consistent  with  Figures  6.4  and  6.5  in  that  surface  displacements  were 
approximately  7  dB  lower  due  to  the  gain  setting,  higher  frequencies  attenuated  faster, 
and  lower  amplitudes  resulted  in  a  more  linear  frequency  response  over  a  greater  band  of 
frequencies.  The  frequencies  between  800  -  2000  Hz  have  attenuated  into  the  noise  floor 
at  X  =  160  cm  for  these  drive  amplitudes  as  shown  by  the  flat  frequency  response  in 
Figure  6.6  (e). 

Figure  6.7  shows  the  surface  displacement  versus  frequency  for  Frequency 
Response  Test  4  (first  iteration).  This  data  was  taken  with  the  power  amplifier  on  Gain 
Setting  2.  The  fundamental,  for  five  drive  amplitudes  at  five  locations,  is  shown  in  the 
graphs  of  Figure  6.7.  The  data  taken  at  x  =  10  cm  (a)  showed  that  doubling  the  drive 
amplitude  doubled  the  surface  displacement  throughout  the  frequency  band  with  the 
exception  of  frequencies  between  1200  Hz  and  1300  Hz.  The  curves  were  relatively 
smooth  with  the  exception  of  this  null  and  the  lower  frequencies  that  remained  near  the 
noise  floor  for  these  drive  amplitudes.  As  the  waves  propagated  in  the  sandbox  the  same 
increase  in  variability  appeared  in  the  higher  frequencies.  This  variability  was  indicated 
by  the  dips  in  the  frequency  response  as  seen  in  the  other  figures. 

There  were  two  results  for  Frequency  Response  Test  4  that  differed  from 
Frequency  Response  Tests  2  and  3.  First,  the  frequency  band  experiencing  nonlinear 
effects  as  the  amplitude  increased  was  different.  As  stated  earlier,  for  Frequency 
Response  Test  2,  100  Hz  -  600  Hz  was  the  band  that  continued  to  double  in  surface 
displacement  as  the  amplitude  doubled  at  x  =  10  cm  (Figure  6.4  (a)).  At  x  =  160  cm  the 


67 


Surface  Displacement  (nm)  Surfara  Displacement  (^m)  Surfa^  Displa^ment  (nm) 


—  0.015625  V 

—  0.031250  V 
.  0.062500  V 

—  0.125000  V 

—  0.250000  V 


'  '■%»  ^  '^5 

VW'<  >  . .  ■ 

!  I  •  !.  jj  ^ 


mf'N"’"' . 


V:Y'-y  'v^^f/;.-  !!!••'•' 


ijlfoM  I  'V'. '  •A'X  I  r^\n 

m  vv-v'^-vM-i 


w  At 


15' -!  ■  V'AS^;v:vm3 


500  1000  1500 

Frequency  (Hz) 


Figure  6.7  -  Frequency  Response  Test 
4  (Gain  Setting  2),  First  Iteration: 
Fundamental  Plotted  for  5 
Amplitudes  (a)  x  =  10  cm  (b)  x  =  20 
cm  (c)  X  =  40  cm  (d)  x  =  80  cm  (e)  x  = 
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frequencies  between  400  Hz  and  600  Hz  had  attenuated  more  and  were  more  variable 
than  those  between  100  Hz  and  400  Hz.  For  Frequency  Response  Test  4  however,  the 
surface  displacements  doubled  with  doubling  drive  amplitude  between  300  Hz  and  1200 
Hz  at  X  =  10  cm  (Figure  6.7  (a)).  At  x  =  160  cm ,  with  the  exception  of  frequencies  lower 
than  400  Hz  which  were  in  the  noise  floor,  the  surface  displacements  for  frequencies 
through  1200  Hz  still  doubled  as  the  drive  amplitude  doubled.  The  curves  showed 
relatively  smooth  attenuation,  without  the  increasing  variability,  for  all  five  amplitudes. 

The  second  difference  in  Frequency  Response  Test  4  was  that  surface 
displacement  increased  between  x  =  10  cm  and  x  =  20  cm  and  between  x  =  20  cm  and  x  = 
40  cm.  This  was  most  likely  due  to  a  property  of  the  sand  at  the  time  the  measurements 
were  taken.  Because  the  increase  in  surface  displacement  was  a  broadband  effect,  it  did 
not  occur  due  to  interference.  This  result  was  inconsistent  with  other  data  and 
inconsistent  with  the  expected  results.  Properties  of  the  sand  that  may  have  caused  this 
include  the  water  table  height,  a  volume  of  sand  that  had  a  different  density,  or  a  volume 
of  sand  that  retained  a  higher  moisture  content. 

Figure  6.8  (a)  shows  the  surface  displacements  versus  frequency  for  Frequency 
Response  Test  3  (first  iteration)  and  Frequency  Response  Test  4  (first  iteration).  Both 
were  measured  at  x  =  10  cm.  They  were  plotted  on  the  same  graph  in  order  to  see  how 
the  five  lower  amplitudes  step  up  into  the  five  higher  amplitudes.  It  is  important  to  note 
that  the  shaker  had  been  moved,  and  the  sand  reconditioned,  between  the  two  tests  so  the 
same  shaker  foot  to  sand  contact  was  not  present  for  all  ten  amplitudes.  This  result  can 
be  observed  by  the  nonlinear  increase  of  the  lower  frequency  range  between  0.25  V  and 
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Figure  6.8  -  Frequency  Response  Tests  3  and  4  (Gain  Setting  2),  First  Iteration: 
Fundamental  Plotted  for  10  Amplitudes  (a)  Radar  measurement  at  x  =  10  cm 
(b)  Accelerometer  measurement  of  center  taken  while  radar  was  at  x  =  1 0  cm 


0.50  V.  Figure  6.8  (b)  shows  the  frequency  response  of  the  shaker  foot,  measured  by  the 
center  accelerometer,  for  the  same  10  measurements.  The  dependence  of  the  frequency 
response  on  the  shaker  foot  to  sand  contact  was  evident  in  this  graph,  particularly  in  the  0 
-400  Hz  band. 

In  order  to  determine  the  degree  of  nonlinearities  present,  the  harmonics  produced 
were  examined.  Figure  6.9  compares  the  surface  displacements  for  the  fundamental  and 
four  harmonics  at  five  locations  ((a)  -  (e))  for  a  constant  drive  amplitude.  The  data  is 
taken  from  Frequency  Response  Test  2  (first  iteration)  with  an  amplitude  of  2.0  V.  These 
results  are  nonlinear  as  shown  by  the  harmonics  generated.  The  fundamental  and 
harmonics  attenuated  as  they  propagated  through  the  sandbox.  The  higher  harmonics 
attenuated  the  fastest,  which  showed  once  again  that  the  higher  frequencies  did  not 
propagate  as  well. 

Figure  6.10  compares  the  surface  displacements  for  the  fundamental  and  four 
harmonics  for  five  drive  amplitudes  ((a)  -  (e))  at  one  location.  The  data  is  taken  from 
Frequency  Response  Test  2  (first  iteration)  at  x  =  40  cm.  At  the  lowest  amplitude  (0.5 
V),  the  first  harmonic  was  above  the  noise  floor  for  the  100  —  1000  Hz  band  and  the 
second  harmonic  was  discemable  above  the  noise  floor  in  the  100  —  700  Hz  band.  When 
the  amplitude  was  1.0  V,  the  first  through  fourth  harmonics  were  generated  in  the  300  - 
600  Hz  band.  At  4.0  V,  at  least  one  harmonic  was  generated  throughout  the  frequency 
band  of  interest.  These  results  showed  that  the  wave  propagation  to  40  cm  was  nonlinear 
even  at  the  lowest  drive  amplitude  for  Frequency  Response  Test  2.  Appendix  B  contains 
a  complete  set  of  data  for  Frequency  Response  Test  2  (second  iteration). 
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Figure  6.9  -  Frequency  Response  Test  2 
(Gain  Setting  1),  First  Iteration; 
Fundamental  and  4  Harmonics  with 
Amplitude  =  2.0  V  (a)  x  =  1 0  cm  (b)  x  =  20 
cm  (c)  X  =  40  cm  (d)  x  =  80  cm  (e)  x  =  160 


Amplitude  Response  Data 


As  with  the  frequency  response  data,  the  amplitude  response  data  was  taken  with 
accelerometers  mounted  on  the  shaker  foot  so  that  any  surface  displacement  read  by  the 
radar  could  be  related  to  some  foot  motion.  During  Amplitude  Response  Test  1,  the  two 
accelerometers  recorded  throughout  the  measurements  were  Kistler  accelerometers. 

They  were  both  mounted  on  the  center  of  the  shaker  foot,  one  for  vertical  acceleration 
and  the  other  for  horizontal  acceleration.  After  processing  the  data,  the  curves  of  surface 
displacement  versus  amplitude  for  the  vertical  accelerometer  were  not  repeating 
themselves  between  measurements  of  the  same  frequency.  The  last  curve  was  rising  at 
about  one-third  the  rate  of  the  first  curve  for  99  Hz.  The  amplitude  response  for  higher 
frequencies  repeated  however. 

This  unexpected  effect  led  to  consultation  with  a  geophysicist.  It  was 
hypothesized  that  propagation  of  lower  frequencies  was  more  dependent  upon  the  solid 
matrix  structure  of  the  sand,  and  that  propagation  of  higher  frequencies  was  more 
dependent  on  the  viscous  forces  of  the  water  content  in  the  sand  [13].  It  was  thought  that 
as  the  measurements  were  taken,  the  foot  was  packing  the  sand  underneath  it.  If  this 
occurred,  the  sand  matrix  structure  was  changing,  thus  changing  the  amplitude  response 
of  the  shaker  foot  for  low  frequencies. 

In  order  to  confirm  that  the  sand  under  the  foot  was  being  packed  during  the 
measurements,  an  accelerometer  was  buried  three  inches  below  the  foot.  Accelerations 
measured  at  this  point  should  have  changed  over  time  as  the  volume  of  sand  effectively 
coupled  to  the  foot  increased  due  to  packing.  Two  PCB  352C67  accelerometers  were 
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used  to  take  the  measurements.  Besides  the  one  buried,  there  was  one  mounted  in  the 
center  of  the  shaker  foot.  Amplitude  Response  Test  2  was  taken  in  this  configuration  but 
did  not  record  results  consistent  with  Amplitude  Response  Test  1.  The  accelerometer 
mounted  on  the  foot  did  not  have  decreasing  amplitude  response  curves  at  99  Hz. 

Instead,  the  curves  were  consistent  throughout  the  measurements  as  originally  expected. 
The  accelerations  measured  by  the  buried  accelerometer  did  not  change  over  time  either. 

This  prompted  an  investigation  of  the  accelerometers.  The  accelerometer 
measurements  for  Amplitude  Response  Test  1  were  ruled  invalid  due  to  a  particular 
Kistler  accelerometer  used  to  record  vertical  foot  acceleration.  For  Amplitude  Response 
Tests  3, 4,  and  5,  PCB  352B22  accelerometers  were  used.  One  was  mounted  on  the 
center  of  the  foot  to  measure  vertical  acceleration  and  one  was  buried  5.5  inches  below 
the  foot  to  see  if  the  shaker  foot  was  packing  the  sand.  The  hypothesis  of  lower 
frequencies  being  dependent  on  the  sand  matrix  structure  and  higher  frequencies  being 
dependent  on  the  viscous  forces  of  the  water  was  never  confirmed  nor  denied.  Neither 
did  burying  an  accelerometer  under  the  shaker  foot  confirm  or  deny  that  the  sand  under 
the  foot  was  packed  over  time.  Amplitude  Response  Tests  3, 4,  and  5  (all  taken  on  Gain 
Setting  1)  are  used  to  present  the  results. 

Figure  6.11  shows  surface  displacement  versus  drive  amplitude  for  Amplitude 
Response  Test  3  (first  iteration).  The  fundamental  for  five  frequencies  is  plotted  in  each 
graph.  Graphs  (a)  though  (e)  are  the  measurements  taken  at  five  different  positions  in  the 
sandbox.  The  portion  of  the  curves  that  is  not  smooth  is  data  that  was  hidden  in  the  noise 
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floor.  For  example,  in  (d),  99  Hz  less  than  10°  and  1584  Hz  less  than  10  ^  are  in  the 
noise. 

Many  of  the  observations  made  for  the  frequency  response  data  were  confirmed 
with  this  data.  The  waves  attenuated  as  they  propagated  in  the  sandbox.  This  was  seen 
by  the  surface  displacement  curves,  for  a  given  frequency,  decreasing  from  one  position 
to  the  next  (from  (a)  to  (e)).  All  of  the  curves  began  in  a  linear  region.  The  curves 
entered  a  nonlinear  region,  as  expected,  when  amplitude  increased.  The  point  at  which 
these  frequencies  began  to  saturate  was  different  for  each  of  them,  showing  the  frequency 
and  amplitude  dependence  of  saturation.  The  higher  frequencies  were  less  predictable 
than  the  lower  frequencies.  For  Figure  6.11,  1584  Hz  showed  a  particularly  large  amount 
of  variability  as  it  propagated  through  the  sandbox.  The  largest  surface  displacements 
occurred  for  396  Hz,  which  was  also  consistent  with  the  frequency  response  tests. 

Figure  6.12  shows  the  surface  displacement  versus  drive  amplitude  for 
Amplitude  Response  Test  3  (second  iteration).  The  fundamental  of  five  frequencies  was 
plotted  for  five  locations  just  as  it  was  in  Figure  6. 12.  Approximately  8.5  hours  elapsed 
between  the  beginning  of  the  first  iteration  and  the  beginning  of  the  second  iteration. 
Because  the  shaker  was  not  moved  between  these  two  iterations,  the  results  were 
repeatable  as  seen  in  the  figures.  The  degree  of  repeatability  was  frequency  dependent 
however.  99  Hz,  198  Hz,  and  396  Hz,  at  all  five  locations  for  the  second  iteration,  were 
similar  to  those  of  the  first  iteration.  792  Hz  was  repeatable  for  x  =  10  cm,  20  cm,  and  40 
cm.  After  40  cm,  the  results  were  not  repeatable  in  the  nonlinear  region.  The  variability 
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of  1584  Hz  was  not  reproducible  in  the  second  iteration.  Although  it  still  varied,  the 
amplitude  where  the  variation  occurred  and  the  extent  of  the  variation  was  not  repeatable. 

Figure  6.13  shows  surface  displacement  versus  drive  amplitude  for  Amplitude 
Response  Test  4  (first  iteration).  The  fundamental  of  five  frequencies  measured  at  five 
locations  was  plotted  just  as  in  the  previous  two  figures.  Although  the  relative  magnitude 
of  surface  displacements  shown  in  Figure  6.13  were  similar  to  those  of  Amplitude 
Response  Test  3,  the  shape  of  the  curves  was  somewhat  different  due  to  the  different 
shaker  foot  to  sand  contact  that  resulted  from  moving  the  shaker  and  reconditioning  the 
sand.  The  differences  are  more  apparent  in  the  higher  frequencies  than  in  the  lower 
frequencies. 

Figure  6.14  shows  surface  displacement  versus  drive  amplitude  for  Amplitude 
Response  Test  5.  The  fundamental  of  five  frequencies  at  five  locations  was  once  again 
plotted.  A  comparison  of  Figure  6. 14  with  either  Figure  6. 1 1  or  Figure  6. 13  confirms  the 
results  stated  above.  The  results  of  Amplitude  Test  5  were  more  similar  to  those  of 
Amplitude  Test  4  however.  These  two  tests  were  done  a  day  apart  whereas  Amplitude 
Test  5  and  Amplitude  Test  3  were  done  two  days  apart.  The  changed  properties  of  the 
sand  were  more  noticeable  in  the  data  that  was  taken  two  days  apart. 

As  previously  mentioned,  data  from  two  accelerometers  was  recorded  throughout 
the  amplitude  response  tests.  Figure  6.15  is  a  side-by-side  comparison  of  the  surface 
displacement  measured  at  some  distance  in  the  sandbox  and  the  acceleration  of  the  shaker 
foot  for  that  measurement.  This  data  was  taken  from  Amplitude  Response  Test  5.  Graph 
(a)  shows  the  surface  displacement  versus  drive  amplitude  for  the  fiindamental  of  five 
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Figure  6.15-  Amplitude  Response  Test  5  (Gain  Setting  1),  Second  Iteration;  Comparison 
of  Radar  and  Center  Accelerometer  Measurements  for  5  Frequencies  (a)  Surface 
displacement  at  x  =  10  cm  (b)  Foot  acceleration  while  radar  was  at  x  =  10  cm  (c)  Surface 
displacement  atx = 40  cm  (d)  Foot  acceleration  while  radar  was  atx = 40  cm 


frequencies  measured  at  x  =  10  cm.  Graph  (b)  was  the  amplitude  response  as  measured 
in  acceleration  by  the  accelerometer  mounted  on  the  center  of  the  shaker  foot  for  the 
same  five  frequencies.  Graph  (c)  showed  the  surface  displacement  versus  amplitude 
measured  at  x  =  40  cm  and  (d)  was  the  corresponding  amplitude  response  of  the  shaker 
foot,  (b)  and  (d)  showed  that  the  amplitude  response  of  the  shaker  foot  was  almost 
identical  from  measurement  to  measurement.  They  also  showed  that  the  shaker  foot 
applied  the  most  force  at  792  Hz  and  1584  Hz.  By  the  time  the  wave  propagated  to  x  = 

10  cm  however,  792  Hz  and  1584  Hz  had  attenuated  to  the  point  that  99  Hz,  198  Hz,  and 
396  Hz  produced  the  largest  surface  displacements.  This  result  emphasized  how  much 
more  high  frequencies  attenuated  than  low  frequencies. 

Just  as  in  the  frequency  response  tests,  the  harmonics  being  produced  were 
examined  to  determine  the  extent  of  nonlinearities  present.  Instead  of  showing  the 
fundamental  and  its  harmonics  as  seen  for  the  frequency  response  tests  however,  this 
section  of  the  results  normalized  the  harmonics  by  the  fundamental  in  order  to  show 
relative  harmonic  generation. 

Figure  6.16  shows  the  first  four  harmonics  of  396  Hz  normalized  by  the  amplitude 
response  of  the  fundamental.  The  data  is  taken  from  Amplitude  Response  Test  5  (first 
iteration).  Each  graph  ((a)  -  (e))  shows  the  amplitude  response  at  a  different  location.  At 
X  =  10  cm,  the  first  harmonic  rose  above  the  noise  floor  when  the  amplitude  was  2e-l  V, 
the  second  harmonic  rose  above  the  noise  floor  when  the  amplitude  was  6e-l  V,  and  the 
third  harmonic  rose  above  the  noise  floor  when  the  amplitude  was  leO  V.  The  fourth 
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Figure  6.16-  Amplitude  Response  Test 
5  (Gain  Setting  1),  First  Iteration;  4 
Harmonics  Normalized  by  the 
Fundamental  for  396  Hz  at  5  Locations 
(a)x=  10cm(b)x=20cm(c)x  =  40cm 
(d)x=80  cm(e)x  =  160  cm 
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harmonic  just  began  to  rise  above  the  noise  at  an  amplitude  of  7e0  V.  The  harmonics 
showed  signs  of  saturation,  but  more  harmonics  were  produced,  as  amplitude  increased. 

Graphs  (b)  and  (c)  showed  that  the  harmonics  attenuated  with  respect  to  the 
fundamental  as  the  wave  propagated  in  the  sandbox.  This  was  expected  because  high 
frequencies  (harmonics  of  396  Hz)  attenuated  faster  than  low  frequencies  (396  Hz).  The 
fourth  harmonic  rose  above  the  noise  floor  at  an  amplitude  of  3e0  V  at  x  =  20  cm 
however.  This  meant  that  at  constant  amplitude  (3e0  V)  and  frequency  (396  Hz),  the 
fourth  harmonic  was  generated  at  20  cm  but  not  at  10  cm.  This  result  showed  that  the 
propagation  path  contributed  to  nonlinearity  between  x  =  10  cm  and  x  =  20  cm.  From  40 
cm  to  160  cm  this  effect  was  not  seen.  The  propagation  was  still  nonlinear,  but  the 
attenuation  of  the  harmonic  frequencies  dominated  the  effect  of  nonlinear  propagation. 

Figure  6.17  shows  four  harmonics  normalized  by  the  fundamental  versus 
amplitude.  The  data  came  from  Amplitude  Response  Test  5  (first  iteration).  All  of  the 
graphs  represent  data  measured  at  x  =  40  cm.  Each  one  was  for  a  different  frequency.  99 
Hz,  198  Hz,  and  396  Hz  all  produced  significant  harmonics  at  this  point  in  the  sand.  The 
frequency  that  had  a  harmonic  rise  above  the  noise  floor  first,  as  the  amplitude  increased, 
was  396  Hz.  The  next  frequency  to  generate  a  harmonic,  as  amplitude  increased,  was 
198  Hz.  This  result  was  consistent  with  the  frequency  response  tests  showing  the  largest 
surface  displacements  in  the  100  Hz  —  600  Hz  band.  The  fist  harmonic  of  792  Hz  rose 
above  the  noise  at  3e0  V  and  no  harmonics  were  generated  at  40  cm  for  1584  Hz.  Once 
again  this  showed  the  effects  of  attenuation  on  the  higher  frequencies.  Appendix  C 
contains  a  complete  set  of  data  for  Amplitude  Response  Test  5  (second  iteration). 
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Nonlinearities  at  the  Source 


An  important  part  of  the  results  was  what  happened  at  the  shaker  for  both  the 
frequency  response  and  amplitude  response  tests.  It  was  shown  that  both  variability  and 
nonlinearity  occur  during  propagation  but  the  source  also  contributed  to  the  overall 
effects  seen  in  the  data.  Of  the  three  shaker  -  shaker  foot  combinations  tested,  the  100 
pound  shaker  with  rectangular  foot  had  the  fewest  modes  of  foot  motion  excited  at  the 
frequencies  and  amplitudes  used,  it  had  fewer  resonances,  and  it  produced  fewer 
harmonics  as  recorded  by  the  accelerometers.  Despite  this,  there  were  still  source 
considerations  to  take  into  account. 

First,  there  was  the  matter  of  the  shaker  foot  to  sand  contact.  Every  time  the 
shaker  foot  was  placed  on  the  sand  the  foot-sand  coupling  was  different.  The  data 
collected  in  Experiment  Two  showed  that  by  leaving  the  shaker  on  the  sand  throughout 
the  measurements,  the  results  repeated  well.  The  frequency  response  changed  slightly 
whenever  the  shaker  was  moved  and  placed  back  on  the  sand. 

This  effect  is  more  noticeable  when  viewed  in  the  time-domain.  Figure  6. 18  (a) 
shows  seven  different  waveforms  from  Experiment  One  (201b  shaker  with  rectangular 
foot)  plotted  on  top  of  each  other.  They  were  all  measured  at  x  =  120  cm  but  the 
amplitudes  increased  linearly  (increment  =  0.04  V)  from  0. 14  V  to  0.38  V.  The  first  six 
amplitudes  were  recorded  without  moving  the  shaker.  The  waveforms  were  almost 
identical  with  the  exception  of  the  increased  amplitude.  Prior  to  the  seventh 
measurement  (drive  amplitude  =  0.38  V),  the  shaker  was  removed,  the  sand  was  watered 
and  packed,  and  the  shaker  was  placed  back  on  the  sand.  The  shape  of  this  waveform 
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Figure  6. 18  -  (a)  Waveform  of  7  amplitudes  (0. 14  V  -  0.38  V)  from  Experiment  One, 
201b  shaker  with  rectangular  foot,  shaker  moved  prior  to  recording  0.38  V  (dotted 
line)  (b)  Waveform  of  2.0  V  and  4.0  V  from  Frequency  Response  Test  2  (first 
iteration)  for  5  locations  (c)  Waveform  of  0.5  V  (scaled  x4)  and  8.0  V  (scaled  xO.5) 
from  Frequency  Response  Test  2  (first  iteration)  for  5  locations 
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(dotted  line)  was  significantly  different  due  to  the  new  shaker  foot  to  sand  contact.  This 
showed  that  it  was  very  important  to  understand  whether  or  not  the  data  being  compared 
in  Experiment  Two  was  taken  with  or  without  moving  the  shaker. 

Careful  inspection  of  the  waveforms  plotted  in  Figure  6.18  (a)  revealed  another 
interesting  feature.  As  the  drive  amplitude  increased,  the  waveforms  were  recorded  later 
in  time  at  the  same  point  in  the  sandbox.  The  cause  of  this  was  either  a  delay  at  the 
source,  a  decreased  propagation  speed,  or  a  combination  of  both.  Another  possibility  was 
that  as  the  drive  amplitude  increased,  the  frequency  content  due  to  harmonic  generation 
changed.  When  the  frequency  spectrum,  which  used  a  linear  assumption,  was  convolved 
with  the  differentiated  Gaussian  and  taken  into  the  time-domain  it  may  have  filtered 
frequencies  that  resulted  in  a  delay  of  the  waveform. 

Figure  6.18  (b)  shows  ten  waveforms  plotted  for  data  recorded  in  Frequency 
Response  Test  2  (first  iteration).  The  two  plotted  at  the  bottom  of  the  graph  were 
recorded  at  x  =  10  cm.  The  two  above  that  were  recorded  at  x  =  20  cm,  then  x  =  40  cm,  x 
=  80  cm,  and  finally  x  =  160  cm  is  at  the  top  of  the  graph.  At  each  position  there  is  a 
waveform  that  had  a  drive  amplitude  of  2.0  V  (solid  line)  and  one  that  had  a  drive 
amplitude  of  4.0  V  (dotted  line).  Using  the  highest  peak  as  a  reference,  there  was 
approximately  0.25  ms  between  the  two  waveforms  at  x  =  10  cm.  At  x  =  80  cm  there 
was  approximately  0.64  ms  between  the  waveforms.  This  indicated  that  the  waveform 
generated  by  the  larger  drive  amplitude  was  propagating  slower  between  these  two 
points. 
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Figure  6.18  (c)  shows  ten  waveforms  plotted  for  data  recorded  in  Frequency 
Response  Test  2  (first  iteration).  The  two  plotted  at  the  bottom  of  the  graph  were 
recorded  at  x  =  10  cm.  The  two  above  that  were  recorded  at  x  =  20  cm,  then  x  =  40  cm,  x 
=  80  cm,  and  finally  x  =  160  cm  is  at  the  top  of  the  graph.  At  each  position  there  is  a 
waveform  that  had  a  drive  amplitude  of  0.5  V  (solid  line)  and  one  that  had  a  drive 
amplitude  of  8.0  V  (dotted  line).  The  waveform  with  the  drive  amplitude  of  0.5  V  was 
scaled  up  by  a  factor  of  4  and  the  waveform  with  the  drive  amplitude  of  8.0  V  was  scaled 
down  by  a  factor  of  2.  This  was  done  so  that  the  two  waveforms  could  be  plotted  on  the 
same  graph  for  comparison.  Using  the  highest  peak  as  a  reference,  there  was 
approximately  1 .0  ms  between  the  two  waveforms  at  x  =  10  cm.  At  x  =  80  cm  there  was 
approximately  2.4  ms  between  the  waveforms.  This  verified  that  the  waveform  with  the 
greater  drive  amplitude  was  propagating  slower  between  these  two  points. 

Figure  6.19  shows  data  recorded  by  the  accelerometer  mounted  in  the  center  of 
the  shaker  foot  for  Amplitude  Response  Test  5  (first  iteration),  (a)  through  (e)  are  the 
five  frequencies  (99  Hz,  198  Hz,  396  Hz,  792  Hz,  and  1584  Hz  respectively)  measured 
while  the  radar  was  at  x  =  10  cm.  The  graphs  show  four  harmonics  normalized  by  the 
fundamental.  This  figure  indicates  that  there  was  a  significant  contribution  of  harmonics 
generated  at  the  source  in  addition  to  that  generated  by  the  propagation  path. 
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Figure  6.19  -  Amplitude  Response  Test  5 
(Gain  Setting  1),  First  Iteration;  4 
Harmonics  Normalized  by  the 
Fundamental  as  Recorded  by 
Accelerometer  Mounted  on  Center  of 
ShakerFootwhiIeRadarisatx=  10  cm  for 
5  Frequencies  (a)  99  Hz  (b)  198  Hz  (c) 
396  Hz  (d)  792  Hz  (e)  1 584  Hz 
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CHAPTER  VII 


CONCLUSIONS 

The  shaker  foot  and  propagation  path  contributed  to  the  nonlinearities  of  the 
investigated  system.  The  source  created  nonlinearities  as  indicated  by  the  harmonics 
recorded  with  accelerometers  mounted  on  the  shaker  foot.  The  propagation  path  created 
nonlinearities  due  to  a  complex,  three-dimensional  crystalline  matrix  with  pockets  of 
varying  amounts  of  water  and  air.  This  was  seen  by  the  number  of  generated  harmonics 
increasing  from  one  point  to  another  in  some  of  the  data,  despite  attenuation  being  the 
dominant  effect.  The  propagation  path  contributed  to  nonlinearities  because  of  the 
changing  solid  particle  wave  paths  and  fluctuating  viscous  and  cohesive  properties. 

The  shaker  foot  to  sand  coupling  was  an  important  contributor  to  the  results 
recorded.  The  results  showed  that  once  the  shaker  foot  to  sand  contact  was  changed  the 
results  were  not  repeatable  with  the  same  degree  of  precision.  When  comparing  sets  of 
data  taken  before  and  after  moving  the  shaker,  it  was  seen  that  the  different  fundamentals 
and  harmonics  for  the  frequency  responses  behaved  similarly  with  respect  to  each  other, 
but  changed  slightly  every  time  the  shaker  was  moved.  This  difference  was  more 
pronounced  in  the  time-domain.  The  surface  displacements  measured  by  the  radar  for 
higher  frequencies,  which  were  already  variable,  showed  the  biggest  changes  after 
changing  the  shaker  foot  to  sand  coupling. 
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When  looking  at  the  frequency  response  of  the  shaker  foot,  as  measured  by  the 
accelerometers,  the  frequency  band  that  was  most  susceptible  to  change  due  to  the  shaker 
foot  to  sand  contact  was  30  -  400  Hz.  This  was  seen  in  Figure  6.8  (b)  where  two  sets  of 
data,  with  different  shaker  foot  to  sand  contact,  were  plotted  on  the  same  graph.  The 
frequency  response  of  the  30  -  400  Hz  band  was  significantly  different.  The  changing 
frequency  response  in  this  band  was  consistent  throughout  the  data.  For  example,  Figure 
5.8  (b)  showed  a  null  around  100  Hz  that  did  not  appear  anywhere  else  in  the  data. 

The  degree  to  which  the  sources  of  nonlinearities  affect  the  propagation  of 
compressional,  surface,  and  shear  waves  was  dependent  upon  the  type  of  shaker  foot,  the 
frequencies,  and  the  amplitudes  utilized.  By  carefully  selecting  these  three  things  a  wide 
variety  of  results  were  produced.  These  results  ranged  from  near-linear  responses  to 
highly  non-linear  responses. 

Different  types  of  shakers  and  shaker  feet  affect  the  results.  Of  the  combinations 
investigated  a  rectangular  foot  with  a  length  to  width  ratio  of  approximately  10:1,  and 
enough  surface  area  to  support  the  shaker  without  burying,  produced  the  most  linear 
results.  The  degree  of  nonlinearity  was  measured  by  the  amount  of  harmonic  generation 
recorded  by  the  accelerometers  mounted  on  the  foot.  Although  the  foot  had  a  square 
cross  section,  the  length  to  thickness  ratio  was  still  such  that  a  bending-about-the-center 
mode  of  vibration  was  excited. 

The  circular  foot  produced  amplitude  responses  which  approached  saturation  and 
then  began  rising  again.  Because  the  amplitude  at  which  this  second  rise  began  changed 
depending  on  where  the  measurement  was  taken,  the  result  was  affected  by  the 
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propagation  path.  This  was  not  seen  on  the  other  two  shaker  feet.  The  likely  cause  of 
this  was  that  the  pressure  wave  and  surface  wave  had  different  saturation  thresholds.  As 
the  dominating  wave  saturated,  the  amplitude  response  had  a  decreased  slope  until  the 
other  wave,  which  was  still  increasing  with  drive  amplitude,  began  to  dominate.  This 
would  also  account  for  some  of  the  dips  in  the  amplitude  responses  of  the  higher 
frequencies. 

The  difference  between  the  round  foot  and  the  two  rectangular  feet  was  the 
surface  area  to  frontal  length  ratio.  If  the  surface  wave  was  dependent  on  the  frontal 
length  while  the  compressional  wave  was  dependent  on  the  surface  area,  the  saturation 
curves  could  behave  as  measured  due  to  the  dominating  wave  changing  from 
compressional  to  surface  during  propagation.  The  other  shaker  feet  would  have  a  surface 
wave  that  dominated  the  curve  from  the  beginning  and  therefore  did  not  produce  this 
two-rise  effect.  The  foot  motion  must  be  well  documented  in  order  to  accurately 
represent  the  source  in  the  computer  model.  It  is  also  important  to  ensure  that  the  power 
amplifier  and  shaker  are  properly  matched  so  that  an  impedance  mismatch  does  not 
increase  the  nonlinearity  of  the  source. 

The  range  of  frequencies  used  to  generate  the  wave  also  impacted  the  results. 
Frequencies  less  than  600  Hz  propagated  well  through  the  sand.  Frequencies  higher  than 
600  Hz  were  highly  vulnerable  to  attenuation,  particularly  once  the  surface  of  the  sand 
dried.  A  flatter  frequency  response,  with  less  variability  of  the  higher  frequencies,  was 
achieved  with  lower  amplitudes  as  seen  in  Figure  6.8  (a).  Doubling  the  amplitude 
doubled  the  surface  displacement  for  the  100  -  1200  Hz  frequency  band.  The  1300  - 
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2000  Hz  frequency  band  also  showed  this  behavior  up  to  x  =  20  cm  where  the  0.25  V 
amplitude  began  to  exhibit  the  variability  of  the  higher  frequencies. 

The  30  Hz  -  2000  Hz  frequency  range  used  for  the  acousto-electromagnetic  mine 
detection  technique  appears  to  be  very  well  suited  for  this  task.  The  generation  of 
harmonics  by  the  lower  frequencies  helps  to  increase  the  surface  displacements  of  the 
higher  frequency  range.  Anything  above  2000  Hz  however,  would  attenuate  so  quickly 
that  no  matter  how  much  contribution  from  lower  frequency  harmonics  was  present  there 
would  not  be  enough  energy  in  these  frequencies  to  propagate  an  appreciable  distance. 
The  higher  frequencies  are  also  less  useftil  in  that  the  variability  of  their  amplitude 
responses  would  produce  nulls  at  unpredictable  locations. 

Increasing  the  drive  amplitude  caused  system  nonlinearities  as  expected.  The 
threshold  of  linearity  changed  as  a  function  of  distance  and  frequency.  Small  amplitudes 
propagated  well  enough  to  be  measured  by  the  radar  at  the  furthest  point  tested.  These 
smaller  amplitudes  had  a  much  flatter  frequency  response,  although  there  was  a  null 
around  1250  Hz  for  the  two  smallest  amplitudes  tested.  This  null  was  due  to  the 
frequency  response  of  the  foot  and  what  appeared  to  be  some  destructive  interference. 

As  the  amplitude  increased  however,  the  surface  displacements  due  to  the  lower 
frequencies  rose  faster  than  surface  displacements  due  to  the  higher  frequencies.  This 
was  important  to  note  when  trying  to  make  the  computer  model  match  the  actual 
experiments.  Increasing  the  drive  amplitude  increased  nonlinearity  by  first  driving  the 
shaker  foot  such  that  harmonics  were  generated  and  also  caused  the  wave  to  propagate  in 
a  nonlinear  way  through  the  unconsolidated  soil  matrix. 
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The  threshold  of  surface  displacement  for  these  experiments  was  approximately 
3000  nm  as  measured  at  x  =  10  cm.  These  displacements  occurred  in  the  200  -  400  Hz 
frequency  band,  despite  the  accelerations  measured  on  the  shaker  foot  being  greatest 
around  1500  Hz.  Regardless  of  the  amplitude  or  frequency  used,  the  radar  never 
measured  any  surface  displacements  greater  than  3000  nm.  The  threshold  did  not  change 
for  different  shakers  and  shaker  feet  combinations.  The  20  pound  shaker  with  small 
rectangular  foot  also  saturated  the  sand  at  this  point  even  though  the  surface  area  of  the 
smaller  shaker  foot  was  almost  3.5  times  less  than  that  of  the  shaker  foot  used  in 
Experiment  Two.  The  larger  shaker  and  shaker  foot  used  more  current  without  burying 
into  the  sand,  but  the  additional  current  was  used  to  drive  the  heavier  foot  and  did  not 
increase  the  magnitude  of  displacement  in  the  sand. 
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CHAPTER  Vin 


RECOMMENDATIONS 

There  are  many  different  shaker  feet  that  could  be  used  in  the  acousto- 
electromagnetic  mine  detection  technique.  A  study  of  these  possibilities  should  be 
conducted.  The  three  possibilities  examined  as  part  of  this  research  produced  very 
different  results  that  indicated  that  the  many  other  possibilities  could  turn  up  a 
configuration  much  more  suited  for  mine  detection.  Shaker  feet  could  also  be  made  that 
did  not  have  resonances  and  had  only  one  mode  of  vibration  excited  for  the  frequencies 
and  drive  amplitudes  used.  A  shaker  foot  very  similar  to  the  one  used  in  Experiment 
Two  could  be  made  with  the  same  surface  area  but  thicker  cross  section  in  order  to 
achieve  this. 

If  an  investigation  of  shaker  feet  was  conducted,  it  should  focus  on  those  with  a 
large  (>  10: 1)  length  to  width  ratio.  The  amplitude  response  of  the  20  pound  shaker  with 
circular  foot  showed  that  the  small  length  to  width  ratio  (1:1  in  this  case)  resulted  in  a 
greater  degree  of  nonlinearity.  Making  a  large  round  foot  for  the  100  pound  shaker 
would  create  the  same  nonlinearites. 

The  way  that  the  shaker  foot  couples  with  the  sand  could  also  be  changed.  The 
only  technique  examined  thus  far  was  placing  the  foot  on  top  of  the  sand  surface  and 
relying  on  the  weight  of  the  shaker  to  keep  shaker  foot  to  sand  contact  steady.  Different 
foot  -  sand  couplings  should  be  investigated  to  determine  if  another  technique  is  more 
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suited  for  the  production  of  surface  waves.  One  example  would  be  something  on  the  foot 
that  penetrates  the  sand,  such  as  nails,  which  is  one  configuration  used  by  researchers  at 
the  University  of  Texas  [4], 

This  research  looked  at  what  was  happening  only  along  the  x-axis  of  the  sandbox. 
At  this  point  it  would  be  beneficial  to  expand  the  research  to  looking  in  two  dimensions. 
Very  little  is  known  about  the  directivity  of  the  various  shaker  -  shaker  feet 
combinations.  This  directivity  changes  depending  on  which  foot  is  used  because  of  the 
unique  shaker  foot  motions,  sizes,  and  shapes.  Since  mine  detection  dwells  in  a  limited 
three-dimensional  space,  which  is  very  large  on  the  surface,  determining  the  directivity  of 
these  sources  will  become  important. 

The  near-field  radiation  pattern  of  the  shaker  is  very  complicated.  The  reason  for 
increasing  surface  displacements  for  the  fundamental,  between  10  cm  and  40  cm,  in 
Frequency  Response  Test  4,  is  still  unknown.  The  low  frequencies  used  and  the  size  of 
the  shaker  foot  resulted  in  a  near-field  of  appreciable  size.  It  could  be  worth  the  effort  to 
try  and  characterize  this  near-field.  If  an  array  of  sources  is  ever  planned  for 
implementation,  determining  what  the  behavior  of  frequencies  in  the  near-field  is  will  be 
even  more  important. 

Some  mine  detection  tests  should  be  conducted  with  much  lower  drive 
amplitudes.  Frequency  Response  Test  4  showed  that  lower  drive  levels  produced  a  flatter 
frequency  response.  There  was  not  as  much  variability  in  the  form  of  frequency  response 
nulls.  The  surface  displacements  due  to  higher  frequencies  were  also  greater  Avith  respect 
to  the  displacements  of  the  lower  frequencies  than  they  were  in  the  other  frequency 
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response  tests.  It  might  be  beneficial,  for  interrogation  of  very  small  objects,  to  lower  the 
drive  amplitude  so  that  higher  frequencies  do  not  become  variable  as  they  do  for  higher 
amplitudes.  If  this  were  done,  the  incident  signal  duration  would  probably  have  to  be 
lengthened  to  improve  the  signal  to  noise  ratio.  This  would  result  in  a  trade-off  of  time  to 
conduct  a  scan  for  better  high  frequency  propagation. 

Finally,  two-dimensional  scans  for  mines  using  an  incident  signal  of  reduced 
bandwidth  should  be  tested  in  the  event  that  time  is  more  important  than  interrogation 
with  frequencies  greater  than  1200  Hz.  Because  of  the  rapid  attenuation  of  higher 
frequencies,  it  may  not  be  worth  using  a  chirp  that  contains  frequencies  between  1200  Hz 
and  2000  Hz.  More  than  a  second  could  be  saved  for  each  measurement  by  utilizing  a  30 
-  1200  Hz  chirp.  Under  the  current  procedure  for  conducting  two-dimensional  scans,  this 
would  reduce  the  9.5  hour  scan  by  about  10  percent.  This  will  become  more  and  more 
important  in  the  future  as  the  research  heads  towards  practical  implementation.  If  neither 
time  nor  frequencies  above  1200  Hz  were  critical,  then  the  same  length  signal  using  a  30 
Hz  to  1200  Hz  chirp  could  be  used  to  improve  the  signal  to  noise  ratio. 
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APPENDIX  A 


EXPERIMENT  ONE  DETAILS 

Design  of  Experiments 

In  order  to  determine  what  type  of  experiments  would  be  the  most  effective  and 
efficient  for  this  research,  some  initial  tests  were  run  on  the  sandbox.  There  were  three 
specific  factors  that  needed  to  be  found.  First,  the  maximum  amplitude  that  the  shaker 
could  be  driven  without  burying  itself  in  the  sand,  and  the  minimum  amplitude  appearing 
above  the  noise  floor  needed  to  be  found  out.  This  would  set  the  upper  and  lower  limits 
of  input  voltages  for  the  experiments.  Second,  the  maximum  duration  a  scan  could  be 
run,  without  experiencing  nonlinearities  due  to  the  sand  drying,  was  needed.  This  would 
determine  how  often  the  scans  needed  to  be  stopped  in  order  to  rewet  and  recompact  the 
sand.  Third,  the  minimum  duration  of  the  input  signal,  while  still  recording  accurate 
data,  needed  to  be  found  in  order  to  minimize  scan  time. 

Shaker  Amplitude  Range 

Total  amplitude  in  this  experimental  setup  was  produced  by  a  combination  of  the 
DAC  and  an  amplifier.  In  order  to  find  the  low  end  of  the  amplitude  range  the  gain  on 
the  amplifier  was  turned  all  the  way  up.  Lower  and  lower  values  of  amplitude  were  then 
entered  into  the  computer  for  the  board  until  a  value  of  0.03  volts  was  found  to  be  the 
smallest  value  that  would  still  register  above  the  noise  floor.  Harmonics  were  not  seen, 
however,  until  the  value  entered  in  the  computer  was  approximately  0.15  volts.  The 
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lower  end  value  of  the  amplitude  range  was  then  chosen  as  0.06  volts  entered  in  the 
computer  and  the  amplifier  set  at  maximum  gain.  This  ensured  that  the  data  collected 
began  in  the  linear  region  (no  harmonics  produced). 

The  goal  for  the  experiment  was  to  use  approximately  24  different  amplitudes  at 
each  point  tested.  About  10%  of  the  amplitudes  on  the  upper  end  of  the  range  would 
cause  the  shaker  foot  to  settle  into  the  ground  during  a  scan.  A  series  of  trials  was  run  in 
order  to  determine  where  the  amplitude  would  have  to  be  set  for  the  shaker  foot  to  settle 
into  the  sand.  It  was  determined  that  if  the  voltage  entered  into  the  computer  was  greater 
than  0.86  V,  with  the  gain  on  the  amplifier  all  the  way  up,  this  occurred.  Therefore,  0.98 
V  was  chosen  as  the  upper  end  of  the  entered  value  of  voltage. 

Maximum  Scan  Duration 

Drying  Test  Number  1  An  experiment  was  conducted  in  order  to  determine  how 
long  the  sand’s  propagation  properties  remained  constant  before  drying  effects  became 
noticeable  in  the  data.  A  program  was  written  to  conduct  a  41  point  scan  (0-120  cm  at 
3  cm  increment)  along  the  x-axis  every  hour.  A  3.5  second  chirp  from  30  Hz  to  2000  Hz 
was  used  as  the  input  signal.  The  sand  was  prepared  for  scanning  and  the  program  was 
executed.  This  data  was  recorded  for  a  72  hour  period  (73  scans). 

The  velocity  of  the  surface  wave  remained  at  approximately  91  m/s  during  the 
entire  72  hour  period.  This  velocity  was  calculated  by  measuring  points  on  the  waterfall 
graphs  so  a  great  deal  of  precision  could  not  be  achieved.  The  velocity  most  likely 
decreased  at  a  rate  that  was  too  small  to  detect  as  the  sand  dried,  however  the  surface 
wave  velocity  did  remain  somewhere  in  the  90.5  m/s  -  92  m/s  range. 
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Frequency  propagation  was  also  studied.  During  the  first  ten  hours  of  drying 
there  was  no  significant  loss  of  frequency  propagation.  After  ten  hours  however, 
frequencies  greater  than  900  Hz  appeared  to  show  a  decrease  in  ability  to  propagate  over 
the  full  120  cm  of  the  scan  region.  At  the  time,  this  was  taken  to  mean  that  the  sandbox 
would  need  to  be  rewetted  and  recompacted  every  eight  to  ten  hours  to  ensure 
propagation  of  the  higher  frequencies. 

The  interpretation  of  the  drying  tests  was  not  precisely  correct  due  to  the  fact  that 
a  chirp  was  being  used  as  the  incident  signal.  This  prevented  a  lot  of  energy  being  placed 
into  any  one  frequency  band  and  the  signal  to  noise  ratio  was  not  as  good  as  it  should 
have  been.  Also,  by  the  time  Experiment  Two  was  conducted,  it  was  realized  that  the 
higher  frequency’s  propagation  ability  actually  drops  off  within  the  first  couple  of  hours 
due  to  drying  and  the  lower  frequencies  remain  able  to  propagate  regardless  of  moisture 
conditions. 

Drying  Test  Number  2  It  was  determined  through  several  sample  data  collections 
that  drying  effects  might  be  affecting  some  of  the  frequencies  when  a  sinusoidal  input  is 
used  instead  of  the  30  -  2000  Hz  chirp.  A  second  drying  test  was  conducted  to  examine 
this.  The  radar  was  positioned  at  point  (40,0).  Every  15  minutes  a  scan  was  taken  from 
100  -  2000  Hz,  at  100  Hz  increments,  and  amplitude  equal  to  0.5  V.  Each  frequency 
input  was  a  3.5  second  sinusoid  with  a  0.5  second  settling  time.  Only  two  seconds  of  the 
3.5  available  was  used  for  data  processing  in  order  to  minimize  the  impact  of  any  start-up 
or  shut-down  transient  signal. 
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Once  this  data  had  been  collected,  the  effects  of  the  sand  drying  were  examined 
by  looking  at  the  frequency  response  of  the  surface  displacements  for  the  fundamental 
and  harmonics.  The  changes  were  compared  over  time.  The  results  were  also  examined 
for  200,  500,  800,  1 100,  1400,  and  1700  Hz  by  comparing  surface  displacements  for  the 
fundamentals  as  a  function  of  time. 

The  results  of  the  second  drying  test  showed  that  certain  frequencies  were  more 
susceptible  to  effects  of  drying  than  others.  For  example,  the  low  range  of  frequencies 
(100  -  500  Hz)  experience  very  little  change  from  drying  effects.  In  the  mid-range  of 
frequencies  (600  -  1500  Hz),  drying  caused  most  of  the  frequency  displacements  to 
diminish.  In  the  high  range  of  frequencies  (1600  -  2000  Hz),  a  variety  of  things  took 
place.  1600  Hz  remained  about  the  same,  1700  and  1800  Hz  increased  as  the  sand  dried, 
and  1900  and  2000  Hz  decreased  as  the  sand  dried.  The  harmonics  behaved  similar  to 
the  fundamental  frequencies  but  were  less  predictable. 

The  displacement  of  the  frequencies  remained  relatively  constant  over  a  ten  hour 
period.  The  largest  changes  occurred  within  the  first  hour  after  preparing  the  sand.  With 
the  exception  of  500  Hz,  the  harmonics  supported  this  observation.  These  results 
indicated  that  if  the  sand  was  reconditioned  about  every  eight  hours,  the  effects  of  drying 
would  be  minimized.  Also,  the  sand  would  be  allowed  to  reach  a  quasi-equilibrium  by 
waiting  one  hour  from  the  time  of  reconditioning  before  data  collection  would 
commence. 
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Minimum  Input  Signal  Duration 


A  program  was  written  to  test  20  different  frequencies  (100-2000  Hz  in  100  Hz 
increments),  with  amplitude  equal  to  0.5  V,  at  two  different  points  (x  =  40  cm  and  x  = 

120  cm).  These  tests  were  run  three  different  times  to  determine  if  the  signal  duration 
had  an  effect  on  the  output.  The  three  time  windows  used  were  4.096  seconds,  2.048 
seconds,  and  1.024  seconds.  The  time  windows  were  composed  of  the  signal  followed  by 
approximately  0.5  seconds  of  settling  time.  The  settling  time  changed  slightly  depending 
on  the  frequency  being  tested  so  that  the  input  sine  wave  ended  after  an  integer  number  of 
periods  each  time. 

The  fundamental  frequencies  and  some  of  the  harmonics  were  seen  using  each 
one  of  these  signal  durations.  However,  the  longer  input  signal  yielded  more  harmonics 
registering  outside  of  the  noise  level.  It  was  suspected  that  the  one  second  signal  would 
yield  the  same  results  as  the  two  and  four  second  signal  with  slight  differences  in  the 
signal  to  noise  ratio.  Because  different  signal  lengths  were  producing  a  different  number 
of  harmonics,  the  data  indicated  that  there  might  be  a  start-up  and/or  a  shut-down 
transient  present  which  was  having  less  of  an  impact  as  it  was  averaged  out  over  the 
longer  signal  duration. 

The  experiment  was  run  again  but  this  time  the  time-domain  data  was  saved.  This 
allowed  the  fast  Fourier  transform  to  be  taken  over  different  time  windows.  This  was 
done  two  different  ways.  First,  a  one  second  time  window,  shifted  a  half  second  at  a 
time,  was  used.  Then  a  half  second  time  window,  shifted  a  half  second  at  a  time,  was 
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used.  In  both  cases  the  results  were  different  for  the  first  time  window  which  verified  the 


presence  of  a  start-up  transient. 

As  a  result  of  this,  the  decision  was  made  to  use  a  4.096  second  time  window. 
Approximately  3.5  seconds  would  be  the  signal  and  the  remaining  would  be  settling  time. 
Only  2.048  seconds  of  the  data  collected  would  be  used  however,  so  that  any  start-up  and 
shut-down  transients  could  be  eliminated  before  processing. 

Data  Collection 

The  data  was  collected  in  two  different  phases.  The  first  phase  was  to  get  the  data 
for  the  frequency  and  amplitude  responses.  The  second  phase  gathered  data  in  order  to 
separate  the  pressure  wave  from  the  surface  wave  so  that  individual  contributions  could 
be  studied.  The  third  phase  was  gathering  information  for  altered  relative  energy 
contents  in  the  pressure  and  surface  waves.  Experiment  One  ended  and  planning  for 
Experiment  Two  began  before  the  third  phase  was  completed. 

Phase  I 

For  this  phase,  a  “scan”  consisted  of  measuring  180  different  frequencies  at  a 
certain  amplitude  and  point.  The  frequencies  ranged  from  33  Hz  to  2002  Hz  by  steps  of 
1 1  Hz.  A  scan  was  taken  for  24  different  amplitudes,  at  each  of  three  different  positions, 
for  a  total  of  72  scans.  The  24  amplitudes  ranged  from  0.06  volts  to  0.96  volts  on  the 
board,  in  steps  of  0.04  volts,  with  the  gain  on  the  amplifier  all  the  way  up. 

The  first  24  scans  were  taken  with  the  radar  at  x  =  40  cm,  y  =  0  cm,  and  z  =  0  cm 
on  the  positioner.  The  distance  from  the  lead  edge  of  the  shaker  foot  to  the  center  of  the 
waveguide  was  actually  71.2  cm.  The  actual  distance  from  the  surface  of  the  sand  to  the 
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bottom  edge  of  the  waveguide  was  2.0  cm.  The  power  meter  on  the  radar  read  -33  dBm 
(+/-  0.8  dBm)  when  it  was  raised  to  z  =  30  cm  on  the  positioner  and  microwave  scattering 
foam  was  placed  under  the  waveguide.  The  power  meter  read  -1.55  dBm  at  the  actual 
position  where  the  data  was  taken. 

The  sand  was  watered  down,  compacted,  and  allowed  to  dry  for  one  hour  prior  to 
beginning  the  data  collection.  The  scans  were  done  in  order  of  increasing  amplitude  so  as 
not  to  disturb  the  sand  under  the  shaker  foot.  15  scans  were  completed  in  8.25  hours  at 
which  point  a  pause  was  taken  to  rewater  and  recompact  the  sand.  Prior  to  starting  again, 
the  radar  power  was  checked.  It  was  —52  dBm  (+/-  4  dBm)  at  z  =  30  cm  over  the  foam 
and  -0.91  dBm  at  the  measuring  position.  Five  more  scans  were  completed  before  the 
rewatering  and  recompacting  procedure  was  once  again  performed.  At  this  point  the 
radar  power  meter  read  —32  dBm  (+/- 1.5  dBm)  at  z  =  30  cm  over  the  foam  and  —0.94 
dBm  (+/- 1.5  dBm)  at  the  measuring  point.  From  this  point  on,  the  data  collection  had  to 
be  stopped  after  every  scan  in  order  to  compact  under  the  shaker  foot  because  the 
amplitude  was  such  that  the  shaker  foot  was  burying  itself  in  the  sand.  The  power  meter 
on  the  radar  read  -2.53  dBm,  -2.59  dBm,  and  -2.74  dBm  at  the  measuring  position  prior 
to  the  last  three  scans. 

Throughout  the  entire  process  described  above,  calibration  scans  were  taken 
before  and  after  each  data  collection  scan.  These  consisted  of  measuring  20  frequencies 
(100  Hz  to  2000  Hz  by  100  Hz  increments)  at  the  same  amplitude  (0.5  on  the  computer 
with  the  amplifier  gain  all  the  way  up)  and  same  position  (x  =  40  cm,  y  =  0  cm,  and  z  —  0 
cm  on  the  positioner)  each  time.  By  measuring  the  exact  same  thing  before  and  after 
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each  scan,  a  comparison  between  the  two  could  be  made  to  see  how  much  effect  drying 
had  during  the  scan.  These  calibration  scans  could  be  used  to  compare  the  condition  of 
the  sand  during  any  scan  regardless  of  when  it  was  taken. 

The  second  24  scans  were  taken  with  the  radar  at  x  =  80  cm,  y  =  0  cm,  and  z  =  0 
cm  on  the  positioner.  The  distance  from  the  lead  edge  of  the  shaker  foot  to  the  center  of 
the  waveguide  was  actually  111.7  cm.  The  actual  distance  from  the  surface  of  the  sand  to 
the  bottom  edge  of  the  waveguide  was  2. 1  cm.  The  power  meter  on  the  radar  read  -32.6 
dBm  (+/-  0.4  dBm)  when  it  was  raised  to  z  =  30  cm  on  the  positioner  and  microwave 
scattering  foam  was  placed  under  the  waveguide.  The  power  meter  read  -2. 19  dBm  (+/- 
0.01  dBm)  at  the  actual  position  where  the  data  was  taken. 

The  sand  was  watered  down,  compacted,  and  allowed  to  dry  for  one  hour  prior  to 
beginning  the  data  collection  just  as  it  had  been  done  for  the  scans  at  x  =  40  cm.  The 
scans  were  again  done  in  order  of  increasing  amplitude  so  as  not  to  disturb  the  sand  under 
the  shaker  foot.  15  scans  were  completed  in  8.5  hours  at  which  point  the  sand  was 
rewatered  and  recompacted.  Prior  to  starting  again,  the  radar  power  was  checked.  It  was 
-36.6  dBm  (+/-  0.5  dBm)  at  z  =  30  cm  over  the  foam  and  -3.00  dBm  (+/-  0.01  dBm)  at 
the  measuring  position.  Three  more  scans  were  completed  before  the  rewatering  and 
recompacting  procedure  was  once  again  performed.  At  this  point  the  radar  power  meter 
read  -3.52  dBm  (+/-  0.01  dBm)  at  the  measuring  point.  One  scan  was  completed  and  the 
reconditioning  procedure  was  repeated  with  the  radar  power  meter  reading  -3.70  dBm 
(+/-  0.01  dBm)  at  the  measuring  point.  Two  more  scans  were  completed  and  then  the 
data  collection  had  to  be  stopped  after  every  scan  in  order  to  compact  under  the  shaker 
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foot  because  the  amplitude  was  such  that  the  shaker  foot  was  burying  itself  in  the  sand. 
The  power  meter  on  the  radar  read  -3.77  dBm,  -3.80  dBm,  and  -3.86  dBm  (+/-  0.01  dBm 
for  each)  at  the  measuring  position  prior  to  the  last  three  scans. 

Throughout  the  entire  process  described  above,  calibration  scans  were  again  taken 
before  and  after  each  data  collection  scan.  The  procedure  for  these  calibration  scans  was 
exactly  like  the  procedure  described  above  for  the  point  at  x  =  40  cm. 

The  final  24  scans  were  taken  with  the  radar  at  x  =  120  cm,  y  =  0  cm,  and  z  =  0 
cm  on  the  positioner.  The  distance  from  the  lead  edge  of  the  shaker  foot  to  the  center  of 
the  waveguide  was  actually  151.5  cm.  The  actual  distance  from  the  surface  of  the  sand  to 
the  bottom  edge  of  the  waveguide  was  1.5  cm.  The  power  meter  on  the  radar  read  -30.5 
dBm  (+/-  0.5  dBm)  when  it  was  raised  to  z  =  30  cm  on  the  positioner  and  microwave 
scattering  foam  was  placed  under  the  waveguide.  The  power  meter  read  +0.69  dBm  (+/- 
0.01  dBm)  at  the  actual  position  where  the  data  was  taken. 

The  sand  was  watered  down,  compacted,  and  allowed  to  dry  for  one  hour  and  20 
minutes  prior  to  beginning  the  data  collection  to  once  again  allow  it  to  reach  a  state  of 
quasi-equilibrium.  Beginning  again  with  the  lowest  amplitude,  8  scans  were  completed 
in  4.5  hours  at  which  point  the  sand  was  rewatered  and  recompacted.  Prior  to  starting 
again,  the  radar  power  was  checked.  It  was  +0.85  dBm  (+/-  0.01  dBm)  at  the  measuring 
position.  Ten  more  scans  were  completed  in  5.75  hours  before  the  sand  had  to  be 
rewatered  and  recompacted.  At  this  point  the  radar  power  meter  read  +2.90  dBm  (+/- 
0.01  dBm)  at  the  measuring  point.  Three  more  scans  were  then  completed  before  the 
data  collection  had  to  be  stopped  after  every  scan  in  order  to  compact  under  the  shaker 
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foot  to  prevent  it  from  burying  itself  in  the  sand.  The  power  meter  on  the  radar  read 
+1.41  dBm,  0.93  dBm,  and  0.79  dBm  (+/-  0.01  dBm  for  each)  at  the  measuring  position 
prior  to  the  last  three  scans. 

Once  again,  calibration  scans  were  taken  before  and  after  each  data  collection 
scan  during  the  entire  process  described  above.  The  procedure  for  these  calibration  scans 
was  exactly  like  the  procedure  described  above  for  the  point  at  x  =  40  cm. 

Phase  II 

The  data  for  this  phase  was  collected  in  the  same  fashion  as  the  data  for  the  points 
at  40,  80,  and  120  cm  in  Phase  I.  The  only  difference  was  the  point  at  which  the  data  was 
collected.  During  this  part  of  the  experiment,  the  point  at  x  =  190  cm  was  used  in  order 
to  allow  the  pressure  wave  and  surface  wave  to  separate  in  time.  This  point  is  at  the  far 
limit  of  the  positioner  in  the  experimental  setup.  By  allowing  the  pressure  wave  and 
surface  wave  to  separate,  the  contributions  of  each  to  the  overall  displacement  were  to  be 
measured  separately.  The  frequencies  and  amplitudes  used  were  the  same  as  those  in 
Phase  I  of  the  data  collection. 

The  distance  from  the  lead  edge  of  the  shaker  foot  to  the  center  of  the  waveguide 
was  actually  220.7  cm.  The  actual  distance  from  the  surface  of  the  sand  to  the  bottom 
edge  of  the  waveguide  was  1 .9  cm.  The  power  meter  on  the  radar  read  -29.9  dBm  (+/- 
0.4  dBm)  when  it  was  raised  to  z  =  30  cm  on  the  positioner  and  microwave  scattering 
foam  was  placed  under  the  waveguide.  The  power  meter  read  +0. 16  dBm  at  the  actual 
position  where  the  data  was  taken. 
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The  sand  was  watered  down,  compacted,  and  allowed  to  dry  for  one  hour  and  five 
minutes  prior  to  beginning  the  data  collection.  10  scans  (amplitudes  0.30  -  0.66)  were 
completed  before  pausing  to  rewater  and  recompact  the  sand.  Prior  to  starting  again,  the 
radar  power  was  checked.  It  was  -28. 15  dBm  (+/-  0. 1  dBm)  at  z  =  30  cm  over  the  foam 
and  +2.29  dBm  (+/-  0.01  dBm)  at  the  measuring  position.  Five  more  scans  (amplitudes 
0.06  -  0.22)  were  completed  before  the  rewatering  and  recompacting  procedure  was  once 
again  performed.  At  this  point  the  radar  power  meter  read  -27.5  dBm  (+/-  0. 1  dBm)  at  z 
=  30  cm  over  the  foam  and  -0.14  dBm  (+/-  0.01  dBm)  at  the  measuring  point.  This  time 
four  scans  (amplitudes  0.26  and  0.70  —  0.78)  were  completed  before  reconditioning  the 
sand.  The  radar  power  meter  then  read  -25. 15  dBm  (+/-  0. 1  dBm)  at  Z  =  30  cm  over  the 
foam  and  +0.96  dBm  at  the  measuring  point.  From  this  point  on,  the  data  collection  had 
to  be  stopped  after  every  scan  (amplitudes  0.82  -  0.98)  in  order  to  compact  under  the 
shaker  foot  because  the  amplitude  was  such  that  the  shaker  foot  was  burying  itself  in  the 
sand. 

Throughout  the  entire  process  described  above,  calibration  scans  were  taken 
before  and  after  each  data  collection  scan.  These  calibration  scans  were  conducted  in  an 
identical  manner  to  those  described  in  Phase  I  of  the  data  collection. 

Results 

As  mentioned  in  Chapter  IV,  the  reoccurring  problems  with  the  data  were  the 
increasing  noise  floor  as  amplitude  increased,  the  erratic  ends  of  the  amplitude  response 
curves  as  the  shaker  buried  or  was  moved,  and  not  enough  amplitudes  measured  to 
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produce  a  smooth  curve  beginning  in  the  linear  region  near  the  noise  floor  and  increasing 
to  saturation. 
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APPENDIX  B 


ADDITIONAL  FREQUENCY  RESPONSE  GRAPHS 

This  appendix  contains  a  complete  set  of  data  for  one  of  the  frequency  response 
measurements.  The  five  figures  come  from  the  second  iteration  of  the  second  frequency 
response  test.  This  data  was  taken  on  Gain  Setting  1 .  Figure  B.  1  shows  the  fundamental 
and  four  harmonics  taken  at  x  =  10  cm.  The  five  graphs  in  this  figure  show  the  results  for 
the  following  five  amplitudes:  0.5  V,  1.0  V,  2.0  V,  4.0  V,  and  8.0  V.  Figures  B.2 
through  B.5  show  data  taken  at  x  =  20  cm,  40  cm,  80  cm,  and  160  cm  respectively.  The 
same  five  amplitudes  were  used  in  each  of  these  figures. 
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Figure  B.l  -  Frequency  Response  Test  2 
(Gain  Setting  1),  Second  Iteration: 
Fundamental  and  4  Harmonics  at  x  =  10 
cm  (a)  Amplitude = 0. 5  V  (b)  Amplitude = 
1.0  V(c)  Amplitude  =  2.0  V  (d)  Amplitude 
=  4.0  V  (e)  Amplitude = 8.0  V 
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Figure  B.3  -  Frequency  Response  Test  2 
(Gain  Setting  1),  Second  Iteration: 
Fundamental  and  4  Harmonics  at  x  =  40 
cm  (a)  Amplitude  =  0. 5  V  (b)  Amplitude = 
1.0  V(c)Amplitude=2.0  V  (d)  Amplitude 
= 4.0  V  (e)  Amplitude = 8.0  V 
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Figure  B.4  -  Frequency  Response  Test  2 
(Gain  Setting  1),  Second  Iteration; 
Fundamental  and  4  Harmonics  at  x  =  80 
cm  (a)  Amplitude  =  0. 5  V(b)  Amplitude = 
1 .0  V  (c)  Amplitude = 2. 0  V  (d)  Amplitude 
= 4.0  V  (e)  Amplitude = 8.0  V 
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Figure  B.5  -  Frequency  Response  Test  2 
(Gain  Setting  1),  Second  Iteration: 
Fundamental  and  4  Harmonics  at  x  =  160 
cm  (a)Amplitude=0.5 V(b)Amplitude= 
1.0  V(c)Amplitude=2.0V  (d)  Amplitude 
= 4.0  V  (e)  Amplitude = 8.0  V 
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ADDITIONAL  AMPLITUDE  RESPONSE  GRAPHS 

This  appendix  contains  a  complete  set  of  data  for  one  of  the  amplitude  response 
measurements.  The  five  figures  come  from  the  second  iteration  of  the  fifth  amplitude 
response  test.  This  data  was  taken  on  Gain  Setting  1.  Figure  C.  1  shows  four  harmonics 
normalized  by  the  fundamental  at  x  =  10  cm.  The  five  graphs  ((a)  —  (e))  are  the  data 
taken  for  99  Hz,  198  Hz,  396  Hz,  792  Hz,  and  1584  Hz  respectively.  Figures  C.2  through 
C.5  show  the  same  information  for  x  =  20  cm,  40  cm,  80  cm,  and  160  cm  respectively. 
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Figure  C.l  -  Amplitude  Response  Test  5 
(Gain  Setting  1),  Second  Iteration;  4 
Harmonics  Normalized  by  the 
Fundamental  at  x  =  10  cm  (a)  99  Hz  (b) 
198  Hz(c)396Hz(d)  792  Hz(e)  1584Hz 
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Figure  C.4  -  Amplitude  Response  Test  5 
(Gain  Setting  1),  Second  Iteration;  4 
Harmonics  Normalized  by  the 
Fundamental  at  x  =  80  cm  (a)  99  Hz  (b) 
1 98  Hz  (c)  3 96  Hz  (d)  792  Hz  (e)  1 5  84  Hz 
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APPENDIX  D 


MATLAB  CODE 

This  appendix  contains  some  of  the  MATLAB  code  used  to  process  the  data 
during  Experiment  Two.  Program  1  is  an  example  of  how  the  Lab  VIEW  files  were  read 
into  MATLAB  for  the  frequency  response  graphs  and  how  they  were  broken  into 
subgroups  and  saved  as  *.mat  files.  Program  2  is  an  example  of  how  the  saved  *.mat 
files  were  used  to  plot  five  different  positions  on  a  graph  of  displacement  versus 
fi'equency  for  a  given  amplitude.  Program  3  is  an  example  of  how  plots  were  generated 
of  displacement  versus  frequency  for  the  fundamental  and  five  harmonics  at  a  given 
position  and  amplitude. 

In  addition  to  the  program  examples  contained  in  this  appendix,  programs  were 
written  to  plot  five  different  amplitudes  on  a  graph  of  displacement  versus  frequency  for 
a  given  position,  plot  the  fundamental  (with  or  without  a  comparison  to  a  second 
experiment)  at  a  given  position  and  amplitude,  plot  the  fundamental  and  harmonics 
normalized  by  the  drive  signal  at  a  given  position  and  amplitude,  and  plot  the  harmonics 
normalized  by  the  fimdamental  at  a  given  position  and  amplitude.  The  same  types  of 
programs  were  written  to  process  the  accelerometer  data.  All  of  these  above  mentioned 
programs  were  written  for  four  different  frequency  response  data  sets.  Similar  programs 
were  also  written  for  the  amplitude  response  data. 
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Program  1 


%  This  program  takes  a  transfer  fiinction  measurement  and  breaks  all  of  the  data  into  five 
%  matrices  (amplitude  groups)  for  the  six  .bin  files  saved. 

%  FIRST  TRANSFUN  MEASUREMENT 


clear  all 

%  Open  and  read  files 

pathname=strcatCc:\blace\datafiIes\F20000 126-1453 14V); 

fileF=strcat(pathname,’parameters.bin’); 

fileG=strcat(palhname,’fund_harm.bin'); 

fileH=strcat(palhname,'noise.bin'); 

fileJ=strcal(pathname,'accel_A.bin'); 

fiIeK=strcat(pathname,'accel_B.bin’); 

fileL=strcat(pathname,'accel_noise.bin^; 

fidF=fopen(fileF,V,'ieee-be'); 

fidG=fopen(fileG,'r’,‘ieee-be'); 

fidH=fopen(fileH,’r’,'ieee-be'); 

fidJ=fopen(fileJ,'r','ieee-be'); 

fidK=fopen(fileK,'f',’ieee-be’); 

fidL=fopen(fileL,'r','ieee-be'); 

F=fi-ead(fidF,‘fioat32’); 

G=fi'ead(fidG,'float32'); 

H=fi'ead(fidH,*float32’); 

J=fi‘ead(fidJ,’float32'); 

K-fi'ead(fidK,‘float32*); 

L=fi^ad(fidL,'float32'); 

%  Initialize  matrices 
parameters0_5 ( 1 :5400, 1  )=0; 
parameters  1_0( 1 : 5400, 1 )=0; 
parameters2_0( 1 :5400, 1 )=0; 
parameters4_0(l  :5400, 1 )=0; 
parameters8_0( 1 :5400, 1 )=0; 
fund_hann0_5(l  :2 1 600, 1 )=0; 
fund_harml_0(l:21600,l)=0; 
fund_harm2_0(l  :2 1600, 1)=0; 
fund_harm4_0(l:21600,l)=0; 
fund_harm8_0(l  :2 1600, 1)=0; 
noise0_5(l:10800,l)=0; 
noisel_0(l:10800,l)=0; 
noise2_0(l:10800,l)=0; 
noise4_0(l:10800,l)=0; 
noise8_0(i:10800,l)=0; 
accel_A0J(l:10800,l)=0; 
accel_AL0(l:10800,l)=0; 
accel_A2_0(l :  10800, 1>=0; 
accel_A4_0(l:  10800,1)-0; 
accetA8_0{l :  10800, 1)=0; 
accel_BOJ(l:10800,l)=0; 
acceLBl_0(l:10800,l)=0; 
accel_B2^0(l:10800,l)=0; 
accel_B4_0(l :  1 0800,  l)-0; 
acceLB8_0(l:10800,l)-0; 
accetnoise0_5(l :  10800, 1  )=0; 
accel_noise  1_0(1 : 1 0800, 1 )=0; 
accel_noise2_0(l :  10800, 1)=0; 
accel_noise4_0(l :  1 0800, 1 )=0; 
accel_noise8_0(l :  10800, 1)=0; 

%  Load  matrices  by  amplitude 
parameter50_5(l:1080,l)-F(l:1080,l); 
parametersO  J(108 1 :2 160, 1  )=F(540 1 :6480, 1); 


%x=10 

%x=20 
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parametersO_5(2161:3240,l)=F(10801:11880,l);  %x=40 

parametere0_5(3241 :4320,1^F(16201:17280,1);  %x=80 

parametereO_5(4321:5400,l)=F(21601;22680,l);  %x=160 

parametersl_0(l:1080,l)=F(1081:2160,l);  %x=10 

parameten;  1_0(108 1 :2 160,1)=F(648 1 :7560, 1);  %x=20 

paramelersl_0(2161:3240,l)=F(11881:12960,l);  %x=40 

parametersl_0(3241:4320,l^F(17281;18360,l);  %x=80 

paranietersl_0(4321:5400,l)=F(22681:23760.1);  %x=160 

parametere2_0(l:1080,l)=F(2161;3240,l);  %x=10 

parameters2_0(108 1  ;2 160,1)=F(756 1 :8640, 1);  %x=20 

parameters2_0(2161:3240,l)=F(12961:14040,l);  %x=40 

pararoeters2_0(3241:4320,l)=F(18361:19440,l);  %x=80 

parainetere2_0(4321;5400,l)=F(23761:24840,l);  %x=160 

parameters4_0(l:1080,l)=F(3241:4320,l);  %x=10 

parainetere4_0(1081:2160,l)=F(8641:9720,l);  %x=20 

parajnetCTs4_0(2161:3240,l)=F(14041:15120,l);  %x=40 

parameters4_0(3241;4320,l)=F(19441:20520,l);  %x=80 

parameters4_0(4321:5400,l)=F(24841:25920,l);  %x=160 

parameters8_0(l:1080,l)=F(4321:5400,l);  %x=10 

paraineters8_0(l 08 1:21 60, 1)=F(972 1:10800,1);  %x=20 

parainetets8_0(2 1 6 1 :3240, 1)=F(1 5 12 1 : 16200, 1);  %x=40 

paraiiieters8_0(3241:4320,l)=F(20521:21600,l);  %x=80 

parameteis8_0(432 1 :5400,1)=F(25921 :27000,1);  %x=160 

fimd_hannO_5(l:4320,l)=G(l:4320,l);  %x=10 

fund_hann0_5(432 1 :8640,1)=G(21601 :25920,1);  %x=20 

lund_hann0_5(8641:12960,l)=G(43201:47520,l);  %x=40 

fund_hann0  5(12961:17280,1)=G(64801:69120,1);  %!c=80 

fimd_hannOl5(17281:21600,l)=G(86401:90720,l);  %x=160 

fund_hannl_0(l:4320,l)=G(4321:8640,l);  %x=10 

fund_harml_0(4321:8640,l)=G(25921:30240,l);  %x=20 

fund_hannl_0(8641:12960,l)=G(47521:51840,l);  %x=40 

fund_hannl_0(12961:17280,l)=0(69121:73440,l);  %x=80 

fund_harml_0(17281:21600,l)=G(90721:95040,l);  %x=160 

fund_haim2_0(l  :4320, 1)=G(8641 : 12960, 1);  %x=10 

fund_harm2_0(4321:8640,l)=G(30241:34560,l);  %x=20 

fund_harm2_0(8641:12960,l)=G(51841:56160,l);  %x=40 

fund_hann2_0(12961:17280,l)=G(73441:77760,l);  %x=80 

£und_harm2  0(17281:21600,1)=G(95041:99360,1);  %x=160 

fund_hann4_0(l  :4320, 1)=G(1296 1 : 17280, 1);  %!c=10 

fiind_harm4_0(4321:8640,l)=G(34561:38880,l);  %x=20 

fund_hann4_0(8641:12960,l)=G(56161:60480,l);  %x=40 

fiind_harm4_0(12961:17280,l)=G(77761:82080,l);  %x=80 

fund_hann4_0(17281:21600,l)=G(99361:103680,l);  %x=160 

fund_harm8_0(l:4320,l)=G(17281:21600,l);  %x=10 

lund_hann8_0(4321:8640,l)=G(38881:43200,l);  %x=20 

fimd_haim8_0(8641:12960,l)=G(60481:64800,l);  %x=40 

£und_haim8_0(12961:17280,l)=G(82081:86400,l);  %x=80 

fund_harm8_0(17281:21600,l)=G(103681:108000,l);  %x=160 
noise0_5(l:2160,l)=H(l:2160,l);  %x=10 

noise0_5(2 161 :4320, 1)=H(1080 1 : 12960, 1);  %x=20 

noise0_5(4321:6480,l)=H(21601:23760,l);  %x=40 

noise0_5(6481 :8640, 1)=H(3240 1:34560, 1);  %x=80 

noise0_5(8641 : 10800, 1)=H(4320 1 :45360, 1);  %x=160 

noisel_0(l:2160,l)=H(2161:4320,l);  %x=10 

noisel_0(2161:4320,l)=H(12961:15120,l);  %x=20 

noisel_0(4321:6480,l^H(23761:25920,l);  %x=40 

noisel_0(6481 :8640, 1)=H(3456 1 :36720, 1);  %x=80 

noisel_0(8641:10800,l)=H(45361:47520,l);  %x=160 

noise2_0(l:2160,l)=H(4321:6480,l);  %xf10 

noise2_0(2161:4320,l)=H(15121:17280,l);  %x=20 

noise2_0(432 1 :6480, 1)=H(2592 1 :28080, 1);  %x=40 

noise2_0(6481:8640,l)=H(36721:38880,l);  %x=80 

noise2_0(8641:10800,l)=H(47521:49680,l);  %x=160 

noise4_0(l:2160,l)=H(6481:8640,l);  %x=10 

noise4_0(2161:4320,l)=H(17281:19440,l);  %x=20 

noise4_0(4321:6480,l)=H(28081:30240,l);  %x=40 
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noise4^0(648 1 :8640, 1  )=H(3 888 1:41 040, 1 ); 
noise4_0(864 1 : 1 0800, 1)=H(4968 1 :5 1 840, 1 ); 
noise8_^0(l:2160,l)=H(8641:10800,l); 
noise8_0(2 16 1 :4320, 1 )=H( 1944 1 :2 1 600, 1 ); 
noise8_0(432 1 :6480, 1 )=H(3024 1 :32400, 1 ); 
noise8_0(6481:8640,l>=H(41041:43200,l); 
noise8_0{864 1 : 10800, 1)=H(5 1 841 :54000, 1 ); 
acceLA0_5(l:2160,l)=J(l:2160,l); 
acce!_A0J(2161:4320,l)=J(10801:12960,l); 
accel_A0J(4321:6480,l)=J(21601:23760,l); 
accel_A0_5(648 1 :8640, 1 )=J(3240 1 :34560, 1 ); 
accel_A0_5(8641 : 10800, 1)=J(4320 1 :45360, 1); 
accel_Al_0(l  :2 160, 1  )=J(2 1 6 1 :4320, 1 ); 
acceLAl_0(2161;4320,l)=J(12961:15120,l); 
accel_Al_0(4321:6480,l)=J(23761:25920,l); 
accel_Al_0(6481:8640,l)=J(34561;36720,l); 
accel_A1^0(864 1 : 10800, 1 )-J(4536 1 :47520, 1); 
accel_A2_0(l  :2 160,1)=J(432 1 :6480, 1); 
accel_A2_0(2161 :4320,1)=J(1 5 121 : 17280, 1); 
accel_A2_0(432 1 :6480,1)=J(25921 :28080, 1); 
accel_A2_0(6481:8640,l)=J(36721:38880,l); 
accel_A2_0(864 1 : 10800, 1 )=J(4752 1 :49680, 1 ); 
accel_A4_0(l  :2 160, 1)=J(648 1 :8640, 1 ); 
accel_A4_0(2161:4320,l)=J(17281:19440,l); 
accel_A4_0(4321:6480,l)=J(28081:30240,l); 
accel_A4_0(6481:8640,l)=J(38881:41040,l); 
accel_A4_0(8641 : 1 0800, 1 )= J(4968 1:51 840, 1 ); 
accel_A8__0(l  :2 160, 1  )=J(864 1 : 1 0800, 1 ); 
accel_A8_0(2161:4320,l)-J(19441:21600,l); 
accel_A8_^0(432 1 :6480, 1 )=J(3024 1 :32400, 1 ); 
accel_A8^0(648 1 :8640, 1 )= J(4 1041:43200,1 ); 
acce!_A8_0(8641: 10800, 1  )=J(5 1 841 :54000, 1); 
accel_B0_5(l  :2 160, 1  )=H(  1 :2 160, 1 ); 
accel_B0  J(2 16 1 :4320, 1 )=K(  1 0 80 1 : 12960, 1 ); 
accel_B0_5(432 1 :6480, 1 )=K(2 1 60 1 :23760, 1 ); 
accel_B0_5(648 1 : 8640, 1 )=K(3240 1 :34560, 1 ); 
accel_B0_5(8641 : 10800, 1 )=K(4320 1 :45360, 1); 
accel_B  1_0(1 :2 1 60, 1  )=K(2 16 1 :4320, 1); 
accel_Bl_0(2161:4320,l)=K(12961:15120,l); 
accel_Bl_0(4321:6480,l)=K(23761:25920,l); 
accel_Bl_0(648 1 :8640, 1)=K(3456 1 :36720, 1); 
accel_Bl_0(8641:10800,l)=K(45361:47520,l); 
accel_B2_0(l:2160,l)=K(4321:6480,l); 
accel_B2_0(2 16 1 :4320, 1  )=K(  15121 : 17280, 1); 
accel_B2_0(432 1 :6480, 1 )=K(2 592 1 :28080, 1); 
acceLB2_0(648 1 : 8640, 1 )=K(3672 1:38880,1); 
accel_B2_0(8641:10800,l)=K(47521:49680,l); 
accel_B4_0(  1 :2 1 60, 1  )=K(648 1 :8640, 1 ); 
4cceLB4_0(2161:4320,l)=K(17281:19440,l); 
accel_B4_0(432 1 :6480, 1)=K(2808 1 :30240, 1); 
acceI_B4__0(648 1 :8640, 1)-K(3888 1 :4 1040, 1); 
acceLB4_0(8641:10800,l)=K(49681:51840,l); 
accel_B8_0(l  :2160,1)=K(8641 : 10800,1); 
accel_B8_0(2 161 :4320, 1)=K(1944 1 :2 1600, 1); 
acceLB8_0(432 1 :6480, 1 )=K(3024 1 :32400, 1); 
accel^B8_0(6481 :8640, 1)=K(41041 :43200,1); 
accel_B8^0(8641:10800,l)=K(51841:54000,l); 
accel_noiseO  J(1 :2 160, 1)-L(1 :2 160, 1 ); 
aocel^noiseO  J(2 16 1 :4320, 1)-L(1080 1 .12960, 1 ); 
accel_noise0_5(432 1 :6480, 1)=L(2 160 1 :23760, 1); 
accel_noise0_5(6481:8640,l)-L(32401:34560,l); 
acce!_noise0J(8641:10800,l)=L(43201:45360,l); 
acceLnoisel_,0(l  :2 160, 1)=1X2 1 61 :4320, 1); 
acceLnoisel  0(2161:4320,1)-L(12961:15120,1); 
accel_noiser0(432 1 :6480, 1)=L(2376 1 :25920, 1); 
accel_noise  1  _0(648 1 :8640, 1  )=L(3456 1 :36720, 1 ); 


%x=80 

%x=160 

%x=10 

%x=20 
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accel_noisel_0(8641:10800,l)=L(45361:47520,l);  %x=160 

accel_noise2_0(  1 :2 160, 1)=L(432 1 :6480, 1);  %x=  10 

accel_noise2_0(2161:4320,l)=L(15121:17280,l);  %x=20 

accel_noise2_0(4321 :6480, 1>=L(25921 :28080,l);  %x=40 

acceLnoise2_0(6481:8640,l)=L(36721:38880,l);  %x=80 

acceLnoise2^0(8641:10800,l)==L(47521:49680,l);  %x=160 

acce!_noise4_0(l:2160,l)=L(6481:8640,l);  %x=10 

accel_noise4^0(2161:4320,l)=L(17281:19440,l);  %x=20 

accel_noise4_0(4321:6480,l>=L(28081:30240,l);  %x=40 

accel_noise4_0(6481:8640,l)=L(38881:41040,iy,  %x=80 

accel_noise4_0(864 1 : 10800, 1)=L(4968 1:51 840, 1 );  %x=l  60 

accel_noise8_0(l:2160,l)=L(8641:10800,l);  %x=10 

acceLnoise8_0(2161:4320,l)=L(19441;21600,l);  %x=20 

acceLnoise8_0(432 1 :6480,1)=L(30241 :32400,1);  %x-40 

accel_noise8_0(6481:8640,l)=L(41041:43200,l);  %x=80 

accel_noise8_0(8641:10800,l)-L(51841:54000,l);  %x=160 


%  Save  matrices  as  .mat  files 
save  parameters0_5 1  parameters0_5; 
save  parameters  1_0 1  parameters  l_0; 
saveparameters2_01  parameters2_0; 
save  parameters4_01  parameters4_0; 
save  parameters8_01  parameters8_0; 
save  fund_harm0_5 1  fond_harm0_5; 
save  fund_harm  1  _0 1  fiind_harm  1  _0; 
save  fiind_harm2  01  fiind_harm2_0; 
save  fund_harm4_01  fund_harm4_0; 
save  fundharm8_0 1  fundharm8_0; 
save  noise0__51  noise0_5; 
save  noisel  Ol  noisel  O; 
save  noise2_01  noise2_0; 
save  noise401  noise4_0; 
save  noise8_01  noise8_0; 
save  accel_A0_51  accel_A0_5; 
save  accel_Al_01  acceI_Al_0; 
save  accel_A2_0 1  accel_A2_0; 
save  accel_A4_01  accel_A4_0; 
save  accel_A8_01  accel_A8_0; 
save  accel_B0_51  accel_B0_5; 
saveacce!_Bl_01  accel_Bl_0; 
save  accel_B2_01  accel_B2_0; 
save  accel_B4_01  accel_B4_0; 
saveaccel_B8_01  accel_B8_0; 
save  accel_noise05 1  accel_noise05; 
save  accel_noisel_01  accel  noise  1_0; 
save  accel_noise2_01  acce!_noise2_0; 
saveaccel_noise4_01  accel_noise40; 
save  acce!noise8_01  accel_noise8_0; 
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Program  2 


%  This  program  will  ask  for  an  amplitude  and  call  up  the  appropriate  transfer  function  data, 

%  It  will  then  plot  the  five  different  positions  for  that  amplitude  on  the  same  graph. 

dispLfreqlO_fund(l :  1 80,1  )=0; 
displ_freq20_fund(  1:180,1  )=0; 
displ_freq40_fund(  1:180,1  )=0; 
displ_freq80_fund(  1 : 1 80, 1  )=0; 
displ_freql  60_fund(  1 : 1 80, 1  )=0; 
freql0(l:180,l)=0; 
freq20(l:180,l>=0; 
freq40(l:180,l>=0; 
freq80(l:180,l)=0; 
freql  60(1:1 80, 1)=0; 

if  existCiteration')==0 

iteration=mput('lst  or  2d  iteration  (type  "1"  or  "2’*):  ','s'); 
end; 

ifiteration=-r 
if  exist(’amplitude')=0 

amplitude=ii^)ut(’ Amplitude  (enter  3  char): 
end; 

if  amplitude=- 0.5’ 

load  parametersOS  1 ; 
load  fund_harm0_5 1 ; 
load  noise05 1; 
load  accel_A0_51; 
load  accel_B0_51; 
load  accel_noise0_5 1 ; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
VoFUNDAMENTAL 

forloop2=0:179  x 

displ  freql0_fund((l+loop2),l)=fund_harm0_5((13+loop2»24),l)+i*fund_hann0_5((14+loop2*24),l); 

displ  freq20_fund((l+loop2Xl)=fund_harm0_5((4333+loop2*24),l)+i*fund_hann0_5((4334+loop2*24),l); 
displ  freq40_fiind((l+loop2),l)=lund_hann0j((8653+loop2*24),l)+i*fiind_hann0_5((8654+loop2*24),l); 

displ  freq80_fund((l+loop2),l)=fiind_hann0_5((12973+loop2»24),l)+i*fund_harm0_5((12974+loop2*24),l); 
displ_freql60_fund((l+loop2),l)=fund_hann0_5((17293+loop2*24),l)+i*fund_hann0_5((17294+loop2*24),l); 

end 

%  Write  out  frequencies  recorded 
freql  0(1 : 1 80,  l)=parameters0_5 ((3 : 6 : 1 080),  1 ); 
freq20(l:  180,l)=parameters0_5((1083:6:2 160),1); 
freq40(l :  180,  l)=parameters0_5((2 1 63:6:3240),  1); 
freq80(l :  180,l)=parameters0_5((3243:6:4320),l); 
freql60(l:180,l)^>arameters0_5((4323:6:5400),l); 
elseif  amplitude='  1 . 0’ 
load  parameters  10 1 ; 
load  fund_harml_01; 
load  noisel_01; 
load  accel_Al_01; 
load  accel_B  1_0 1 ; 
load  accel_noisel_0 1; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
%FUNDAMENTAL 

forloop2=0:179  . 

displ_freql0_fund((l+loop2),l)=fund_hannl_0((13+loop2»24),l)+i»aind_hannl_0((14+loop2*24),l); 
displ  freq20  fund((l+loop2),lHund_hannl_0((4333+loop2*24),l)+i*fund_hannl_0((4334+loop2*24),l); 
displ  freq40''fund((l+loop2),l)=fund_hannl_0((8653+loop2*24),l)+i*lund_harml_0((8654+loop2»24),l); 
displ  freq80  fund((l+loop2),l)=fund_haiml_0((12973+loop2»24),l)+i*fund_hannl_0((12974+loop2*24),l); 
disprfreql60_fund((l+loop2),l)=lund_hamil_0((17293+loop2*24),l)+i*fund_harml_0((17294+loop2*24),l); 

end 

%  Write  out  amplitudes  recorded 
fi^ql0(l :  180,l)F=parameters  1_0((3 :6: 1080),  1); 
freq20(l :  1 80, 1  )=parameters  1  _0(( 1083:6:21 60),  1 ); 
freq40(l :  1 80,1  )=parametersl_0((2 163 :6:3240),1); 
freq80( 1 : 1 80, 1 )=parameters  1_0((3243 :6:4320),  1 ); 
freql60(l:180,l)=parametersl_0((4323:6:5400),l); 
elseif  amplitude=- 2.0’ 
load  paramelers2_0 1; 
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load  fund_harm2_0 1 ; 
load  noise2_01; 
load  accel_A2_0 1 ; 
load  accel_B2_01; 
load  accel_noise2_0 1 ; 

%  Combine  real  and  imaginaiy  parts  of  displ_freq 

%FUNDAMENTAL 

for  loop2-0:179 

displ_freq  1 0_fund((  1 +loop2),  1  )-fund_harm2_0((l  3  +loop2*24),  1  )+i*fund  hann2_0((  1 4+loop2*24),  1 ); 

displ_fi‘eq20_fiind((  1  +loop2),  1  )=fund_harm2_0((4333 +loop2*24),  1  )+i*iund_harm2_0((4334+loop2*24),  1); 

displ_freq40^fund((l+loop2),lHund_harm2_0((8653+loop2*24),l)+i*fiind_hann2_0((8654+loop2*24),l); 

displ_freq80_fund((l+loop2Xl)==lund_harm2_^0((12973+loop2*24),l)+i*fiind_harTn2_0((12974+loop2*24),l); 

displ_freql60_fund((l+loop2),l)=fund_hanTi2_0((17293+loop2*24),l)+i»fiind_hann2_0((17294+loop2*24),l); 

end 

%  Write  out  amplitudes  recorded 
freqlO(l :  1 80,  l)=parameters2_0((3 :6: 1080),  1); 
fireq20( 1 : 1 80, 1 )=parameters2_0(( 1 0 83 :6 :2 1 60),  1 ); 
freq40( 1:180,1 )=parameters2_0((2 163:6:3 240),  1 ); 
freq80( 1:180,1 )=parameters2_0((3243  r6 :4320),  1 ); 
freql60(l:180,l)^)arameters2_0{(4323:6:5400),l); 
elseif  amplitude=-4.0' 
load  parameters4_0 1 ; 
load  fund_harm4_0 1 ; 
load  noise4_01; 
load  accel_  A4_0 1 ; 
load  accel_B4_01; 
load  accel_noise4_01; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

for  loop2=0:179 

displ_freq  1 0_fund((  1  +loop2),  1  )=fimd_harm4_0((  1 3 +loop2  *24),  1  )+i*fiind_harm4_0((  1 4+loop2  *  24),  1 ); 

displ_freq20_fund((l+loop2),l)=fiind_harm4_0((4333+loop2*24),l)+i*fund_harm4_0((4334+loop2*24),l); 

displ_freq40_fund((  1  +loop2),  1  )=ftind_harm4_0((86 53+loop2  *24),  l)+i*fundharm4_0((8654+loop2*24),  1); 

displ_freq80_fund((l+loop2y)=fund_harm4_0((12973+loop2*24),l)+i*fund__harm4_0((12974+loop2*24),l); 

displ_freql60_fun^(l+loop2),l)=flind_harm4_0((17293+loop2*24),l)+i*fund_harm4_0((17294+loop2*24),l); 

end 

%  Write  out  amplitudes  recorded 
freql0(l  :180,l)^)arameters4_0((3:6:1080),  1); 
freq20( 1:180,1 )=parameters4_0(( 1083:6:2160),!); 
freq40(l :  1 80,  l)=parameters4_0((2 163 :6:3240),  1); 
freq80(l :  1 80,l)=parameters4_0((3243 :6:4320),  1^, 
freql60(l:180,l)=parameters4_0((4323:6:5400),l); 
elseif  amplitude^- 8.0’ 
load  parameters8_0 1 ; 
load  fund_harm8_0 1 ; 
load  noise8_01; 
load  accel_  A8_0 1 ; 
load  accel_B8_01; 
load  accel_noise8_01; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

forloop2==0:179 

displ_freql0_fund((l+loop2),l)=fund_harm8_0((13+loop2*24),l)+i*fund_harm8_0((14+loop2*24),l); 

displ_freq20_fund((  1  +loop2),  1  )=fund_harm8_0((4333+Ioop2  *24),  l)+i*fund_harm8_0((4334+loop2*24),  1); 

displ_freq40^fund((l+loop2),l)-fund_harm8_0((8653+loop2*24),l)+i*fund^hann8_0((8654+ioop2*24),l); 

displ_freq80Jiind((l+loop2y)=fund_harm8_0((12973+loop2*24),l)+i*fund_harm8_0((12974+loop2*24),l); 

displ_freql60_fund((l+loop2),l)=fund_harm8_0((17293+loop2*24),l)+i*fund_hann8_0((17294+loop2*24),l); 

end 

%  Write  out  amplitudes  recorded 
freql  0(  1 : 1 80, 1  )=parameters8_0((3 :6: 1 080),  1 ); 
freq20(l:180,l)=parameters8_0((1083:6:2160),l); 
freq40(l :  1 80,  l)=parameters8_0((2 163 :6:3240),iy, 
frcq80(l :  1 80, 1  )=parameters8_0({3243 :6:4320),  1 ); 
freql60(l:180,l)^)arameters8_0((4323:6:5400),l); 
end; 

elseif  iteration==7’ 

if  existCamplitude')==0 

ainplitude=inputC Amplitude  (enter  3  char):  ','s'); 
end; 
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if  amplitude=-0.5' 
load  parameters0_52; 
load  fund_hami0_52; 
load  noiseO_52; 
load  accel_A0_52; 
load  accel_B0_52; 
load  acce!_noise0_52; 

%  Combine  real  and  imaginajy  parts  of  displ  freq 

%FUNDAMENTAL 

for  loop2=0:179 

displ_freql0_fund((l+loop2),l)=fund_harm0_5((13+loop2*24),l)+i*fund_hann0_5((14+loop2*24),l); 

displ_freq20_fund((  1  +loop2),  1  )=fund_harm0_5((4333+loop2*24),  1  )+i*fund_hann0_5((4334+loop2*24),  1 ); 

displ_freq40_fiind((l+loop2),l)=fiind3ann0_5((8653+loop2*24),l)+i*fund_hann0_5((8654+loop2*24),l); 

displ_freq80_fund((l+loop2),l)=fund_harm0_5((12973+loop2*24),l)+i*fund_harm0_5((12974+loop2*24),l); 

displ_freql60_fun^(l+loop2),l)=fund_harm0_5((17293+loop2*24),l)+i*fund_harm0_5((17294+loop2*24),l); 

end 

%  Write  out  frequencies  recorded 
freql0(l:180,l)=parameters0_5((3:6:1080),l); 
freq20(  1 ;  1 80, 1  )=parameters0_5(( 1083:6:21 60),  1); 
freq40(l :  1 80,  l)=parameters0_5 ((2 163:6:3240),  1); 
freq80(l :  1 80,  l)=parameters0_5((3243 :6:4320),  1); 
freql60(l:180,l)=parameters0_5((4323:6:5400),l); 
elseif  amplitude^- 1.0' 
load  parameters  1_0  2; 
load  fund_harml_02; 
load  noisel_02; 
load  accel_Al_02; 
load  accel_Bl_02; 
load  accel_noisel_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

for  loop2-0: 179 

displ_freql0_fund((l+loop2),l)=fund_harml_0((13+loop2*24),l)+i*fund_harml_0((14+loop2*24),l); 

disp!_freq20_flind((  1  +loop2),l)=fund_harm  l_0((4333+loop2*24),  l)+i*fund_hann  1  _0((4334+loop2*24),  1); 

displ_freq40_fund((l+loop2),l)=fund_harml_0((8653+loop2*24),l)+i*fund_harml_0((8654+loop2*24),l); 

displ_freq80_fund((l+loop2),l)=fund_hannl_0((12973+Ioop2*24),l)+i*fUnd_harml_0((12974+Ioop2*24),l); 

displ__freql60_fund((l+loop2),l)=fund_harml_0((17293+loop2*24),l)+i*fiind_harml_0((17294+loop2*24),l); 

end 

%  Write  out  amplitudes  recorded 
freq  10(1:180,1  )=parameters  1_0((3 :6 : 1 080),  1 ); 
freq20(l:180,l)=parametersl_0((1083:6:2160),l); 
freq40( 1:180,1 )=parameters  1  _0((2 1 63 :6 :3 240),  1 ); 
freq80( 1 : 1 80, 1  )=parametersl_0((3243 :6:4320),  1 ); 
freq 160( 1 : 1 80, 1  )=parameters  1  _0((4323 :6 : 5400),  1 ); 
elseif  amplitude— 2.0' 
load  parameters2_02; 
load  fund_harm2_02; 
load  noise2_02; 
load  accel_A2_02; 
load  accel_B2_02; 
load  accel_noise2_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

forloop2=0:179 

displ_freql0_fund((l+loop2),l)=fund_harm2J((13-floop2*24),l)+i*fund_harm2_0((14+loop2*24),l); 

displ_freq20_^fund((l+loop2),l)=fund_harm2^0((4333+loop2*24),l)+i*fund_harm2_0((4334+loop2*24),l); 

displ_freq40_fiind((l +loop2),  1  )=fiind_harm20((865 3+loop2*24),  1  )+i*fimd_harm2_0((8654+Ioop2*24),  1 ); 

displ_freq80_fund((l+loop2),l>=fund_harm2_0((12973+loop2*24),l)+i*fund_harm2^0((12974+loop2*24),l); 

displ_freql60_jaind((l+Ioop2),l)=fiind_hann2_0((17293+loop2*24),l)+i*fund_hann2_0((17294+loop2*24),l); 

end 

%  Write  out  amplitudes  recorded 
freql0(  1 : 1 80, 1  )=parameters2_0((3 :6 : 1 080),  1 ); 
freq20(  1:180,1 )=parameters2_0(( 1083:6:2160),!); 
freq40( 1:180,1 )=parameters2_0((2 163:6 :3240),  1 ); 
freq80(l :  1 80,  l)=parameters2_0((3243 :6:4320),  1 ); 
freql60(l:180,l)=parameters2_0((4323:6:5400),l); 
elseif  amplitude^- 4.0* 
load  parameters4_02; 
load  fimd_hann4_02; 
load  noise4_02; 
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load  accel_A4_02; 
load  accel_B4_02; 
load  accel_noise4_02; 

%  Combine  real  and  imaginary  parts  of  displ  freq 

%FUNDAMENTAL 

for  loop2=0:179 

displ_freql0_lund((l+loop2),l)=fund_harm4_0((13+loop2*24),l)+i*fiind_harm4_0((14+Ioop2*24),l); 
displ_freq20_fund((l+loop2),l)=fiind_hann4_0((4333+loop2*24),l)+i*fund_hann4_0((4334+loop2*24),l); 
displ_freq40_fiind((l+loop2),l)=fund_hann4_0((8653+Ioop2*24),l)+i*fund_hann4_0((8654+loop2*24),l); 
displ  freq80_lund((l+loop2),l)=fund_hann4_0((12973+loop2*24),l)+i»lund_hann4_0((12974+loop2»24),l); 

displ_freql60_fUnd((l+loop2),l)=fund_haim4_0((l7293+loop2*24),l)+i»fund_hann4_0((17294+loop2»24),l); 

end 

%  Write  out  amplitudes  recorded 
freq  1 0(  1 : 1 80, 1  )=parameters4_0((3 :6 : 1 080),  1); 
freq20(l :  180,  l)=parameters4_0((1083:6:2 160),1); 
freq40( 1 : 1 80, 1 )=parameters4_0((2 1 63 :6 :3240),  1); 
freq80(l :  1 80,l)=parameters4_0((3243 :6:4320),  1); 
freq 1 60(1 : 1 80,1 )=parameters4_0((4323 :6 :5400),  1 ); 
elseif  amplitude=- 8.0' 
load  parameters8_02; 
load  fiind_harm8_02; 
load  noise8_02; 
load  accel_A8_02; 
load  accel_B8_02; 
load  accel_noise8_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

for  loop2=0:179 

displ_freql0_fund((l+Ioop2),l)=ftind_hann8_0((13+loop2*24),l)+i*fund_harm8_0((14+loop2*24),l); 
displ_freq20_lund((l+loop2),l)=fund_hann8_0((4333+loop2»24),l)+i»fund_hann8_0((4334+Ioop2*24),l); 
displ  freq40Jund((l+loop2il)=fund_hann8_0((8653+loop2*24),l)+i»fund_hann8_0((8654+loop2»24),l); 
displ  freq80_fund((l+loop2),l)=fund_hann8_0((12973+loop2*24),l)+i*fund_harm8_0((12974+loop2*24),l); 
displ_freql60_fund((l+loop2),l)=fund_hann8_0((17293+loop2»24),l)+i*fund_hann8_0((17294+loop2*24),l); 
end 

%  Write  out  amplitudes  recorded 

freq  1 0(  1 : 1 80,  l)=parameters8_0((3 :6: 1 080),  1 ); 

freq20( 1 : 1 80, 1 )=parameters8_0((  1 083 :6:2 1 60),  1 ); 

freq40( 1 : 1 80, 1 )=parameters8_0((2 1 63 :6:3240),  1 ); 
freq80(l:180,l)=parameters8_0((3243:6:4320),l); 
freql60(l:180,l)=parameters8_0((4323:6:5400),l); 
end; 
end; 

%  Plot  results 
figure(l) 

semilogy(freq  1 0,abs(displ_freq  1  Ofixnd),'-') 
hold  on 

seniilogy(freq20,abs(displ_freq20_fund),’--’) 

semilogy(freq40,abs(displ_freq40_fiind),’-.’) 

semilogy(freq80,abs(displ_freq80_fund),':') 

semilogy(freq  1 60,abs(displ_freql  60_fund),'-') 

title  l='Displacement  vs  Frequency'; 

title2-  Amplitude  = 

title2“Strcat(title2,amplitude,'Volts'); 

titles-  Iteration  =  '; 

title3=strcat(title3,iteration); 

title jJata-char(  {title  I,title2,title3  }); 

title(title_data); 

ylabel(T)isplacement') 

xlabelCFrequency  (Hz)*) 

legendflO  cm','20  cm', '40  cm', '80  cm','160  cm') 

hold  off 

orient  landscape 

figure(2) 

plot(freq  1 0,abs(displ_freq  1 0_fund),'-') 
hold  on 

plot(freq20,abs(displ_freq20_fund),'--') 

plot(freq40,abs(displ_freq40_fund),'-.’) 

plot(freq80,abs(displ_freq80_fund),':') 
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plot(freq  1 60,abs(ciispl_freq  1 60_fiind),'- ) 

title  1- Displacement  vs  Frequency*; 

title2-  Amplitude  =  '; 

title2=strcat(title2, amplitude, ’V  cits’); 

titles-  Iteration  =  '; 

title3=strcat(title3,iteration); 

title_data=char(  {title  1  ,title2,title3  }); 

title(title_data); 

ylabelCDisplacement*) 

xlabelCFrequency  (Hz)*) 

legendlciO  cmV20  cm’, ’40  cm','80  cm’,’160  cm') 

hold  off 

orient  landscape 

figure(3) 

loglog(freq  1 0,abs(displ_fi*eq  1  Ofund),'-*) 
hold  on 

Ioglog(freq20,abs(displ_freq20_fund),'--’) 

loglog  (freq40,abs(displ_freq40  fund),’-.') 

loglog  (freq80,abs(displ_freq80_fund),':') 

loglog  (freql  60,abs(displ_freql  60_fund),'-') 

title  1- Displacement  vs  Frequency; 

title2-  Amplitude-'; 

title2^strcat(title2,amplitude,'Volts’); 

titles-  Iteration  = 

titles  =strcat(title3,  iteration); 

titledata=char({title  1  ,title2,title3  }); 

title(title_data); 

ylabelCDisplacemenf) 

xlabelfFrequency  (Hz)’) 

legendi^lO  cm','20  cm’,’40  cm',’80  cm',' 160  cm’) 

hold  off 

orient  landsc^ 


%  This  program  will  ask  for  an  amplitude  and  a  position  and  call  up  the  appropriate 
%  transfer  function  data.  It  will  then  plot  the  transfer  function  for  that  location 
%  and  amplitude  showing  the  fundamental  and  five  harmonics  on  the  same  graph, 

%  Initialize  matrices 
disptfreql  0_fund(  1:180, 1)=0; 
displ_fi^20^fund(  1 : 1 80, 1)=0; 
displ_fi-eq40lfiind(l :  1 80, 1)=0; 
displ_fi-eq80  fund(l:180,l)=0; 
displ_freql60_fund(l :  180,1)=0; 
displ_fi'eql  0__harml(  1 : 1 80, 1  )=0; 
displ_fi^20_harml  ( 1 : 1 80, 1  )=0; 
displ_freq40_harml(l :  180, 1)=0; 
displ  Jfreq80_hann  1(1:1 80, 1)=0; 
displ_fi"eql  60_harm  1(1:1 80, 1)=0; 
displ_fi‘eql0_hann2(l :  180,1)=0; 
displ_fi-eq20_harm2(  1 ;  1 80, 1)=0; 
displ_fi’eq40_hami2(  1 : 1 80, 1  )=0; 
displ_fi-eq80_harm2(  1 : 1 80, 1  )=0; 
displ_fi-eql60_harm2(l :  1 80, 1)=0; 
displ_fi-eql0_harm3(l :  1 80, 1)=0; 
displ_fi-eq20_harm3(  1 : 180, 1)=0; 
displ_freq40_harm3(l :  1 80, 1  )=0; 
displ_fi‘eq80_harm3(l :  1 80, 1)=0; 
displ_fi‘eq  1 60_harm3(l :  1 80, 1  )=0; 
displ_freql0_harm4(l :  1 80, 1)=0; 
displ_fi^20_hann4(l :  1 80, 1)=0; 
dispLfi‘eq40_hann4(  1 : 1 80, 1  )=0; 
displ_fi-eq80_hann4(l :  1 80, 1)=0; 
displ_fi-eql60_harm4(l :  1 80, 1)=0; 
displ_fi'eql0_harm5(l :  1 80, 1)=0; 
dtspl_freq20_harm5(l :  1 80, 1)=0; 
displ_fi^eq40_harm5(  1 : 1 80, 1  )=0; 
displ_freq80_harm5(  1 : 1 80, 1  )=0; 
displ_fi‘eql60_harm5(l :  180,1)=0; 
fi‘eql0(l:180,l)=0; 
fi-eq20(l:180,l)=0; 
fi-eq40(l:180,l)=0; 
fi-eq80(l:180,l)=0; 
fi‘eql60(l:180,l)-0; 

if  exist(’iteration')==0 

iteration=input('lst  or  2d  iteration  (type  "1"  or  "2"):  Vs'); 
end; 

if  iteration=- r 
if  exist(’amplitude')==0 

amplitude=input(* Amplitude  (enter  3  char):  Vs'); 
end; 

if  amplitude=-0.5' 

load  parameters0_5 1 ; 
load  fund_harm0_5 1 ; 
load  noiseOSl; 
load  accel_A0_5 1 ; 
load  accel_B0_51; 
load  accel_noise0_5 1; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

forloop2=0:179 

(Iispl_freql0_fund((l+loop2),l)=fund_hann0_5((13+loop2*24),l)+i*fund_harm0_5((14+loop2*24),l); 

displ_freq20_fund((l+loop2y^nd_hann0_5((4333+loop2*24),l)+i*£und_harm0_5((4334+loop2»24),l); 

displ_freq40_fund((l+loop2),lHund_hann0_5((8653+Ioop2*24),l)+i*fiind_hann0_5((8654+loop2»24),l); 

displ_freq80_fund((l+loop2),l)=fund_harm0J((12973+loop2*24),l)+i»fund_harm0_5((12974+loop2»24),l); 

displ_freql60_fund((l+loop2),l)=fund_hann0_5((17293+loop2»24),l)+i*fund_hann0_5((17294+loop2*24),l) 
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end 

%HARMONIC  1 


for  Ioop2=0:179 

displ  freql0_haniil({l+loop2),l)=fiind_hann0_5((15+loop2*24),l)+i*fiin(i_harm0_5((16+loop2*24),l); 
displ_freq20_haiml((l+loop2),l)=fund_harm0_5((4335+loop2*24),l>+i*fund_harm0_5((4336+loop2*24),l); 
displ  freq40  hannl((l+loop2),l)=ftind_hann0_5((8655+loop2*24),l)+i*fund_harm0_5((8656+loop2*24),l); 
displ  fteq80  hamil((l+loop2),l)=fund_hann0_5((12975+loop2*24),l)+i*fund_harm0_5((12976+loop2*24),l); 

displ_freql60_hannl((l+loop2),l)=fund_hann0_5((17295+loop2*24),l)+i‘fund_harm0_5((17296+loop2»24),l); 

end 

%HARMONIC  2 
for  Ioop2=0:179 

displ  freqlO  hann2((l+loop2),l)-fund^hann0_5((17+loop2*24),l)+i*fund_hann0_5((18+loop2»24),l); 
displ  freq20  hann2((l+loop2),l)=fund_harm0_5((4337+loop2*24),l)+i*fund_hann0_5((4338+loop2*24),l); 

displ  freq40  harm2((  1  +loop2), I)=fund_hann0_5((8657+loop2*24), l)+i’’‘fund_hami0_5((8658+loop2*24),l); 

disprfreq80”harm2((l+loop2),lHund_hann0_5((12977+loop2*24),l)+i*fund_hann0J((12978+loop2*24),l); 

displ_freql  60_harm2((  1 +loop2),  1  )=fiind_hann0_5(( 1 7297+loop2*24),  1  )+i*iUnd_harmO_5(( 17298+loop2*24),  1); 
end 

%HARMONIC  3 


for  loop2=0:179 

displ  freql0_hann3((l+loop2),l)=fiind_hann0_5((19+loop2»24),l)+i»fund_hann0_5((20+loop2*24),l); 
displ  freq20_harm3((l+loop2),l)=fund_harm0_5((4339+loop2*24),l)+i*fiind_hami0_5((4340+loop2*24),l); 
displ  freq40  hann3((l+loop2),l)=lund_liarm0_5((8659+loop2*24),l)+i»fiind_hann0_5((8660+loop2»24),l); 
displ  freq80  hann3((l+loop2),l)=fund_harm0_5((12979+loop2»24),l)+i»fund_hann0_5((12980+loop2*24).l 
disprfreql60_hann3((l+loop2),l)=fund_haim0_5((17299+loop2*24),l)+i»fund_harm0_5((17300+loop2*24),l); 

end 

%HARMONIC  4 


,1); 


forloop2=0;179  , 

dUpl  freql0_harm4((l+loop2),l)=ftind_hann0_5((21+loop2*24),l)+i*fund_hami0_5((22+loop2*24),l); 
disprfreq20_hann4{(l+loop2),l)=fund_harm0_5((4341+loop2*24),l)+i*fund_hann0_5((4342+loop2*24),l); 
displ  freq40  harm4((l+loop2),lHiind_hann0_5((8661+loop2*24),l)+i*fund_harm0_5((8662+loop2»24),l); 
displ  freqSO  hann4((l+loop2),l)=fund_harm0_5((12981+loop2*24),l)+i*fund_hann0_5((12982+loop2»24).l); 
displ_freql60_hann4((l+loop2),l)=fund_harm0_5((17301+loop2»24),l)+i»fund_hann0_5((17302+loop2*24),l); 

end 

%HARMONIC  5 


forloop2=0:179 

displ  freql0_harm5((l+loop2),l)=fund_hann0_5((23+loop2»24),l)+i*fund_hann0_5((24+loop2»24),l); 
displ  freq20_hann5((l+loop2),l)=fund_hann0_5((4343+loop2»24),l)+i»£und_hanii0_5((4344+loop2*24),l); 
displ  freq40_hann5((l+loop2),l)=fund_hanti0_5((8663+loop2*24),l)+i*fund_harm0_5((8664+loop2*24),l); 
displ  freqSO  hann5((l+loop2).l)=fund_hann0_5((12983+loop2*24),l)+i*fund_harm0_5((12984+loop2*24).l); 
displ_freql60_hann5((l+loop2),l)=fimd_harm0_5((17303+loop2*24),l)+i*lund_hann0_5((17304+loop2*24),l); 

end 

%  Write  out  frequencies  recorded 
freqlO(l :  1 80,  l)=parameters0_5((3 :6: 1 080),  1); 
freq20(l :  1 80,l)=parameters0_5((l083:6:2 160),  1); 
freq40(l :  180,l)=parameters0_5((2163;6:3240),l); 
freq80(l :  1 80, 1  )=parameters0_5((3243 :6:4320),  1); 
freql60(l:180,l)=parameters0_5((4323:6:5400),l); 
elseif  aniplitude=~  1.0' 
load  parameters  1_01; 
load  fiindharml  _0 1 ; 
load  noisel_01; 
load  accel_Al_01; 
load  accel_Bl_01; 
load  accel_noise  1  _0 1 ; 

%  Combine  real  and  imaginary  parts  of  displ  freq 

VoFUND  AMENT  AL 

forloop2=0:179 

displ  freql0_fund((l+loop2),l)=fund_hannl_0((13+loop2*24),l)+i*fiind_hiirml_0((14+loop2*24),l); 

displ_freq20_fund((l+loop2),l)=fund_haitnl_0((4333+loop2*24),l)+i*fund_hannl_0((4334+loop2»24),l); 

displ_freq40_fund((l+loop2),l)=fimd_hannl_0((8653+loop2*24),l)+i*fund_haniil_0((8654+loop2*24),l); 
displ  freq80_fund((l+loop2),l)=fund_hannl_0((12973+loop2*24),l)+i‘fund_harml_0((12974+loop2»24).l); 
displ_freql60_fund((l+loop2),l)=fund_harml_0((17293+loop2*24),l)+i*fund_harml_0((17294+loop2*24),l); 

end 

%HARMONIC  1 
forloop2’=0:179 

displ_freql0_harml((l+loop2),l)=fund_harml_0((15+loop2*24),l)+i*fund_harml_0((16+loop2*24),l); 
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displ  freq20  harml((l+loop2),l)=lund_hannl_0((4335+loop2*24),l)+i*fund_harml_0((4336+loop2*24),l); 
displ  freq40  hannl((l+loop2),lHund_hannl_0((8655+loop2*24),l)+i»fund_hannl_0((8656+loop2  24), 1); 
displ  freqSO  hannl((l+loop2),l)=fund_harml_0((12975+loop2*24),l)+i*fund_harml_0((12976+loop2*24)  l); 
displ_freql60_hannl((l+loop2),l)=fund_hannl_0((17295+loop2»24),l)+i*fund_hannl_0((17296+loop2  24), 1); 

end 

•/oHARMONIC  2 

dkprfreql0_hann2((l+loop2),l)=fund_hannl_0((17+loop2*24),l)+i»fund_hannl_0((18+loop2»24),l); 

displ  freq20  harm2((l+loop2),l)=fund_harnil_0((4337+loop2*24), l)+i*fund_hannl_0((4338+Ioop2*24),l); 
dispffreqdO  harm2((l+loop2).l)=fund_hannl_0({8657+loop2*24),l)+i*£und_hannl_0((8658+loop2»24),l); 
displ  freqSO  harm2((l+loop2),l)=fund_hannl_0((12977+loop2*24),l)+i*fund_harml_0((12978+loop2*24),l); 
displ_freql60_hann2((l+loop2),l)=fund_hannl_0((17297+loop2*24),l)+i»fund_hannl_0((17298+loop2»24),l); 

end 

%HARMONIC  3 


displ  freqlO  hann3((l+loop2),l)=fund_hannl_0((19+loop2*24),l)+i‘fund_hannl_0((20+loop2*24),l); 
dispffreqJO  haiTn3((l+loop2),l)=fund_hannl_0((4339+loop2*24),l)+i*fund_hannl_0((4340+loop2»24),l); 
dispffreqdO  harm3((l+loop2),l)=fiind_harml_0((8659+loop2*24),l)+i»fund_hamil_0((8660+Ioop2«24),l); 
displ  freqSO  harm3((l+loop2),l)=fund_haiml_0((12979+loop2»24),l)+i*fiind_harml_0((12980+loop2*24),l); 
disprfreql60_hann3((l+loop2),l)=fund_harml_0((17299+loop2*24),l)+i»fund_hannl_0((17300+loop2  24),1); 

end 

VoHARMONIC  4 

for  loop2=0:179 

displ  freqlO  hann4((l+loop2),l>=fund_hannl_0((21+loop2»24),l)+i*fund_harml_0((22+loop2  24),1); 
displ  freq20  hann4((l+loop2),l)=fund_hannl_0((4341+loop2»24),l)+i»fund_harml_0((4342+loop2»24),l); 
displ  freqdO  harm4((l+loop2),l)=fund_haiml_0((8661+Ioop2*24),l)+i*fiind_harml_0((8662+loop2*24),l); 
displ  freqSO  hann4((l+loop2),l)=fund  hannl_0((12981+loop2*24),l)+i»fund_harml_0((12982+loop2*24),I); 

disprfreql60_hann4((l+loop2),l)=fund_harml_0((17301+loop2»24),l)+i*fund_hannl_0((17302+loop2»24),l); 

end 

“/oHARMONIC  5 


displ  freqlO  hann5((l+loop2),l)=fund_hannl_0((23+loop2»24).l)+i»fund_hannl_0{(24+loop2*24).l); 
displ  freq20  haim5((l+loop2),l)=fiind_hannl_0((4343+loop2*24),l)+i*fiind_harml_0((4344+loop2»24),l); 
displ  freq40  hann5((l+loop2),l)=fiind_hannl_0((8663+loop2»24),l)+i*fund_hannl_0((8664+loop2*24),l); 
displ  freqSO  hann5((l+loop2Xl)=&nd_hannl_0((12983+loop2»24),l)+i»fund_hannl_0((12984+loop2»24),l); 
displ_freql60_hann5((l+loop2),l)=fiind_hannl_0((17303+loop2»24),l)+i*fund_hannl_0((17304+loop2*24),l); 

end 

%  Write  out  frequencies  recorded 
freqlO(l :  180,l)=parametersl_0((3:6;  i080),l); 
freq20( 1:180,1 )=parameters  1  _0(( 1 0 83 : 6 : 2 1 60),  1 ); 
freq40(l :  180, 1  )=parametersl_0((2 163 :6:3240),  1); 
freq80(l:180,l)-parametersl_0((3243:6:4320),l); 
freql60(l :  1 80,  l)^)arametersl_0((4323:6:5400),  1); 
elseif  amplitude^- 2.0’ 
load  parameters2_0 1 ; 
load  fund_hann2_0 1 ; 
loadnoise2_01; 
load  accel_A2_01; 
loadaccel_B2_01; 
load  accel_noise2_0 1 ; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
%FUNDAMENTAL 


displ_freqlO_fund((l+Ioop2),l)=lund_harm2  0((13+loop2*24),l)+i*fund_harm2_0((14+loop2*24),l); 
displ  freq20  fund((l+loop2).l)=fund_liann2_0((4333+loop2*24),l)+i*fiind_hann2_0((4334+loop2*24),l); 
displ  fieq40  fUnd((l+loop2),l)=fund_hann2_0((8653+loop2»24),l)+i*fund_harm2_0((8654+loop2*24),l); 
displ  freqSO  fund((l+loop2),l)=fiind_hann2_0((12973+loop2“24),l)+i*fund_harm2_0((12974+loop2“24),l); 
displ_freql60_fund((l+loop2),l)=fund_hanii2_0((17293+loop2»24),l)+i»fund_harm2_0((17294+loop2»24),l); 

end 

%HARMONIC  1 


for  loop2=0: 179 

displ  freqlO  harml((l+loop2),l)=fund_hann2_0((15+loop2*24),l)+i*fund_hann2_0((16+loop2»24),l); 
displ  freq20_hannl((l+loop2).l)=fund_haiin2_0((4335+loop2*24).l)+i*fund_hann2_0((4336+loop2“24),l); 

displ  freq40  harml((l+loop2),l)=fund  harm2  0((8655+loop2»24),l)+i»fund_hann2_0((8656+loop2*24),l); 

displ  freqSO  hannl((l+loop2),l)=fundlhann2_0((12975+loop2*24),l)+i*fund_hann2_0((12976+loop2*24),l); 
displ_fieql60_hannl((l+loop2),l)=fund_harm2_0((17295+loop2*24),l)+i*fund_harm2_0((17296+loop2*24),l); 
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end 

%HARMONIC  2 
for  Ioop2=0:179 

displ  freql0_hann2((l+loop2),l)=fund_hanii2_0((17+loop2»24),l)+i*fund_harm2_0((18+loop2»24),l); 

displ_freq20_hartn2((l+loop2),lHund_hann2_0((4337+loop2»24),l)+i*fund_harm2_0((4338+loop2*24).l); 

displ  fi«q40_hanii2((l+loop2il)=fund_hann2_0((8657+loop2*24),lHi*fund_hann2_0((8658+loop2»24),l); 
displ  freq80_hann2((l+loop2il)=fund_hanii2_0((12977+loop2*24),l)+i*fiind_hann2_0((12978+loop2*24),l); 

displ_freql60_hann2((l+loop2),l)=fund_hann2_0((17297+loop2»24),l)+i*fl>nd_hann2_0((17298+loop2*24),l); 

end 

%HARMONIC  3 
for  loop2=0:179 

displ_freql0_hann3((l+loop2),l)=lund_hann2_0((19+loop2*24),l)+i*fund_hann2_0((20+loop2*24),l); 

displ_freq20_hami3((l+loop2il)=fund_hami2_0((4339+loop2*24Xl)+i*fund_liarm2_0((4340+loop2»24),l); 

displ_freq40_hann3((l+loop2),lHund_lianii2_0((8659+loop2»24),l)+i*iund_harm2_0((8660+loop2*24),l); 

displ_fi«q80_hann3((l+loop2),lHund_harm2_0((12979+loop2*24),l)+i*fund_liarm2_0((12980+loop2»24),l); 

displ_freql60_harm3((l+loop2),l)=fund_hann2_0((17299+loop2»24).l)+i»fiind_hann2_0((17300+loop2*24),l); 

end 

%HARMONIC  4 
forloop2=0:179 

displ_freql0_hann4((l+loop2),l)=fund_hanii2_0((21+loop2»24),l)+i*fund_hann2_0((22+loop2*24),l); 

displ_fteq20_hann4{(l+loop2Xl^nd_hann2_0((4341+loop2»24),l)+i*fund_harm2_0((4342+loop2*24),l); 

displ_freq40_hami4((l+loop2Xl^nd_hann2_0((8661+loop2*24),l)+i*fund_hann2_0((8662+loop2»24),l); 

displ_freq80_hann4((l+loop2),l)=fijnd_hann2_0((12981+loop2*24),l)+i*fund_haiin2_0((12982+loop2*24),l); 

displ_freql60_hann4((l+loop2),l)=fund_hann2_0((17301+loop2*24),l)+i*fund_hann2_0((17302+loop2»24),l); 

end 

VoHARMONIC  5 
for  loop2=0:179 

displ_freql0_hann5((l+loop2),l)=fund_hann2_0((23+loop2*24),l)+i*fUnd_hann2_0((24+loop2*24),l); 

displ_freq20_harm5((l+loop2),l)=fund_hanii2_0((4343+loop2»24),l)+i*fund_hann2_0((4344+loop2*24),l); 

displ_freq40_hann5((l+loop2),l^nd_hann2_0((8663+loop2»24),l)+i*fund_harm2_0((8664+loop2»24),l); 

displ  freq80_hann5((l+loop2y^nd_harm2  0((12983+loop2*24),l)+i*fund_hann2_0((12984+loop2*24),l); 

displ_freql60_hann5((l+loop2),l)=fund_hann2_0((17303+loop2*24),l)+i»fund_hann2_0((17304+loop2*24),l); 

end 

%  Write  out  frequencies  recorded 
freql0(l:180,l)=parameters2_0((3:6:1080),l); 
freq20( 1:180,1 )=parameters2_0(( 1083:6:21 60),  1 ); 
freq40(l:180,l)=parameters2_0((2163:6:3240),l); 
freq80(l:180,l)=parameters2_0((3243:6:4320),l); 
freql60(l:180,l)=parameters2__0((4323:6:5400),l); 
elseif  amplitude=='4.0' 
load  parameters4_0 1 ; 
load  fund_harm4_0 1; 
load  noise4  01; 
load  accel_A4_0 1 ; 
load  accel_B4_0 1 ; 
load  accel_noise4_0 1; 

%  Combine  real  and  imaginary  parts  of  displ  freq 

%FUNDAMENTAL 

for  loop2=0:179 

displ_freql0_fund((l+loop2),l)=fund_harm4_0((13+loop2*24),l)+i*fund_harm4_0((14+loop2*24),l); 

displ_fieq20_£und((l+loop2^1)=fund_hann4_0((4333+loop2*24y)+i*fund_bann4_0((4334+loop2»24),l); 

displ_fleq40_fund((l+loop2),l)=fund_hann4_0((8653+loop2»24),l>i-i*£und_harm4_0((8654+loop2*24),l); 

displ_freq80_fimd((l+loop2).l)=fund_hann4_0((12973+loop2»24),l)+i*fiind_harm4_0((12974+loop2*24),l); 

displ_freql60_fund((l+loop2),l)=fund_hann4_0((17293+loop2»24),l)+i*fund_harm4_0((17294+loop2*24),l); 

end 

%HARMONIC  1 
forloop2=0:179 

displ_freql0_hannl((l+loop2),l)=fUnd_ham4_0((15+loop2»24),l)+i*fund_hann4_0((16+loop2»24),l); 

displ_freq20_hannl((l+loop2Xl)=fiind_harm4_0((4335+loop2»24),l)+i*fund_hann4_0((4336+loop2»24),l); 

displ_freq40_haiml((l+loop2yHund_hann4_0((8655+loop2*24),l)+i»fund_hann4_0((8656+loop2*24),l); 

displ_freq80_hannl((l+loop2y)=ftind_hann4_0((12975+loop2»24),l)+i»fund_hann4_0((12976+loop2»24),l); 

displ_freql60_harml((l+loop2),l)=fund_hami4_0((17295+loop2*24),l)+i*fi>nd_hann4_0((17296+loop2*24),l); 

end 

%HARMONIC  2 
forloop2-0:179 

displ_freql0_harm2((l+loop2),l)=fUnd_harm4_0((17+loop2*24),l)+i*fund_hann4_0((18+loop2*24),l); 
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displ_freq20_hann2((l+loop2),l)=fund_hann4_0((4337+loop2»24),l)+i*fund_harm4_0((4338+loop2*24),l); 

displ_freq40_hami2((l +loop2),  l^nd_hann4  0((8657+loop2  *24),  1  )+i*fund_hann4_0((865 8+loop2*24),  1 ); 
displ_freq80_hann2((l+loop2),l^nd_harm4_0((12977+loop2*24),l)+i»fund_hann4_0((12978+loop2*24),l); 
displ_freql60_haiTO2((l+loop2),l)=fiind_hann4_0((17297+loop2*24),l)+i»fund_hann4_0((17298+Ioop2*24),l); 
end 

%HARMONIC  3 
forloop2=0:179 

displ_freql0_hanTi3((l+loop2),l)=fund_hann4_0((19+loop2*24),l)+i*fond_harm4_0((20+loop2*24),l); 

displ  jfreq20_harm3((l+loop2),l)=fUnd_harm4_0((4339+loop2*24),l)+i*fUnd_harm4_0({4340+loop2*24),l); 

disprfreq40_^harm3((l+loop2y)=fund_harm4_0((8659+loop2*24),l)+i*fund_harm4_0((8660+loop2*24),l); 

displ_freq80_hann3((l+loop2),l)=fund_harm4_0((12979+loop2*24),l)+i*fund_hann4_0((12980+Ioop2*24),l); 

displ_freql60_harm3((l+loop2),l)=fundjiarm4_0((17299+loop2*24),l)+i*fund_harm4_0((17300+loop2*24),l); 

end 

%HARMONIC  4 
forloop2N):179 

displ_freql0_harm4((l+loop2),l)=fund_harm4_0((21+loop2*24),l)+i*fund_hann4_0((22+loop2*24),l); 

displ_freq20_harm4((l+Ioop2),l)=fund_harm4_0((4341+loop2*24),l)+i*fiind_hann4_0((4342+Ioop2*24),l); 

displ_freq40_harm4((l+loop2),l)=fiind_hann4_0((8661+loop2*24),l)+i*fiind_hann4_0((8662+loop2*24),l); 

displ_freq80_harm4((l+loop2),l)=aind_hami4_0((12981+loop2*24),l)+i*fund_harm4_0((12982+loop2*24),l); 

displ_freql60_harm4<(l+loop2),lHiind_hann4_0((17301+loop2*24),l)+i*fund_harm4_0((17302+loop2*24),l); 

end 

%HARMONIC  5 
for  loop2=0:179 

dispI_freql0_hami5((H-ioop2),l)=fund_hann4_0((23+loop2*24),l)+i»fund_hann4_0((24+loop2*24),l); 

displ  freq20_hami5((H-loop2),lHund_harm4_0((4343+loop2*24),l)+i*fund_harm4_0((4344+loop2*24),l); 

disprfreq40_hanii5((l+loop2Xl)=fund_harm4_0((8663+loop2*24),l)+i*fund_harm4_0((8664+loop2*24),l); 

displ_freq80_hann5((l+loop2),l)=fijnd_hami4_0((12983+loop2*24),l)+i*fund_harm4_0((12984+loop2*24),l); 

displ_freql60Jiarm5((l+loop2),l)=fiind_hann4_0((17303+Ioop2*24),l)+i*fund_harm4_0((17304+loop2»24),l); 

end 

%  Write  out  frequencies  recorded 
freq  1 0(  1 : 1 80, 1  )=paramelers4_0((3 : 6 : 1 0 8 0),  1 ); 
freq20(l :  180,l)=parameters4_0((1083:6:2 160),1); 
freq40( 1:180,1 )=parameters4_0((2 1 63 :6 :3 240),  1 ); 
freq80(  1:180,1 )=parameters4_0((3 243 :6 :43 20),  1 ); 
freql 60( 1 : 1 80,  l)=parameters4_0((4323 :6 :5400),  1); 
elseif  amplitude=- 8.0' 
load  parametersSO  1 ; 
load  fundharm8_0 1; 
load  noise8_01; 
load  accel_A8_01; 
load  accel_B8_01; 
load  accel_noise8_01; 

%  Combine  real  and  imaginary  parts  of  displ  freq 

%FUNDAMENTAL 

for  Ioop2=0:l79 

displ_freql0_flind((l+loop2),l>=fund_hann8_0((13+loop2*24),l)+i*fund_hann8_0((14+loop2*24),l); 

displ_freq20_fiind((  1 +loop2),  1  )=fund_harm8_0((433 3 +loop2  *24),  1  )+i  *fund_hajTn8_0((43 34+loop2  *24),  1 ); 

displ_freq40Jund((l+loop2y)=fund_harm8_0((8653+loop2*24),l)+i*fund__harm8_^0((8654+loop2*24),l); 

displ_freq80_fund((l+loop2y)=fiind_hann8^0((12973+loop2*24),l)+i*fund_harm8_0((12974+Ioop2*24),l); 

displ_freql60_fun^(l+loop2),l)=fund_harm8_0((17293+loop2*24),l)+i*fund_harm8_0((17294+loop2*24),l); 

end 

%HARMONIC  1 
forloop2=0:179 

displ_freql0_harml((l+loop2).l)=fund_harm8_0((15+loop2*24),l)+i*fund^harm8^0((16+Ioop2*24),l); 

displ_freq20  harml((l+loop2),l)=fund_hann8_0((4335+loop2*24Xl)+i*fiind_harm8_0((4336+loop2*24),l); 

displ_freq40_harml((l+loop2),l)=fiind_harm8_0((8655+loop2*24),l)+i*fund_harm8_0((8656+loop2*24),l); 

displ_freq80_hannl((l+loop2y)=fund_harm8_0((12975+loop2*24),l)+i*fiind_harm8_0((12976+loop2*24),l); 

displ_freql60_haiml((l+loop2),l)=fund_hann8^0{(17295+loop2*24),l)+i*fund_harm8_0((17296+loop2*24),l); 

end 

%HARMONIC2 

forloop2=0:179 

displ  freql0_hann2((l+loop2),l)=fund_hann8_0((17+loop2*24),l)+i*fund_hann8_0((18+loop2*24),l); 

disprfreq20_hann2((l+loop2Xl)=fund_hann8_0((4337+loop2«24),l)+i*fund_hann8_0((4338+loop2*24),l); 

displ_freq40_hann2({l+loop2y)=&nd_hami8_0((8657+loop2*24il)+i*fund_harm8_0((8658+Ioop2*24),l); 

displ_freq80_hann2((l+loop2),l>=fund_hann8_0((12977+loop2*24),l)+i*fund_hann8_0((12978+loop2»24),l); 

displ_freql60_hann2((l+loop2),l)=fund_harm8_0((17297+loop2*24),l)+i»fund_hann8_0{(17298+loop2*24),l); 
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end 

%HARMONIC  3 


dispi  freqio  hann3((l+loop2),lHund_hann8J((19+loop2*24)4)+i‘fUnd_haTO8_0((20^^^^  HL  n- 

dispffreq20"hanti3((l+Ioop2),l)=fUnd_harm8_0((4339+loop2»24),l)+i*fund_hanT,8_0((4340+  oop2  24)  ). 

disprfieq40"hanii3((l+Ioop2),l)=fund_hann8_0((8659+loop2»24Xl)+i»fund_harm8_0((8660+loop2  24)1); 
yrfre380~hann3((l+loop2yHi'nd_hann8_0((12979+loop2*24)J)+i*fund  h^8  0((12980+loop2^^ 

disprfreql60_hanii3((l+loop2),l)=fund_hann8_0((17299+loop2‘24).l)+i*ft*nd_hann8_0((17300+loop2  24).!). 

end 

%HARMONIC  4 

SrKHj^h^4((l+loop2Xl)=fUnd_ham8_0((21+loop2»24).l)+i*fUnd_ham8_0((22+l^^^^^ 
disprfreq20'hann4((l+IcK>p2)a)=fund_haim8J((4341+loop2»24)4)+i*fund_^8_0((«42+bop2m 
dispffi440"hann4((l+loop2),l)=^nd_hann8_0((8661+loop2*24),l)+i*fUnd_harm8_0((8662+I^»<V2  24)1) 

disjrfre380"hann?(l+loop2)4)=*nd_harm8_0((12981+loop2»24),l)+i»&nd  h^8  0((12982+1^2  24y 

disprfreql60_han^(l+loop2),l)=fund_harm8_0((17301+loop2*24),l)+i‘fund_harm8_0({17302+loop2  24),1), 

end 

“/oHARMONIC  5 


dispi  freqio  hann5((l+loop2)4)=fUnd_hann8_0((23+loop2*24)4)+i*fund_hann8_0((24+to^  24)  1); 

disprfreq20  hann5((l+loop2)4)=fUnd_hann8_0((4343+loop2n4)4)+i*&nd_ham8_0  4344+  oop2  24)^,  )^ 

di^rfreq40"harm5((l+loop2),l)=*nd_harm8_0((8663+loop2*24),l)+i*fi.nd_h^8  0((8664+loop2  24)1) 
dii“freq80“hann5((l+loop2),l)=li.nd_hann8_0((12983+Ioop2»24),l)+i*fund  h^8  0((12984+1^2  24)  1) 
dispffreqieo  hann5((l+loop2),l)=fund_hami8_0((17303+loop2»24),l)+i*fimd_hann8_0((17304+loop2  24),1), 


end 

%  Write  out  frequencies  recorded 

freql  0(  1:180,1  )=paranieters8_0((3 :6 : 1 080),  1 ); 

freq20(l :  180,  l)=parameters8_0((1083  ;6:2 1 60),  1); 

freq40( 1 : 1 80, 1 )=parameters8_0((2 1 63 :6 :3240),  1 ); 

freq80(l :  1 80,  l)=parameters8_0((3243:6:4320),  1); 

freql60(l:180,l)=parameters8_0((4323:6:5400),l); 


end; 


elseif  iteration^- 2* 

if  existCamplitude')==0 

aniplitude=inputC Amplitude  (enter  3  char):  ',V); 
end; 

if  anq)litude='0.5' 

load  parameters0_52; 
load  fund_harm0_52; 
load  noise0_52; 
load  accel_A0_52; 
load  acceI_B0_52; 
load  accel_noise0_52; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
%FUNDAMENTAL 


dispi  freqio  fund((l+loop2)4)=*nd>ami0J((13+loop2*24)4)+i*fund>aim0J((14+l^^^^^  24),1); 
dispffreq20"fund((l+loop2)4)=*nd_harm0_5((4333+loop2*24)4)+i*fond_ham0_5  4334 joop2  24  4^ 
di™rfreq40"fiind((l+loop2),l)=fUnd_hannO_5((8653+loop2»24),l)+i*&nd_hannO_5((8654+loop2  24)1), 
disprfreq80"fund((l+loop2Xl)=fund_hann0_5((12973+loop2*24),l)+i*fi.nd  h^O  5((12974+Wp2  24)1); 
disprfieql60_fund((l+loop2),l)=fund_hann0_5((17293+loop2*24),l)+i»fund_hann0_5((17294+loop2  24),1), 


end 

%HARMONlC  1 


dispi  LilO  harail((l+loop2),l)=*nd_ha™0J((15+loop2*24),l)+i*fund_ha™0_5(^^^^^^^  24).l); 

disprfreq20”haml((l+loop2),l)=*nd_hann0J{(4335+Ic)op2*24),l)+i*fund_harm0_5((4336+oop2  24).  ) 
disprfreq40"hannl((l+Ioop2),l)=fund_harm0J((8655+loop2*24),l)+i*fond  h^0  5((8656+1^^^ 

disprfreq80"hannl((l+Ioop2Xl)=fund_hann0J((12975+loop2‘24),l)+.*fund  h^O  5((12976+loop2  24)  1) 
d4rfi4l60_harail((l+loop2),l)=fund_hann0_5((17295+loop2*24),l)+i»fund_hann0_5((17296+loop2  24),1), 

end 

%HARMONIC2 

dUprfreq7o\ann2((l+loop2),l)=fund_hannOJ((17+loop2*24),l)+i*fund_hannOJ((18+^^^^^^^ 
dis5rfreq20"hann2((l+loop2),l)=fund_hann0_5((4337+Ioop2*24),l)+i»lund_harm0_5((4338+ bop2  24  .  ); 

disprfreq40"hann2((l+loop2ll)=fund_hann0_5((8657+loop2»24),l)+i»fund_haim0_5((8658+loop2*24),l); 


139 


dispi  freq80  hann2((l+loop2)4)=fiind_hann0J((12977+Ioop2‘24)4)+i*fund  h^O  5((12978+l^p2*24^^^^^ 
displ>eql60_hann2((l+loop2),lHund_harm0_5((17297+loop2*24),l)+i*fund_harm0_5((17298+loop2  24),1), 

end 

•/oHARMONIC  3 

Sl^5^((l+loop2)>l)=&"d_ham0J((19+loop2»24)4)+i»innd_l^^ 

dispffreq20  haitii3((l+Ioop2)4Hiind_haim0_5((4339+loop2»24)4)+i*&nd_hann0_5((«40+ oop2  24),1), 
dispffreqdO  hami3((l+loop2),l)=fund_harm0_5((8659+loop2*24).l)+i»fUnd_harm0_5((8660+loop2  24)  1); 
disprfreq80  hann3((l+Ioop2Xl)=fUnd_hann0_5((12979+loop2»24),l)+i*fiind_hann0_5((12980+loop2»24)  l); 
d5rfreql60_hann3((l+loop2),l)=fund_hann0_5((l7299+loop2»24),l)+i*fimd_hann0_5((17300+loop2  24), 1); 

end 

VoHARMONIC  4 

dkpl'TLH)lhLi4((l+loop2),l)=fund_haini0_5({21+loop2*24),l)+i*fund_harm0_5((22+loop2»24),2; 

disprfreq20_hann4((l+loop2),l)=fUnd_hann0_5((4341+loop2»24),l)+i*fund_hann0_5((4342+  oop2  24)  ); 

dispi  freq40  hann4<(l+loop2),l)=fiind_harm0_5((8661+loop2*24),l)+i*fund_hann0_5((8662+l^oop2  24), 1), 
disprfreq80:hann4{(l+loop2),l)=fund_haim0_5((12981+loop2»24),l)+i*fund_harm0_5((12982+loop2  24)1); 
disprfreql60_hann4((l+loop2),l)=fiind_hann0_5((17301+loop2»24),l)+i*fund_harm0_5((17302+loop2  24),1); 

end 

%HARMONlC  5 

dispi  freql0_hann5((l+loop2),l)=fund_hann0_5((23+loop2»24),l)+i*fund_hann0_5((24+l(wp2*24),l); 
disprfreq20_harm5((l+loop2Xl)=fiind_hann0_5((4343+loop2»24),l)+i*fund_harm0_5((4344+  oop2  24)  ); 

dispi  freq40  hann5((l+lcx)p2),l)=fund_hann0_5((8663+loop2»24),l)+i*fund_hann0_5((8664+lciop2  24),  1), 
di^ffrejso  hann5((l+loixi)=fUnd_ham0_5((12983+loop2»24),l)+i*fund_harm0_5((12984+l<wp2  24)1); 

dUprfreql60_hann5((l+loop2),l)=fund_liann0_5((17303+loop2*24),l)+i»fund_harm0_5((17304+loop2  24),1); 

end 

%  Write  out  frequencies  recorded 

freql  0(1 : 1 80,  l)=parameters0_5((3:6: 1 080),  1); 

freq20( 1 : 1 80,  l)=paraineters0_5((  1 083 :6:2 1 60),  1 ); 
freq40(l :  180,l)=paranieters0_5((2 163 :6:3240),  1); 
fi^q80( 1:180,1 )=parameters0_5((3 243 :6 :4320),  1 ); 
freql60(l:180,l)=parameters0_5((4323:6:5400),l); 
elseif  aniplitude=- 1.0’ 
load  parameters!  02; 
load  fund_harml_02; 
load  noisel_02; 
load  accel_Al_02; 
load  accel_Bl_02; 
load  accel_noisel_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 

%FUNDAMENTAL 

for  loop2=0:179 

dispi  freql0_fund((l+loop2),l)=fund_hannl_0((13+loop2*24),l)+i*fund_harml_0((14+loop2*24),l); 

dispi  freq20  fimd((l+loop2),l)=fund_hannl_0((4333+loop2*24),l)+i»fund_hannl_0((4334+loop2»24)  1); 
disprfreq40_£und((l+loop2),l)=fund  hannl_0((8653+loop2*24),l)+i*fiind_hannl_0((8654+loop2*24),l); 
disprfreq80  fund((l+loop2),l)=fUnd_hannl_0((12973+loop2»24),l)+i*fund_harml_0((12974+loop2  24),1); 
displ_freql60_fund((l+loop2),l)=fund_harml_0((17293+loop2»24),l)+i*fund_hannl_0((17294+loop2  24),1); 

end 

VoHARMONIC  1 

dbirfr^lO  harml((l+loop2),l)=fund_hannl_0((15+loop2*24),l)+i*fund_liamil_0((16+loop2*24)  l); 
dispi  freq20  hannl((l+loop2),l)=fund_hannl_0((4335+loop2«24),l)+i»fund_hamil_0((4336+ bop2  24),  ); 
disprfreq40_hannl((l+loop2),l)=fund_hannl_0((8655+loop2»24Xl)+i‘lund_harml_0((8656+loop2»24),l); 
dispi  freq80  hannl((l+loop2),l)=fund_hannl_0((12975+loop2*24),l)+i*fund_hannl_0((12976+loop2  24),1); 
disprfieql60_hannl((l+loop2),l)=fiind_harml_0((17295+loop2»24),l)+i*fund_harml_0((17296+Ioop2  24),1); 

end 

"/.HARMONIC  2 

dUrf°frMl0  hann2((l+loop2),l)=fiind_hannl_0{(17+loop2*24),l)+i*fund_hannl_0((18+loop2*24),l); 
dispi  freq20  hann2((l+loop2),l)=fund_hannl_0((4337+loop2»24),l)+i‘fund_harml_0((4338+loop2*24),l); 
dispi  freq40  hann2((l+loop2),l)=fund_hannl_0((8657+loop2*24),l)+i*fund_hannl_0((8658+loop2  24),1); 
dUpffreq80  hann2((l+loop2),l)=fund_hannl_0((12977+loop2*24),l)+i»fund_harml_0((12978+loop2  24),1); 
displ_freql60_hann2((l+loop2),l)=fUnd_hannl_0((17297+loop2*24),l)+i*fund_hannl_0((17298+loop2»24),l); 

end 

"/.HARMONIC  3 
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disprfreql0_hanti3((l+loop2),l)=fund_harml_0((19+loop2*24),l)+i*fUnd_hannl_0((20+loop2»24),l); 
displ  freq20  hann3((l+loop2),l)=fiind_hannl_0((4339+loop2*24),l)+i*fi<nd_harml_0((4340+loop2*24),l); 
disprfreq40  harm3((l+loop2),l)=fund_hannl_0((8659+loop2*24),l)+i*fiind_hannl_0((8660+loop2»24),l); 
displ  freqSO  harm3((l+Ioop2),l)=fund_hannl_0((12979+Ioop2»24),l)+i*fund_hannl_0((12980+loop2*24),  1); 
displ_freql60_hann3((l+loop2),l)=fund_hannl_0((17299+loop2*24),l)+i*fund_hannl_0((17300+loop2*24),l); 

end 

%HARMONIC  4 


for  loop2=0:179 

displ  freql0_hann4((l+loop2),lHund_hannl_0((21+loop2*24),l)+i*fund_hannl_0((22+loop2*24),l); 

displ  freq20_harm4((l  +loop2),  1  )=fund_harm  1_0((434 1 +loop2*24),  1  )+i*fund_harml _0((4342+loop2 *24),  I ); 
displ  fi«q40  hann4((l+loop2),lHund_haniil_0((8661+loop2*24),l)+i*fund_hannl_0((8662+loop2»24),l); 
displ  freqSO  hann4((l+loop2),l)=fijnd_hannl_0((12981+loop2»24),l)+i*fiind_hannl_0((12982+loop2*24),l); 
displ_fieql60_hann4((l+loop2).l)=fund_haiml_0((17301+loop2*24),l)+i*fund_hannl_0((17302+loop2*24),l); 

end 

%HARMONIC  5 
for  loop2=0:179 

displ  fteql0_hann5((l+loop2),lHund_harml_0((23+loop2*24).l)+i*fund_hannl_0((24+loop2  24),1); 
displ  freq20  hann5((l+loop2),l)=fund_hannl_0((4343+loop2*24),l)+i*fiind_harml_0((4344+loop2»24),l); 
displ  freq40  hann5((l+loop2),l)=fund_hannl_0((8663+loop2»24),l)+i*fiind_harml_0((8664+loop2*24),l); 
displ  freq80  hann5((l+loop2),l)=fund_hamil_0((12983+loop2»24),l)+i*fund_harml_0((12984+loop2*24),l); 
displ_freql60_hann5((l+loop2),l)=fund_liarml_0((17303+loop2»24),l)+i*fund_harml_0((17304+loop2*24),l); 

end 

%  Write  out  frequencies  recorded 
freql0(l:180,l)^)arametersl_0((3:6:1080),l); 
freq20( 1 : 1 80, 1 )=parameters  1_0(( 1 0 83 :6 :2 1 60),  1 ); 

freq40(l :  1 80,  l)=parametersl_0((2 163:6:3240),  1); 
freq80(l :  1 80,  l)=parametersl_0((3243 :6:4320),  1); 
freql60(l:180,l)==parainetersl_0((4323:6:5400),l); 
elseif  amplitude=='2.0’ 
load  parameters2_02; 
load  fund_hann2_02; 
load  noise2_02', 
load  accel_A2_02; 
load  acceI_B2_02; 
load  accel_noise2_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
%FUNDAMENTAL 

for  loop2=0: 179  ^ 

displ  freql0_fund((l+loop2),l)=lund_hann2_0((13+loop2*24),l)+i*fiind_harm2_0((14+loop2*24),l); 
displ  freq20_fimd((l+loop2),l)=fiind_hanii2_0((4333+loop2*24),l)+i»fund_hann2_0((4334+loop2*24),l); 
displ  freq40  £und((l+loop2Xl^nd_haiTn2_0((8653+loop2*24),l)+i*fund_harm2_0((8654+loop2*24),  1); 
displ  fieqSO  fund((l+loop2),l)=fund_hann2_0((12973+loop2*24),l)+i*fUnd_hann2_0((12974+loop2*24),l); 
displ_freql60_fund((l+loop2),l)=fund_harm2_0((17293+loop2»24),l)+i*fiind_harm2_0((17294+loop2*24),l); 

end 

%HARMONIC  1 
for  loop2=0:179 

displ_freql0_harml((l+loop2),l)=fund_harm2^0((15+loop2*24),l)+i*fund_hann2_0((16+loop2*24),l); 

displ  freq20  hannl((l+loop2),l)=fund_hann2_0((4335+loop2*24),l)+i*fiind_hann2_0((4336+loop2»24),l); 
displ  freq40_hannl((l+loop2),l)=fund_hann2_0((8655+loop2*24),l)+i*ftind_hann2_0((8656+loop2*24),l); 
displ  freq80  hannl((l+loop2),l)=fund_harm2_0((12975+loop2*24),l)+i*fund_harm2_0((12976+loop2»24),l); 
displ_freql60_hannl((l+loop2),l)=fund_harm2_0((17295+loop2*24),l)+i*fund_hann2_0((17296+loop2*24),l); 

end 

%HARMONlC2 


displ_freql0_harm2((l+loop2),l)=fund_harm2_0((17+loop2*24),l)+i*fund_harm2_0((18+loop2*24),l); 

displ_freq20_hann2((l+loop2),l)=fund_hann2_0((4337+loop2»24),l)+i*fiind_harm2_0((4338+loop2*24),l); 

displ_freq40_hami2((l+loop2),l)=fund_hami2_0((8657+loop2»24),l)+i*fund_hann2_0((8658+loop2*24),l); 

displ  freq80_hann2((l+loop2),l)=fund_hann2_0((12977+loop2»24),l)+i*fund_hann2_0((12978+loop2»24),l); 
displ_freql60_harm2((l+loop2),l)=fund_hann2_0((17297+loop2*24),l)+i*lund_hann2_0((17298+loop2*24),l); 

end 

%HARM0N1C3 

for  loop2=0:179 

displ  freql0_hanii3((l+loop2),l)=fund_hann2_0((19+loop2»24),l)+i*fUnd_haiin2_0((20+loop2»24),l); 

displ_freq20_hann3((l+loop2),l)=fund_hann2_0((4339+loop2»24),l)+i*fund_harm2_0((4340+loop2»24),l); 

displ_freq40_hann3((l+loop2),l)=fund_hann2_0((8659+loop2»24),l)+i*lund_harm2_0((8660+loop2*24),l); 
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displ  freqSO  hann3((l+loop2),lHund_hann2_0((12979+loop2»24),l)+i*fiind_hann2_0((12980+loop2‘24).l); 
displ_freql60_hann3((l+loop2),l)=fund_hami2_0((17299+loop2*24),l)+i*fi«nd_hann2_0((17300+loop2*24),l); 

end 

%HARMONIC  4 

for  loop2=0: 1 79  ,  ^ 

displ  freql0_harm4((l+Ioop2),l)=fund_hann2_0((21+loop2»24),l)+i*fund_harm2_0((22+loop2»24),l); 

displ  fieq20  hann4{(l+loop2),l)=lund_harm2_0((4341+loop2*24),l)+i*fiind_harm2_0((4342+loop2*24),l); 
disprfieq40_harm4((l+loop2),l)=fund_hanii2_0((8661+loop2»24),l)+i*fund_harm2_0((8662+loop2‘24),l); 
displ  freqSO  hann4((l+loop2),l)=fund_haitn2_0((12981+loop2*24),l)+i*fund_harm2_0((12982+loop2*24),l); 
displ_freql60_harm4((l+loop2),l)=fund_haim2_0((17301+loop2*24),l)+i*fund_hann2_0((17302+loop2*24),l); 
end 

%HARMONIC  5 
forloop2=0:179 

displ  freql0_hann5((l+loop2),l)=fund_ham)2_0((23+loop2»24),l)+i»fund_harm2_0((24+loop2*24),l); 
displ  fieq20  hann5((l+loop2),l)=fund_hami2_0((4343+loop2*24),l)+i*fi*nd_hann2_0((4344+loop2*24),l); 
displ  freq40  ha™5((l+loop2),l)=fund_hann2_0((8663+loop2»24),l)+i»fund_hann2_0((8664+loop2»24),l); 
displ  freq80_hann5((l+loop2),l)=fund_hann2_0((12983+loop2»24),l)+i»fund_hann2_0((12984+loop2*24),l); 
displ_freql60_hann5((l+loop2),l)=fund_harm2_0((17303+loop2*24),l)+i*fund_hann2_0((17304+loop2»24),l); 

end 

%  Write  out  frequencies  recorded 

freqlO(l :  180,l)=parameters2_0((3:6: 1080),  1); 

freq20( 1 : 1 80, 1 )=parameters2_0((10 83 :6:2 1 60),  1 ); 

freq40( 1:180,1 )=parameters2_0((2 1 63 :6 :3240),  1 ); 
freq80(l :  180,  l)=parameters2_0((3243 :6:4320),  1); 
freql60(l  :180,  l)^arameters2_0((4323:6:5400),  1 ); 
elseif  ainplitude=='4.0’ 
load  parameters4_02; 
load  fund_harm4_02; 
load  noise4_02; 
load  accel_A4_02; 
load  acce!_B4_02; 
load  accel_noise4_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
%FUNDAMENTAL 

for  loop2=0:179  * 

displ  freql0_fund((l+loop2),l)=fund_harm4_0((13+loop2*24),l)+i*fimd_harm4_0((14+loop2*24),l); 
displ  freq20_fund((l+loop2),l)=lund_hann4_0((4333+loop2»24),l)+i*fund_hann4_0((4334+loop2*24).l); 
displ_freq40_fund{(l+loop2),l)=fiind_harm4_0((8653+loop2*24),l)+i*fund_harm4_0((8654+loop2*24),l); 
displ  freq80  fund((l+loop2).l)=fiind_hann4_0((12973+loop2*24),l)+i*fund_harm4_0((12974+loop2»24),l); 
displ_freql60_ftind((l+loop2),l)=fund_harm4_0((17293+loop2»24),l)+i*fund_hann4_0((17294+loop2*24).l); 

end 

%HARMONIC  1 
forloop2=0:179 

displ  freql0_hannl((l+loop2),l)=fiind_hann4_0((15+loop2»24),l)+i*&nd_harm4_0((16+loop2*24),l); 
displ  freq20_hannl((l+loop2),l)=fund_hann4_0((4335+loop2»24),l)+i*fiind_hann4_0((4336+loop2*24),l); 
displ  freq40_haniil((l+loop2),l)=fund_hann4_0((8655+loop2*24),l)+i*fi*nd_hami4_0((8656+loop2*24),l); 

disprfreq80_hannl((l+loop2),l)=fund_harm4_0((12975+loop2*24),l)+i*fiind_hann4_0((12976+loop2»24),l); 

displ_freql60_hannl((l+loop2),l)=flind_hann4_0((17295+loop2*24),l)+i*fund_hann4_0((17296+loop2*24),l); 

end 

%HARMONIC  2 
forloop2-0:179 

displ  freql0hann2((l+loop2),l)=fund_hann4_0((17+loop2»24),l)+i*fund_hann4_0((18+loop2*24),l); 

displ  freq20  hann2((l+loop2),l)=fund_hann4_0((4337+loop2»24),l)+i*fi>nd_harm4_0((4338+loop2»24),l); 

displ_freq40_hanti2((l+loop2),l)=fund_hann4_0((8657+loop2*24),l)+i*fund_hami4_0((8658+loop2*24),l); 

disprfreq80_hann2((l+loop2),lHund_hann4_0((12977+loop2*24),l)+i»fund_harm4_0((12978+loop2»24),l); 

displ_freql60_hami2((l+loop2),l)=fund_hami4_0((17297+loop2*24),l)+i»lund_hann4_0((17298+loop2»24),l); 

end 

%HARMONIC3 
for  loop2=0:179 

displ_fi«ql0_haiTn3((l+loop2),l)=fund_hann4_0((19+loop2*24),l)+i*fund_hann4_0((20+loop2*24),l); 
displ  freq20_hann3((l+loop2),l)=fund_hann4_0((4339+loop2*24),l)+i*fund_hann4_0{{4340+loop2*24),l); 

displ_freq40_hanii3((l+loop2),l)=fiind  hami4_0((8659+loop2*24),l)+i*fund_hann4_0((8660+loop2*24),l); 

disprfreq80_hann3((l+loop2),l)=fiind~harm4_0((12979+loop2*24),l)+i»fiind_harm4_0((12980+loop2»24),l); 

displ_freql60_hann3((l+loop2),l)=fiind_harm4_0((17299+loop2»24),l)+i*fimd_hann4_0((17300+loop2»24),l); 

end 

%HARMONIC  4 
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for  loop2  0.179  A  Viormd  +1ood2*24^  harm4_0((22+loop2*24),l), 

displ_freql0_hann4((l+bop2).  ^nd_«_0((21  1^^  24)  ) 

dispLfreq20_harm4((l+loop2),l^nd_to4_0((43  p2^2  )  ^  0((8662+loop2»24),l); 

dispLfreq40_hann4((l+  oop2 ,  ^nd_^4_0  8^  )  -  0((i2982+loop2»24),l); 

end 

•/oHARMONIC  5 

displ_fieq20_hann5((l+loop2Xl^nd_^  ((4343  p  ^  hann4_0((8664+loop2»24),l); 

displ_fi«l40_ham5((l+  oop2X  ^nd_ham4_0  8663+U^2  2^)  ; 

end 

%  Write  out  frequencies  recorded 
freql0(l :  1 80,  lHJaraineters4_0((3:6: 1 080),  1); 
freq20(l :  1 80,l)=parameters4_0((1083:6:2 160),  1); 
freq40(l :  1 80,  l)=parameters4_0((2 163 :6:3240),  1); 
freq80( 1 : 1 80, 1 )=parameters4_0((3243 :6 :43 20),  1 ); 

freql60(l;180,l)=parameters4_0((4323:6:5400),l); 
elseif  ainplitude=- 8.0' 
load  parameters8_02; 
load  fund_hann8_02; 
load  noise8__02; 
load  accel_A8_02; 
load  accel_B8_02; 
load  accel_noise8_02; 

%  Combine  real  and  imaginary  parts  of  displ_freq 
%FUND  AMENT  AL 

dUpLfreq20>nd((l+loop2),l^nd_h^8_0  (4333  oog^24)4^)^^ 

displ_fireq40>nd((l+  oop2y^nd_h^8_0  (8653  ^p2^) 

end 

%HARMONIC  1 

end 

%HARMONlC2 

displ_freq20_harm2((l+oop2)  ^nd_^8_0((433yo  p 

displ_freq40_harm2((l+ oop2),l^nd_h^  ((  657  1^  p  ^  0((12978+loop2»24),l); 

end 

%HARMONIC  3 

Sg;S:Ln3((K.oop2,l^^nd_^^^^^^^ 

dispLfreq20_hann3  +  oop2  . 

displ_freq40_harm3((l+  oop2X  ^nd_^  ((8659  P  P^^j2ggQ^,^p2»24^^ 

end 

“/oHARMONIC  4 
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displ  freqSO  hann4((l+loop2)4Hund_hami8_0({12981+loop2*24)J)+i*lund_harm8_0((12982+loop2*24)  lX 

dispffieql60_hami4((l+loop2),l)=fund_hajm8_0((17301+loop2*24).l)+i*fund_hann8_0{(17302+loop2  24).l); 

end 

•/oHARMONIC  5 

displ  ^10_harm5((l+loop2),l)=fund  harm8_0((23+loop2*24),l)+i*fund_harm8_0((24+llwp2*24Xl); 

disprfieq20  hann5((l+loop2).l)=fund_hann8_0((4343+Ioop2»24),l)+i*fund_harm8_0((4344+ oop2  24)  ); 

disprfreq40  harm5((l+loop2),l)=fund_hann8_0((8663+loop2*24),l)+i»fund_hann8_0((8664+loop2  24)1); 
dispffreqSO  hann5((l+loop2),l)=fund_hann8J((12983+loop2»24).l)+i»fund_ham8_0((12984+loop2  24)1X 

dispr^l60_hann5((l+loop2),l)=fund_hann8_0((17303+loop2*24),l)+i‘fiind_hann8_0((17304+loop2  24),1); 

end 

%  Write  out  frequencies  recorded 

freql  0(  1 : 1 80,  l)=parameters8_0((3 :6 : 1 080), 1 ); 

freq20(l :  1 80,l)=parameters8_0((1083 :6:2 160),  1); 

freq40{l :  1 80, 1  )=paranieters8_0((2 1 63 :6:3240),  1 ); 

freq80(l :  1 80,  l)=parameters8_0((3243 :6:4320),  1); 

freql60(l:180,l)=parameters8_0((4323:6:5400),l); 


%  Plot  results 

if  existCpositionT^^O 

position=input('Position  in  cm  (enter  5  char):  ’,'s'); 

end; 

if  position=- 10,00' 
figure(l) 

semilogy(freql  0,abs(displ_freq  1 0_fund),'-') 
hold  on 

semilogy(freql  0,abs(displ_freql  0_harm  1 ),'—') 

semilogy(freq  1 0,abs(displ_freq  1 0_hann2),'-.') 

semilogy(freq  1 0,abs(displ_freq  1 0_harm3  ),';') 

semilogy(freq  1 0,abs(displ_freq  1 0_harm4),'-') 

semilogy(freql  0,abs(displ_freq  1 0_harm5),'--’) 

title l='DispIacement  vs  Frequency*; 

title2='  Amplitude  = 

title3-  Position  =  '; 

title2=strcat(title2,amplitude,'Volts'); 

title3=strcat(title3,position,'cm'); 

title4='  Iteration  =  *; 

title4=strcat(title4,iteration); 

title^data=char({titlel,title2,title3,title4}); 

title(title_data); 

ylabel(Displacement') 

xlabel(Trequency  (Hz)') 

legend(Tund','Harm  1  VHann  2','Harm  3’, 'Harm  4','Harm  5*) 
hold  off 
orient  landscape 
figure(2) 

plot(freql  0,abs(displ_freql  0  Jund),'-') 
hold  on 

plot(freqlO,abs(displ_freql  Oharm  I),'-') 
plot(freq  1 0  ,abs(displ_freq  1 0_harm2),'-.') 
plot(freqlO,abs(displ_freq  1 0_harm3  ),':’) 
plol(freq  1 0,abs(displ_freq  1 0harm4),'-') 
plot(freql  0,abs(displ_freql  OharmS),'-') 
title  l=T)isplacement  vs  Frequency*; 
title2-  Amplitude  =  '; 
title3='  Position  =  ’; 
title2=strcat(title2,ampIitude,'Volts'); 
title3 =strcat(title3  ,position,'cm'); 
title4-  Iteration  =  ’; 
title4==strcat(title4, iteration); 
title_data=char((titlel,title2,title3,title4}); 
title(title_data); 
ylabelCDisplacement’) 
xlabelCFrequency  (Hz)*) 
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legendCFund’/Harm  r,Tiarm  2VHann  3','Harm  4VHann  5’) 
hold  off 

orient  landscape 
elseif  position— 20.00' 
figure(l) 

semilogy(freq20,abs(displ_freq20_fund),'-*) 

hold  on 

semilogy(freq20  ,abs(displ_fireq20_hann  1 

seniilogy(freq20,abs(displ_freq20_hann2),'-.’) 

semilogy(freq20,abs(displ_freq20_harm3),':') 

semilogy(freq20,abs(displ_freq20_haiTn4),'-') 

semilogy(freq20,abs(disp!_freq20_hann5),'— ') 

titlel- Displacement  vs  Frequency*; 

title2-  Amplitude-’; 

titles-  Position  =  '; 

title2=strcat(title2,ampIitude,'Volts’); 

title3=strcat(title3  ,position,*cm’); 

title4-  Iteration  =  '; 

title4=strcat(title4,iteration); 

title_data=char(  {title  1  ,title2,title3,title4  }); 

title(title_data); 

yIabel(Displacement*) 

xlabelCFrequency  (Hz)*) 

legend(TundVHarm  I'/Harm  2',’Harm  3VHarm  4','Harm  5') 
hold  off 

orient  landscape 
figure(2) 

plot(freq20,abs(displ_freq20_fund),’-’) 
hold  on 

plot(freq20,abs(displ  _freq20_harm  1 
plot(freq20,abs(displ_freq20_harm2),'-.') 
plot(freq20,abs(displ_freq20_hann3 
plot(freq20,abs(displ_freq20_hann4),'-') 
plot(freq20,abs(displ_freq20_harm5),'“') 
title  1- Displacement  vs  Frequency*; 
title2-  Amplitude  =  '; 
titles-  Position  = 

title2=strcat(title2,ampIitude,’Volts'); 
title3=strcat(title3,position,’cm’); 
title4='  Iteration  = '; 
title4=strcat(title4,iteration); 
title_data=char({title  1  ,title2,title3  ,title4}); 
title(title_data); 
ylabelCDisplacement') 
xlabelCFrequency  (Hz)*) 

legen<iCFund','Harm  r,’Harm  2VHarm  3’,’Harm  4’, 'Harm  5*) 
hold  off 
orient  landscape 
elseif  position=- 40.00' 
figure(l) 

semilogy(freq40,abs(displ_freq40_fund),'-') 
hold  on 

semilogy(freq40,abs(displ_freq40_harm  1 

semilogy(freq40,abs(displ_freq40_harm2),'-.') 

semilogy(freq40,abs(displ_freq40_harm3),*:') 

semilogy(freq40,abs(displ_freq40_harm4),’-') 

semilogy(freq40,abs(dtspl_freq40_harm5),'-') 

titlel- Displacement  vs  Frequency*; 

title2='  Amplitude  =  '; 

titles-  Position-'; 

title2=strcat(title2,amplitude,*Volts'); 

title3=strcat(title3  ,position,'cm'); 

title4='  Iteration  = '; 

title4=strcat(title4,iteration); 

titledata=char(  {titlel  ,title2,title3  ,title4}  ); 

title(titledata); 
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ylabelCDisplacement') 
xlabelf  Frequency  (Hz)*) 

legendCFund’,’Hann  I '/Harm  2VHarm  3','Harm  4VHarm  S') 
hold  off 
orient  landscape 
figure(2) 

plot(freq40,abs(displ_freq40jund),’-') 
hold  on 

plot(freq40,abs(displ_freq40_hannl),'‘-') 
plot(freq40,abs(displ_freq40__harm2);-.') 
plot(freq40,abs(displ_freq40_harm3),':') 
plot(freq40,abs(displ_fireq40_harm4),'-') 
plot(freq40,abs(displ_fireq40_hann5),'— *) 
title  l='Displacement  vs  Frequency*; 
title2='  Air^litude  = '; 
titles-  Position  =  ’; 

title2=strcat(title2,aniplitude,’Volts'); 
title3=strcat(title3,position,'cm'); 
title4=’  Iteration  =  '; 
title4=strcat(title4, iteration); 
title_dala=char(  {title  1  ,title2,title3  ,title4  }); 
title^itledata); 
ylabel(T)isplacement') 
xIabel(Trequency  (Hz)*) 

Iegend(Tund’,’Harm  IVHaim  2’,'Hami  3',’Hann  4’,’Harm  S') 
hold  off 
orient  landscape 
elseif  position-- 80.00' 
figure(l) 

semilogy(freq80,abs(displ_freq80_fUnd),'-’) 

hold  on 

semilogy(freq80,abs(displ_freq80_hann  1 ),'—') 
seniilogy(freq80,abs(displ_freq80_ham[i2),  ) 
semilogy(freq80,abs(displ_fi‘eq80_harm3  )/ :') 
semilogy(freq80,abs(displ_freq80__hami4),'-*) 
semilogy(freq80,abs(displ_freq80_hamiS),  —') 

title  1- Displacement  vs  Frequency*; 
title2='  Amplitude  = '; 
titles-  Position  =  '; 

title2-strcat(title2,amplitude, 'Volts'); 
title3=strcat(title3,position,'cm'); 
title4='  Iteration  =  '; 
title4=strcat(title4,iteration); 
title_data=char({title  1  ,title2,title3,title4 } ); 
title(title_data); 
ylabel(Displacement') 
xlabel(Trequency  (Hz)') 

legend(Tund',’Harm  l’,'Harm  2’,'Hann  3','Hann  4','Harm  S') 
hold  off 
orient  landscape 
figure(2) 

plot(freq80,abs(displ_freq80_fund),'-’) 

hold  on 

plot(freq80,abs(displ_freq80_harm  1  ),'-*) 

plot(jfreq80,abs(displ_freq80_harm2),'-.') 

plot(freq80,abs(displ_freq80_hann3),':') 

plot(freq80,abs(displ_freq80_harm4),'-*) 

plot(freq80,abs(displ_freq80_harm  S  ),'—') 

title  I  ='I)isplacement  vs  Frequency  ; 

title2=’  Amplitude  =  '; 

titles-  Position  ==’; 

title2=strcat(title2,amplitude,' V  olts'); 

title3=strcat(title3,position,'cm'); 

title4-  Iteration  =  ’; 

title4=strcat(title4,itcration); 

title  data-char( (title  1  ,title2,title3,title4 } ); 
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title(title_data); 
ylabelfDisplacemenO 
xlabelfFrequency  (Hz)*) 

legendCFund'/Harm  IVHarm  2';Hami  3VHarm  4VHarm  5') 
hold  off 
orient  landscape 
elseif  position— ’160.0’ 
figure(l) 

seniilogy(freq  1 60,abs(displ_freq  1 60  fund),'-') 
hold  on 

semilogy(freql  60,abs(displ_freql  60_harm  1 ),'—') 

semilogy(freq  1 60,abs(displ_freq  1 60_harm2),'-. ) 

semilogy(freql  60,abs(displ_freql  60_hann3),':') 

semilogy{freq  1 60,abs(displ_freq  1 60_harm4),'-^ 

semilogy(freql  60,abs(displ_freq  1 60harm5),'— *) 

title  l='Displacement  vs  Frequency*; 

title2=^'  Amplitude  = '; 

title3-  Position  = 

title2=strcat(title2, amplitude, 'Volts'); 

title3=strcat(title3,position,'cm'); 

title4='  Iteration  =  '; 

title4=strcat(title4,iteration); 

title_data=char({titlel,title2,title3,title4}); 

title(title_data); 

ylabelCDisplacemenf) 

xlabelCFrequency  (Hz)') 

legend(Tund','Hann  r,'Hann  2','Harm  3', 'Harm  4','Hann  5') 
hold  off 
orient  landscape 
figure(2) 

plot(freq  1 60,abs(displ_freq  1 60_fund),'-') 
hold  on 

pIot(freq  1 60,abs(displ_freq  1 60_harm  1 ),'-') 
plot(freql  60,abs(displ_freq  1 60_hami2),'-.') 
plot(freql60,abs(displ_freql60_harm3  ),':') 
piot(freq  1 60,abs(displ_freq  1 60_harm4),'-') 
plot(freq  1 60,abs(displ_freq  1 60_harm5  ),'—') 
title  1- Displacement  vs  Frequency'; 
title2-  Amplitude  =  '; 
title3-  Position  =  '; 
title2=strcat(title2,amplitude,'Volts'); 
title3=strcat(title3  ,position,'cm'); 
title4-  Iteration  =  '; 
title4=strcat(title4,iteration); 
title_data=char(  {title  1  ,title2,titl  e3  ,title4  }  ); 
title(title_data); 
ylabelCDisplacemenf) 
xlabelCFrequency  (Hz)') 

legendCFund','Harm  i’,’Harm  2',’Harm  3',’Harm  4','Harm  5*) 
hold  off 
orient  landscape 


end 


APPENDIX  E 


LAB  VIEW  CODE 

Lab  VIEW  programming  was  done  by  Dr.  Gregg  D.  Larson,  a  research  engineer 
assigned  to  the  investigation  of  acousto-electromagnetic  mine  detection  at  the  Georgia 
Institute  of  Technology.  The  software  evolved  from  existing  code  and  continued  to 
evolve  over  the  year  that  data  was  collected  for  this  thesis.  This  appendix  includes  only  a 
couple  of  the  numerous  subroutines  utilized  for  taking  data  in  Experiment  2. 

The  program  found  on  pages  149  and  150  was  the  major  program  collecting  the 
radar  data.  This  code  drove  several  subroutines  and  was  driven  itself  by  other  programs. 
The  overall  controlling  code  is  found  in  the  second  program  (page  151).  From  this  code, 
all  of  the  data  for  Experiment  2  was  collected.  This  included  the  four-accelerometer 
measurements  and  the  alternating  two-accelerometer/radar  measurements  for  both  the 
frequency  response  tests  and  the  saturation  curve  tests. 
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